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PREFACE 


The  High-Speed  Research  Program  and  NASA  Langley  Research  Center  sponsored  the  NASA 
High-Speed  Research  Program  Aerodynamic  Performance  Workshop  on  February  25-28, 
1997.  The  workshop  was  designed  to  bring  together  NASA  and  industry  High-Speed  Civil 
Transport  (HSCT)  Aerodynamic  Performance  technology  development  participants  in  areas  of 
Configuration  Aerodynamics  (transonic  and  supersonic  cruise  drag  prediction  and  minimiza- 
tion), High-Lift,  Flight  Controls,  Supersonic  Laminar  Flow  Control,  and  Sonic  Boom  Predic- 
tion. The  workshop  objectives  were  to  (1)  report  the  progress  and  status  of  HSCT  aerodynamic 
performance  technology  development;  (2)  disseminate  this  technology  within  the  appropriate 
technical  communities;  and  (3)  promote  synergy  among  the  scientist  and  engineers  working 
HSCT  aerodynamics.  In  particular,  single-  and  multi-point  optimized  HSCT  configurations  and 
HSCT  high-lift  system  performance  predictions  were  presented  along  with  executive  summa- 
rizes for  all  the  Aerodynamic  Performance  technology  areas. 

The  workshop  was  organized  in  three  sessions  as  follows: 

Session  I Plenary  Session 

Session  II  Independent  Session 

Session  III  Executive  Summaries 

The  proceedings  are  published  in  two  volumes: 

Volume  I,  Parts  1 and  2 Configuration  Aerodynamics 

Volume  II  High  Lift 

Conference  Chairmen:  Daniel  G.  Baize  and  Robert  L.  Calloway 

NASA  Langley  Research  Center 
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Distributed  Computing  for  CFD 
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Grid  Splitting  Strategy  for  CFL3Dhp 
TCA  wing/body  Navier-Stokes  grid,  1.5  Million  Grid  points 
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Load  Balancing 
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Load  Distribution  for  the  TCA  W/B/N/D  Euler  Grid 
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CFL3Dhp  Applications  on  the  IBM  SP-2 
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Residual  History  Comparison  (Cray  C-90  vs.  IBM  SP-2) 

CFL3D  (C-90)  and  CFL3Dhp  (SP-2)  Navier-Stokes  (B-L),  3 level  multigrid,  12  or  16  blocks 
TCA  W/B,  1 .5  million  grid  points,  M„=2.4,  Re=4x106/ft 
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Drag  Convergence  Comparison  (Cray  C-90  vs.  IBM  SP-2) 

CFL3D  (C-90)  and  CFL3Dhp  (SP-2),  Navier-Stokes  (B-L),  3 level  multigrid,  12  or  16  blocks 
TCA  W/B,  1.5  million  grid  points,  M^=2.4,  Re=4x106/ft 
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Pressure  Drag  Coefficients  (CDp)  for  the  TCA  W/B/E  Configuration 

CFL3D  Euler,  M =2.4 
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CFL3D  Navier-Stokes  Grid  Refinement  Study 
TCA  W/B  Configuration,  M„=2.4,  Re=6.36x106,  <x=3.518 
Solution  Performed  on  IBM  SP-2 
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Predicted  Drag  Polar  for  the  TCA  Wing/Body  Configuration 

CFL3Dhp  Navier-Stokes  solution,  1^=2. 1,  Re=4x106/ft 
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onclusions  and  Future  Work 
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Full  Configuration  Force  and  Moment  Calculations  using 
Multiblock  CFL3D  on  HSCT  Configurations 

Grant  L.  Martin  and  Robert  P.  Narducci 
McDonnell  Douglas  Corporation 
Long  Beach,  California  90807-5309 

During  the  past  year,  the  McDonnell  Douglas  Corporation  (MDC)  has 
made  large  strides  in  Computational  Fluid  Dynamics  (CFD)  analysis  of 
increasingly  complex  HSCT  configurations  using  both  serial  and  parallel 
computational  platforms.  While  tools  for  grid  generation  and  analysis  on 
serial  computers  have  remained  relatively  unchanged,  a new  gridding 
strategy  has  been  employed  to  obtain  Navier-Stokes  analyses  of  HSCT 
configurations  which  include  the  wing,  body,  nacelles,  diverters,  and 
empennage.  Additionally,  with  the  promising  efficiency  of  parallel 
machines,  MDC  has  contributed  to  the  development  of  CFL3Dhp,  a parallel 
version  of  CFL3D  for  the  IBM  SP-2. 

Presented  herein  are  full  configuration  Euler  and  Navier-Stokes  solutions 
obtained  using  CFL3D  on  the  NAS  C-90  and  IBM  SP-2.  With  the 
objectives  of  validating  CFL3D  for  supersonic  cruise  calculations  on 
several  platforms,  CFD  results  for  the  Reference  H and  Technology 
Concept  Airplane  (TCA)  configurations  are  presented  in  a build-up 
fashion.  The  build-up  fashion  entails  analyzing  the  simplest  of 
configuration  first,  the  wing/body  (W/B)  followed  by  the  additional 
complexity  of  the  empennage  (W/B/E),  then  nacelles  and  diverters 
(W/B/N/D),  and  finally  the  entire  configuration  (W/B/N/D/E).  A thorough 
build-up  has  been  performed  on  the  Reference  H configuration,  while  the 
TCA  build-up  work  is  still  in  progress.  To  assist  in  the  validation,  a 
number  of  comparisons  are  made  to  available  experimental  data  from  the 
NASA  Langley  Unitary  Plan  Wind  Tunnel  (UPWT); 
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Configuration  Aerodynamics  Workshop 
February  25-28,  1997 

NASA  Langley  Research  Center,  Hampton,  VA 


In  this  report,  we  investigate  flow  over  full  configuration  (W/B/N/D/E)  HSCT  configurations  using 
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Summary  and  conclusions 
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Objectives 
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Flow  Analysis  and  Platforms 
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Flow  Analysis  and  Platforms 
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Clutter  charts  used  for  skin-friction  drag  estimates  to 
complement  Euler  drag  predictions 
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Reference  H Build-up 
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- Euler  solutions  for  W/B,  W/B/N/D,  W/B/E,  & W/B/N/D/E 

- Navier-Stokes  solutions  for  W/B  & W/B/E 
Empennage  and  nacelle/diverter  drag  increments 


1.675% -Scale  HSR  Reference  H Modular  Controls  Model 

Extended  Aft  Body  with  Tails 

NASA  Langley  Unitary  Plan  Wind  Tunnel,  Section  #1 
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1.675%  Scale  HSR  Reference  H Modular  Controls  Model 

Extended  Aft  Body  with  Tails 

NASA  Langley  Unitary  Plan  Wind  Tunnel,  Section  #1 
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W/B  and  W/B/E  Supersonic  Grids 
1.675%  Reference  H Modular  Controls  Model 
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W/B/N/D  and  W/B/N/DE  Supersonic  Euler  Grids 
1.675%  Reference  H Modular  Controls  Model 
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Ref.  H Full  Configuration  Supersonic  Euler  Grid 
.675%  Modular  Wind-Tunnel  Model,  W/B/N/D/E  Configuration 
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Ref.  H Full  Configuration  Supersonic  Euler  Grid 
1 .675%  Modular  Wind-Tunnel  Model,  W/B/N/D/E  Configuration 


Comparison  of  Predicted  and  Experimental  Lift  Curves 
1.675%  Reference  H Modular  Controls  Model 

CFL3D  & Langley  UPWT  1812,  M^.l,  Re=3xlQ6/ft 
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Comparison  of  Predicted  and  Experimental  Lift  Curves 

1.675%  Reference  H Modular  Controls  Model 

CFL3D  & Langley  UPWT  1812, 1^=2. 1,  Re=3x106/ft 
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Angle-of-Attack,  a(°) 


Comparison  of  Predicted  and  Experimental  Drag  Polars 
1 .675%  Reference  H Modular  Controls  Model 
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Comparison  of  Predicted  and  Experimental  Drag  Polars 
1.675%  Reference  H Modular  Controls  Model 

CFL3D  & Langley  UPWT  1812, 1^=2.!,  Re=3xl06/ft 
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Comparison  of  Predicted  and  Experimental  Lift-to-Drag  Ratio 
1.675%  Reference  H Modular  Controls  Model 

CFL3D  & Langley  UPWT  1812,  Moo=2.1,  Re=3xl06/ft 
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Comparison  of  Predicted  and  Experimental  Lift-to-Drag  Ratio 

1 .675%  Reference  H Modular  Controls  Model 

CFL3D  & Langley  UPWT  1812,  M.=2.1 , Re=3x106/ft 
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Comparison  of  Predicted  and  Experimental  Pitching  Moments 
1.675%  Reference  H Modular  Controls  Model 

CFL3D  & Langley  UPWT  1812,  Moo=2.1,  Re=3xl06/ft 


D 

Dh 

G 

D 

D 

X 

73 

(D 

G 

D 

oo 

CD 


G 

0) 

s 

o 

S 


o 

£ 


G 

G 

73 

G 

G 

CD 

£ 

» i— H 

Uh 

d 


D 

X 


D 

X 


bJQ 

.S  O 

X ~ 
o 73 

.ti  a) 


CD  CD 

* D 

X 


D 


box 

G ?-> 


X 

O 


G 


bJQ 

G 

Ui 

73 


a.S 


G 

CD 

£ 

o 

a 

C/3 

G 

O 

G 


D 

bJO 

S-i 

G 


D 

bO 


^ O 

3 ° 

t5  W 

CD  X 

D*  LJ 

cdg: 


£ CQ 

is 

o 

cj  73 
G G 
- G 


G 

O 


73 

CD 


3 Q 
JS  Uh 

G (J 

o 


o 
o 

Cd  bO 

G.S 


(D 


X 

o 


Q 

Pd 

U 

Vh 

o 


Q 

Z 

OQ 


CD' 

c/5 

CD 

Q 

73 

(D 


D 

X 


CD 


CD 

£ p 


G UJ 
O 

X CQ 


r <P 


G 

<D 

D 

X) 

C/5 

G 

X 


CQ 


G 

C 

CD 

<D 

X 

C/5 

G 

X 


o 

G 


73 

G 

G 

G 

D 

6 

• «— 4 

U< 

CD 

CD 

X 

(D 

G 

D 

<D 

£ 

D 

X 

D 

O 

G 

D 

D-t 

D 


G 

C/5 

G 


5— ! 

,o 


>> 

* % 

i 

G X 

3 

.2  > 

a 

X Uh 

° 

U 

a o 

o 

* ^ D— i 

73 

C/3 

D G 

• 

D-i  O 

X 

CD  X 

^ w 

bD 

G 

Ui 

Q 


G 

D 


D 

X 


C/5 

D 

C/5 

G 

D 

u 

O 

c 


D 

3 

PQ 


CJ 

D 

u* 

O 

a 


D 

D 

u 

bJQ 

G 

73 

O 

c 

bJQ 


D 

73 

3 

w 

• 

C 

bJD 

G 


68 


Comparison  of  Predicted  and  Experimental  Pitching  Moments 
1 .675%  Reference  H Modular  Controls  Model 

CFL3D  & Langley  UPWT  1812,  M_=2.1 , Re=3x106/ft 
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Comparison  of  Predicted  and  Experimental  Drag  Polars 
1.675%  Reference  H Modular  Controls  Model 

CFL3D  & Langley  UPWT  1812,  M^.l,  Re=3xlQ6/ft 
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Comparison  of  Predicted  and  Experimental  Drag  Polars 

1 .675%  Reference  H Modular  Controls  Model 

CFL3D  & Langley  UPWT  1812,  Mto=2.1,  Re=3x106/ft 
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Comparison  of  Predicted  and  Experimental  Drag  Polars 
2.7%  Reference  H,  W/B  Model 

CFL3D,  TLNS3D,  & Ames  9x7  Tunnel,  Moo=2.4,  Rec=7xlQ6 
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Comparison  of  N-S  CFD  and  Experimental  Results:  CL  vs.  C 
2.7%  Ref.  H Wing/Body,  M j=2.4,  Rec=7x1 0® 
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Nacelle/Diverter  and  Empennage  Effect  on  Drag 
1.675%  Reference  H Modular  Controls  Model 

CFL3D  & Langley  UPWT  1812,  Moo=2.1,  at  CL=0.12,  Re=3xl()6/ft 
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Component  Drag  Increments 
1.675%  Reference  H Modular  Controls  Model 
CFL3D  & Langley  UPWT  1812,  M =2.1,  C,=0.12,  Re=3x10 
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ExpEulN-S  ExpEulN-S  ExpEulN-S 

W/B/E  W/B/N/D  W/B/N/D/E 


Empennage  Effect  on  Pitching  Moment 
1.675%  Reference  H Modular  Controls  Model 

CFL3D  & Langley  UPWT  1812,  Moo=2.1,  at  CL=0.12,  Re=3xl()6/ft 
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Component  Pitching  Moment  Increments 
1.675%  Reference  H Modular  Controls  Model 
CFL3D  & Langley  UPWT  1812,  M =2.1,  C=0.12,  Re=3x106/ft 


00 

CD 

LT) 

CO 

CM 

T — 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

d 

o 

86euu0diu3  o\  eng  ^qv 


77 


W/B-»W/B/E  W/B/N/D— >W/B/N/D/E 


HSR  Technology  Concept  Airplane  (TCA) 
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TCA  Build-Up 
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closed  aftbody 


TCA  W/B/E  Grid  Structure,  iH=0° 
12  Zones,  3.1  Million  Grid  Points 
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HSR  Technology  Concept  Airplane  (TCA)  W/B/E,  i 
12  Zones,  3.1  Million  Grid  Points 


Comparison  of  Predicted  and  Experimental  Lift  Curves 
1.675%  TCA  Model,  W/B  Configuration 

CFL3D  & Langley  UPWT  1671,  M^.4,  Re=4xl06/ft 
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Comparison  of  Predicted  and  Experimental  Lift  Curves 
TCA  Wing/Body  Configuration 

Langley  Test  1671,  UPWT-2,  Mto=2.4,  Re=4x106/ft 
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Angle-of-Attack,  a(°) 


Comparison  of  Predicted  and  Experimental  Drag  Polars 
1.675%  TCA  Model,  W/B  Configuration 

CFL3D  & Langley  UPWT  1671,  M00=2.4,  Re=4xl06/ft 
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Comparison  of  Predicted  and  Experimental  Drag  Polars 
TCA  Wing/Body  Configuration 

Langley  Test  1671,  UPWT-2,  Moo=2.4,  Re=4x106/ft,  Trip  Drag  Removed 
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Comparison  of  Predicted  and  Experimental  L/D  Ratios 
1.675%  TCA  Model,  W/B  Configuration 

CFL3D  & Langley  UPWT  1671,  M00=2.4,  Re=4xl06/ft 
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Comparison  of  Predicted  and  Experimental  Lift-to-Drag  Ratios 

TCA  Wing/Body  Configuration 

Langley  Test  1671,  UPWT-2,  M_=2.4,  Re=4x106/ft,  Trip  Drag  Removed 
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Comparison  of  Predicted  and  Experimental  Pitching  Moments 
1.675%  TCA  Model,  W/B  Configuration 

CFL3D  & Langley  UPWT  1671,  Moo=2.4,  Re=4xl06/ft 


2 

H C3 

b| 

0 <5 

Oh 

• k>t 

1 S 

o <D 


c 6 
o U 
u ro 

'u 

x>  §> 

to  O 


£ bO 
O oc 

x a 


> .a 

b t3 

^ U. 

Cl 

C c/: 
<u 

e 

o ° 

£ 00 
M i 


-C 

a z 

'5.-0 

a 

CQ  cd 


.2  W 

H Q 


cn 

»J  03 

ft  ts 

U T3 


88 


Comparison  of  Predicted  and  Experimental  Pitching  Moments 

TCA  Wing/Body  Configuration 

Langley  Test  1671,  UPWT-2,  M_=2.4,  Re=4x10“/ft 
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Comparison  of  Predicted  and  Experimental  Lift  Curves 
1.675%  TCA  Model,  W/B/N/D  Configuration 

CFL3D  & Langley  UPWT  1671,  M„=2.4,  Re=4xlO<Vft 
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Comparison  of  Predicted  and  Experimental  Lift  Curves 
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Angle-of-Attack,  a(°) 


Comparison  of  Predicted  and  Experimental  Drag  Polars 
1.675%  TCA  Model,  W/B/N/D  Configuration 

CFL3D  & Langley  UPWT  1671,  Moo=2.4,  Re=4xl06/ft 
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Comparison  of  Predicted  and  Experimental  Drag  Polars 
TCA  Wing/Body/Nacelle/Diverter  Configuration 

Langley  Test  1671 , UPWT-2,  M„=2.4,  Re=4x106/ft,  Trip  Drag  Removed 
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Comparison  of  Predicted  and  Experimental  L/D  Curves 
1.675%  TCA  Model,  W/B/N/D  Configuration 

CFL3D  & Langley  UPWT  1671,  M^.4,  Re=4xlQ6/ft 
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Comparison  of  Predicted  and  Experimental  Lift-to-Drag  Ratios 
TCA  Wing/Body/Nacelle/Diverter  Configuration 

Langley  Test  1671,  UPWT-2,  M_=2.4,  Re=4x106/ft,  Trip  Drag  Removed 
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Comparison  of  Predicted  and  Experimental  Pitching  Moments 
1.675%  TCA  Model,  W/B/N/D  Configuration 

CFL3D  & Langley  UPWT  1671,  M^.4,  Re=4xlQ6/ft 
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Comparison  of  Predicted  and  Experimental  Pitching  Moments 

TCA  Wing/Body/Nacelle/Diverter  Configuration 

Langley  Test  1671 , UPWT-2,  Mw=2.4,  Re=4x106/ft 
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Moo  2.4,  CL=0. 1 , Re=4x  1 Q6/ft. 
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Lift  Coefficients  for  the  TCA  W/B/E  Configuration 

CFL3D  Euler,  M^=2.4 
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Lift  C 
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Pressure  Drag  Coefficients  (CDp)  for  the  TCA  W/B/E  Configuration 

CFL3D  Euler,  Moo=2.4 
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Pitching  Moments  for  the  TCA  W/B/E  Configuration 

CFL3D  Euler,  M«,=2.4 
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Pitching  Moment  for  the  TCA  Wing/Body/Empennage  Configuration 

CFL3D  Euler,  M =2.4 
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Pitching  Moment,  C 


Pressure  Contours  on  the  TCA  Empennage 
CFL3D  Euler,  Moo=2.4 
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Pressure  Contours  on  HSR  Technology  Concept  Airplane  (TCA)  W/B/E 
CFL3D  Euler,  12  Zones,  3.1  Million  Grid  Points 


Pressure  Contours  on  the  TCA  Empennage 
CFL3D  Euler,  M„=2.4 
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Computer  Resources  Used  for  Calculations  on  Various  Configurations 

at  Supersonic  Conditions 
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Computer  Resources  Used  for  Calculations  on  Various  Configurations 

at  Supersonic  Conditions 


Memory 

Req. 

(Mw) 

40 

68 

58 

85 

O) 

104 

85 

212 

06 

CPU 

Time 

(hrs) 

1.5 

o 

CM 

LO 

co 

3.5 

CO 

r 3.5 

Platform 

C-90 

J-90 

C-90 

SP-2 

C-90 

C-90 

C-90 

C-90 

o 

G> 

1 

O 

Grid  Size 
(Millions) 

60 

1.5 

1.5 

CO 

2.2 

CM 

CM 

CM 

o 

CD 

2.6 

Turb. 

Model 

N/A 

_i 

i 

m 

V/N 

V/N 

B-L 

B-B 

V/N 

B-L 

V/N 

Solution 

Type 

Euler 

N-S 

Euler 

Euler 

CO 

i 

N-S 

Euler 

N-S 

Euler 

Config. 

W/B 

W/B 

W/B/E 

W/B/E 

W/B/E 

W/B/E 

W/B/N/D 

W/B/N/D 

W/B/N/D/E 

in 


The  flow  over  two  full  configuration  (W/B/N/D/E)  HSCT  configurations  were  analyzed  using 
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Summary  & Conclusions 
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and  wind-tunnel  for  TCA  W/B  and  W/B/N/D 
TCA  full  configuration  analysis  is  underway 


Supersonic  Cruise  Point  Design  Optimization  of  TCA 

Eric  R.  Unger 
Robert  P.  Narducci 
James  O.  Hager 
Geojoe  Kuruvila 
Peter  M.  Hartwich 
Shreekant  Agrawal 

McDonnell  Douglas  Corporation 
Long  Beach,  California  90807-5309 

Since  July  of  1996,  McDonnell  Douglas  (along  with  other  teams  from 
NASA  Ames  and  Boeing  Commercial  Aircraft  Group),  has  been  working  on 
a second  series  of  optimizations  for  the  TCA  configuration.  The  approach 
used  at  MDC  was  conservative  in  terms  of  acceptable  geometric  qualities 
that  were  allowed  to  appear  in  the  final  Cycle  2 design.  The  hope  was  that 
any  final  outcome  would  be  more  robust  and  raise  the  least  amount  of 
uncertainties  from  other  technology  disciplines.  The  downside  of  this 
approach  was  the  inability  to  fully  maximize  the  possible  L/D  gains  within 
the  given  time  and  within  these  strict  geometric  guidelines. 

This  paper  presents  an  overview  of  MDC’s  final  Cycle  2 configuration. 
First  a brief  introduction  and  highlights  of  the  new  design  are  given  along 
with' some  geometric  details.  Second,  a look  at  the  configuration’s  overall 
performance  and  pressure  field  details  will  be  given.  Next,  some  details  of 
the  design  constraints  that  were  used  during  optimization  will  be  described. 
And  finally,  the  paper  will  close  with  a summary  of  the  Cycle  2 
configuration  and  a look  ahead  to  the  immediate  future. 
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Supersonic  Cruise  Point 
Design  Optimization  of  TCA 
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Listed  below  are  the  contents  of  the  paper  to  be  discussed. 
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TCA  Nonlinear  Point  Design 
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- Three  very  different  final  outcomes 
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TCA  Cycle  2 Highlights 
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7.4%  increase  in  fuel  volume 
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McDonnell  Douglas  Wing/Body 
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McDonnell  Douglas  Wing/Body 


Waterline  Contours  for  the  Baseline  TCA  Configuration 


g 

.2? 

O 

a 

CD 

Uh 

00 

CD 

*00 

o 

CD 

OJO 

<D 

’t— • 

4= 

G 

XJ 

<N 

CD 

00 

£ 

g 

cd 

4= 

a 

jD 

=3 

.2 

7j 

00 

• 

%— * 
cd 

CD 

G 

5-h 

U 

"G 

g 

C/0 

CD 

a 

.cd 

C+H 

U 

=3 

0£> 

tfl 

o 

X) 

cd 

O 

G 

<D 

G 

00 

o 

a 

CD 

*-H 

G 

cd 

cd 

G 

CD 

00 

cd 

• k-H 
£ 

o 

00 

d 

i-H 

CD 

43 

PQ 

J-h 

00 

op 

» 

< 

u 

CD 

£ 

o 

<d 

+-» 

*00 

<D 

"G 

+-» 

cd 

44 

4-» 

Jg  c 

■5  cS 

J— I ^ 

O g 

4h  Oh 
00  2^ 

rr  m 


G C 

. r— ( ^ 


2 a? 

.S  c 


o c 

-C  a> 


<U  Oh 

'O  c 


3 3 
3 § 
S 8 
g B 

03  oi) 

j3  3 

o 44 

44  O 

g o 


G C4 

c 2 

o 

4h  b 
G O 
< c3  o 
f • ▼— < 

Qh  +s 


129 


Waterline  Contours  for  the  TCA  Cycle  2 Configuration 
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Illustrated  below  is  a comparison  of  the  Baseline  and  Cycle  2 twist  distributions.  The  Cycle 
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Wing  Twist  Distribution 

Baseline  and  Cycle  2 TCA  Confiqurations 
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Camber  Distribution 

Baseline  and  Cycle  2 TCA  Configurations 
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selage  Station  (in.) 


Baseline  and  Cycle  2 Maximum 
Thickness-to-Chord  Ratio  Distributions 
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Maximum  Thickness-to-Chord  Ratio  Distribution 

Baseline  and  Cycle  2 TCA  Configurations 
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Span  Station  (in) 


TCA  Baseline  and  Cycle  2 Fuselage 
Cross-Sectional  Area  Distributions 
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TCA  Fuselage  Cross-Sectional  Area  Distributions 
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1000  2000  3000  4000 

Fuselage  Station  (in.) 


TCA  Baseline  and  Cycle  2 Wing/Body 
Cross-Sectional  Area  Distributions 
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TCA  Cross-Sectional  Area  Distributions 

(Wing/Body  transition  is  based  on  the  grid) 
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1000  2000  3000  4000 

Fuselage  Station  (in.) 


TCA  Baseline  and  Cycle  2 W/B/N/D  Lift  Comparisons 
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Lift  for  the  TCA  W/B/N/D  Configurations 

CFL3D,  Navier-Stokes  (Baldwin-Lomax),  Mj=2.4,  Rec=6.36x10 
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Angle-of-Attack,  a (°) 
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TCA  Baseline  and  Cycle  2 W/B/N/D 
Pitching  Moment  Coefficient  Comparisons 
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Pitching  Moment  for  the  TCA  W/B/N/D  Configurations 

CFL3D,  N-S  (B-L),  Mm=2.4,  Rec=6.36x1 06,  Moment  about  50%  MAC 
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TCA  Baseline  W/B/N/D  Pressure  Distributions 
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TCA  Cycle  2 W/B/N/D  Pressure  Distributions 

CFL3D,  N-S  (B-L),  M „=2.4,  CL=0.10  (cc=3.45°),  Rec=6.36x10 
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Component  Drag  for  the  Baseline  and  Cycle  2 
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TCA  Baseline  and  Cycle  2 Pressure  Coefficient 
and  Airfoil  Cut  Comparisons 
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TCA  Baseline  and  Optimized  Configurations 

Span  Station  156  (19.9%  Semi-span) 

CFL3D  Navier-Stokes  (Baldwin-Lomax),  1.675%  W/B/N/D  Configuration,  Moo=2.4,  a=3.518°,  Re  =6.36x10' 
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Chord  Fraction,  x/c 
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TCA  Baseline  and  Optimized  Configurations 

Span  Station  256  (32.6%  Semi-span) 

CFL3D  Navier-Stokes  (Baldwin-Lomax),  1.675%  W/B/N/D  Configuration,  M„=2.4,  o=3.518°,  Re 


Chord  Fraction,  x/c 


TCA  Baseline  and  Cycle  2 Pressure  Coefficient 
and  Airfoil  Cut  Comparisons  (cont.) 

This  slide  shows  the  pressure  and  airfoil  variations  at  a wing  semi-span  of  41.3%.  Once 
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TCA  Baseline  and  Optimized  Configurations 

Span  Station  325  (41 .3%  Semi-span) 

CFL3D  Navier-Stokes  (Baldwin-Lomax),  1 .675%  W/B/N/D  Configuration,  M„=2.4,  a=3.518°,  Rec=6.36x10 
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Chord  Fraction,  x/c 


TCA  Baseline  and  Cycle  2 Pressure  Coefficient 
and  Airfoil  Cut  Comparisons  (cont.) 

This  slide  shows  the  pressure  and  airfoil  variations  at  a wing  semi-span  of  67.9%  an 
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TCA  Baseline  and  Optimized  Configurations 

Span  Station  535  (67.9%  Semi-span) 


Chord  Fraction,  x/c 
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Wing  Constraints 
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see  figures 


Front  Spar  (Wl)  Sensitivity  Study 
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Front  Spar  Thickness  (W1 ) Sensitivity  Study 
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TCA  Cycle  2 Landing-Gear  Bay  Layout 

(Gear  Bay  Shown  in  Red) 


This  figure  shows  a graphical  representation  of  the  fuel  volume  layout  on  the  Cycle  2 
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Body  Constraints 
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Fuselage  Constraints  (cont.) 
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Body  Constraints 
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Fuselage  Constraints  (cont.) 
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Body  Constraints 
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Fuselage  Station  (in.) 
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TCA  Fuse 
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Straight  streamlines 
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Improvements  to  the  MDC 
Nonlinear  Aerodynamic  Design  Tools 

James  O.  Hager,  Peter  M.  Hartwich,  Eric  R.  Unger, 
Geojoe  Kuruvila,  Robert  P.  Narducci,  and  Shreekant  Agrawal 
McDonnell  Douglas  Corporation 
Long  Beach,  California  90807-5309 


Nonlinear  aerodynamic  optimization  is  considered  a key  technology 
required  to  develop  a High  Speed  Civil  Transport  (HSCT).  Within  the 
High  Speed  Research  (HSR)  program,  McDonnell  Douglas  is  developing 
nonlinear  optimization  tools  to  be  able  to  support  the  launch  of  an  HSCT 
program  at  the  end  of  HSR  II.  This  paper  presents  recent  improvements  to 
the  tools. 

The  first  set  of  improvements  were  made  to  be  able  to  optimize  the 
Technology  Concept  Aircraft  (TCA).  The  TCA  presented  some  grid- 
generation  issues  because  it  is  a true  low-wing  configuration.  In  addition, 
several  constraints  were  required  to  maintain  a realistic  design. 

Second,  the  geometry  modeling  capability  was  improved  to  move  toward 
full-configuration  modeling.  Empennage  effects  have  been  modeled,  and 
wing/body/flaps  configurations  can  be  modeled.  Efforts  were  also  made  to 
produce  and  improve  tools  required  for  integrated  wing/body/nacelle/ 
diverter  modeling. 

Third,  alternate  gradient  evaluation  techniques  are  being  examined  to 
replace  the  finite-difference  calculations  currently  being  used.  ADIFOR 
was  applied  to  CFL3D  and  demonstrated  for  a 100+  design-variable 
problem.  Also,  an  adjoint  module  is  being  created  for  TLNS3D. 

Finally,  a transition  is  being  made  to  a modular  design  environment  to 
facilitate  improvements  and  the  addition  of  new  codes. 
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Improvements  to  the  MDC  Nonlinear 
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NASA  Langley,  25-28  February  1 997 
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patched-grid  capability  will  allow  flexibility 
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Improved  Design  Capability: 

Design  Variables 

Two  types  of  design  variables  were  modified.  The  wing  sinusoidal  shape  functions 
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Improved  Design  Capability: 
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Improved  Geometry  Capability: 
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Flap  Optimization  Grid 
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Improved  Geometry  Capability: 

Integrated  W/B/N/D  Modeling 

This  chart  outlines  the  procedure  that  will  be  used  to  perform  integrated  W/B/N/D 
odeling.  The  next  several  charts  will  fill  in  some  of  the  details. 
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Baseline  Nacelle/Diverter  Reintegrated  to  a New  Wing  Geometry 

This  slide  gives  a demonstration  of  the  nacelle/ diverter  (N/D)  re-integration  routine. 
The  baseline  N/D  has  been  re-integrated  into  a new  wing  geometry.  The  steps 
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Baseline  Nacelle/Diverter  Reintegrated  to  a New  Wing  Geometry 


Diverter  (black)  is  regenerated  to  fit  new  wing 


Improved  Geometry  Capability: 
Automatic  W/B/N/D  Surface  Grid  Generation 
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Improved  Geometry  Capability: 
Automatic  W/B/N/D  Surface  Grid  Gene 
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Automatic  W/B/N/D  Surface  Grid  Generation 
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Improved  Geometry  Capability: 
W/B/N/D  Volume  Grid  Perturbation  (FlexMesh 
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Improved  Geometry  Capability: 
W/B/N/D  Volume  Grid  Perturbation  (FlexMesh 
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Improved  Geometry  Capability: 
W/B/N/D  Volume  Grid  Perturbation  (FlexMesh) 
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Multiblock  Grid  Perturbation  Demonstration 
Integrity  of  Grid  in  the  Presence  of  Multiple  Body  Surfaces 
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Multiblock  Grid  Perturbation  Demonstration 

Integrity  of  Grid  in  the  Presence  of  Multiple  Body  Surfaces 
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Automated  Surface-Grid  Regeneration 
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Improved  Geometry  Capability: 
Automated  Surface-Grid  Regeneration 
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□ ■ "non-slave"  grid  points 
• ■ baseline  grid 
O □ perturbed  grid 


Automated  Surface  Grid  Projection 
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Improved  Geometry  Capability: 
Efficiency  of  FlexMesh 


q 5 
< 


o> 

- S fe 


o $0 


do 


co  _ 
o £ 
io 

5 °- 

ft.  CM 
<X)  O 

CM  1_ 


O CD 


(0  c o 

o o 
0 0 
(0  CO 

co  o> 
^ uo 


CO  0 

E 

Q O 
Q 0 

< CD 


w o o T p 

■o  *=  g 5 to 

O ZJ  .2  Z5 

Z W "S  & ^ 

O O 0 c < 

. C C "O 

O O 0 O '5= 

Z Z CD  O CD 


232 


(CFL3D/FlexMesh)  (40/21)  (96/55)  (105/68) 

timings  are  for  Cray  C-90  computer,  single  processor 
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Alternate  Gradient  Capability 
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Alternate  Gradient  Capability: 
ADIFOR  for  CFL3D 
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Alternate  Gradient  Capability: 
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This  work  was  performed  by  Prof.  A.  Taylor  and  Dr.  A.  Oloso  through  Eagle 
Aeronautics,  Inc.,  Newport  News,  VA. 
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Alternate  Gradient  Capability: 
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CFL3D.ADII  - Accuracy 

For  the  same  benchmark  test  discussed  in  the  previous  chart,  this  table  gives  the  ratios 
of  the  sensitivities  as  computed  with  automatic  differentiation  (labeled  ADII)  and  with 
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Alternate  Gradient  Capability: 
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Alternate  Gradient  Capability: 
CFL3D.ADII  - Design  Demonstration 
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Alternate  Gradient  Capability: 
CFL3D.ADII  - Desian  Demonstrat 
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Alternate  Gradient  Capability: 
CFL3D.ADII  - Design  Demonstration 
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Alternate  Gradient  Capability: 
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- # of  processors  = # of  design  variables 

- about  12  hours  CPU  time/design  cycle 


Alternate  Gradient  Calculation: 
MD03D/CFL3D.ADII  Coupling 
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Alternate  Gradient  Calculation: 
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Design  Output  File 
Design  Restart  File 
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TCA  Wing  Twist  Sensitivity 


Finite-difference:  dCL/dv  = 0.930E-3,  dCD/dv  = 0.100E-3 
CFL3D.ADII:  dCL/dv  = 0.936E-3,  dCydv  = 0.101  E-3 
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Alternate  Gradient  Capability: 
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R (U,X) : Euler  equations;  F (U,X) : Cost  function; 
X (U,X) : Adjoint  variables;  u (X) : Flow  variables; 
x (b) : Grid;  b : Design  variables 
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Summary 


The  implementation  of  the  geometric  modeling  capability  for  use  with  an 
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Future  Work 
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TCA  Nacelle  Installation  Assessment  and 

Design  Studies 
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- Baseline  nacelle/diverter  alignment  for  improved  inlet  flow  quality 

- Boundary  layer  height  assessment 

- Nacelle/diverter  camber  modifications 


Drag  Polars  for  the  TCA  W/B  Configuration 
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Drag  Coefficient,  C 


Wing/Body/Nacelle/Diverter  Navier-Stokes  Grid  Topology  for 
the  TCA  Baseline  Configuration  (Wing/Body  Blocks) 
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Wing/Body/Nacelle/Diverter  Navier-Stokes  Topology 
for  the  Baseline  TCA  Configuration 
(19  Blocks,  6.01  Million  Grid  Points,  Wing/Body  Blocks) 
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Wing/Body/Nacelle/Diverter  Navier-Stokes  Topology  for  the 
Baseline  TCA  Configuration  (Nacelle_Back  Blocks) 
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Wing/Body/Nacelle/Diverter  Navier-Stokes  Topology 
for  the  Baseline  TCA  Configuration 
(19  Blocks,  6.01  Million  Grid  Points,  Nacelle_Back  Blocks) 
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Drag  Polars  for  the  TCA  Baseline  W/B/N/D  Configuration 
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0.0100  0.0150 
Drag  Coefficient,  Cn 


Drag  Polars  for  the  Baseline  TCA  W/B  and  W/B/N/D 

Configurations 
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Pressure  Distributions  for  the  TCA  Baseline  W/B/N/D 

Configuration 
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Pressure  Distributions  for  the  Baseline  TCA 
W/B/N/D  Configuration 

CFL3D,  Moo=2.4,  0=3.518°,  Re  =212x106 
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Lower  Surface  Pressure  Contours  for  the  Baseline  W/B/N/D 

TCA  Configuration 
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Euler  N-S  (Baldwin-Lomax) 


Mach  Contours  for  the  TCA  Isolated  Inboard  Nacelle 

CFL3D  N-S  Solution 
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Mach  Contours  for  the  TCA  Isolated  Inboard  Nacelle, 
Side  and  Top  Views  with  Inlet  Aligned  with  Freestream, 
CFL3D  N-S  Solution,  Moo=  2.4,  Rec=  21 2 xl  0$ 

(3  Block  H-0  Grid,  Baldwin-Lomax) 
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Pressure  Drag  Breakdown  for  the  TCA  Configuration 
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Skin-Friction  Drag  Breakdown  for  the  TCA  Configuration 

M = 2.4,  C,  = 0.1 
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Skin-Friction  Drag  Breakdown  for  the  TCA  Configuration 
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Flat-Plate  N-S  Flat  Plate  N-S  Flat  Plate  N-S 


Nacelle  Installation  Drag  for  the  TCA  Configuration 
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Euler  N-S  Euler  N-S 


Nacelle  Interference  Drag  for  the  TCA  Configuration 
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Re  = 212  x 10B  Re  = 6.36  x 10 


Comparisons  of  Nacelle  Inlet  Mach  Number  Distributions 

TCA  W/B/N/D  Configuration 
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Comparisons  of  Nacelle  Inlet  Mach  Number  Distributions 

TCA  W/B/N/D  Configuration 
(CFL3D  Analyses,  Moo=2.4,  a=3.5°) 
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Inlet  Flow  Angle  Variations  (w.r.t.  to  the  angle  at  the  inlet  center) 

TCA  W/B/N/D  Configuration 
(CFL3D  Analyses,  Moo=2.4,  a=3.5°) 
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Nacelles  and  Diverters  Alignment  Assessment  for  the  Baseline 
Wing/Body/Nacelle/Diverter  TCA  Configuration 
CFL3D  Euler  M = 2.4,  a = 3.5° 
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Alignment  of  Nacelles/Diverters 
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Lower  Surface  Pressure  Contours  for  the  TCA  W/B/N/D 
Configuration  Before  and  After  Nacelle/Diverter  Alignment 


CO 


c o 
0 

co  0 

C~  f—  . . 


c 

o 


CD 

o 

E 

CO 

o 
o 
sz 

CO 

CD 
O) 

“O 

CD 

D> 

C 

TD 

03 

J0 

0 

0 

> 

TD 

“O 

0 

o 

£ 

=3 

O 

CD 

x: 

<u 

;g 

"co 

CO  cii 
3 (0 

o cd 
> o 

CD  -Q 

Q-c 

cd  a) 


c c _ 
cd  £ cd 
aPr 

S 8 “ 
Po 
o o ’l 

1?  s 

>,  X 

Cu  J=  _ 
■n  ra  C 
C *o  Cd 

o c - 

o O CO 

a)  o-o 

“ S 9 
® « * 


<D 

t 

CD 

> 

"a 

"O 

k. 

03 

O 

_Q 

=3 

O 

CD 


C 

CD 

E 

c 

O) 

IB 

a> 


0-2 
7%  O 


CD 


CD 


C — 

as 

■s-s 

cd  o 
O _c 
JD  in 

.£  CD 
CD  £ 


_ a> 

.E  co  c5  p 

•a  cd 
cd 
c 

.2  o 

c 

CD 

E 


w O)  2 

-C  c -Q 

^ £ 

£ So 


CO 

< 


CD 

E 

E 

>* 

CO 


c *a 

CD  0 

CO 

_Q 

O 


0 

sz 

o 

c 

c 2* 

o o 

CO  O 
C JC 
0 CO 
Q_  ^ 
X 0 
0 O) 

£§ 

=3 

00  W 

E « 

CD  JC 

V—  4—* 

e-'i 

03  CO 
O T3 
C 
CO  (D 

E T3 

&§ 

-p-o 

§ % 

S ° 

o »- 

Cd 

:ii 

CO 

cd  cd 
£ o 
o cd 

E Q. 


291 


Lower  Surface  Pressure  Contours  for  the  Wing/Body/Nacelle/Diverter 
TCA  Configuration  Before  and  After  Nacelle/Diverter  Alignment 

CFL3D  Euler  , = 2.4,  a = 3.5° 
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Baseline  Aligned 
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Comparison  of  Boundary  Layer  Profiles  from  W/B  and  W/B/N/D 
CFL3D  N-S  solutions,  Baldwin-Lomax  Turbulence  Model 
ot=  3.518° , Moo=  2.4,  Rec  = 212  x 10<s 
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Estimated  Boundary  Layer  Heights 
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Estimated  Boundary  Layer  Heights  at  Diverter  Leading-Edge 
Locations  on  the  TCA  Baseline  Configuration 
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Nacelle/Diverter  Modification  Study 
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Diverter  Cambering  Procedure  for  the  TCA 

Divl  Casel 
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Procedure  for  the  TCA 
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Inboard  Nacelle  Cambering  Results  for  the  Aligned  TCA 

W/B/N/D  Configuration 
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Inboard  Nacelle  Cambering  Results  for  the  Aligned  TCA 

W/B/N/D  Configuration 
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Conclusion  and  Work  in  Progress 
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Recent  Wind  Tunnel  Tests 
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Recent  Wind-Tunnel  Tests  to  Evaluate  Nozzle  Drag 
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Several  (14+)  2-D  nozzles  representative  of  transonic  flight 
for  HSCT  tested  in  April  1 996 


Comparison  of  Predicted  and  Measured  Drag  Polars 
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Convergence  Histories  for  the  Ref.  H Installed  Axisymmetric 

Supersonic  Nozzle  Configuration 
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Convergence  Histories  for  the  Ref.  H Installed  Supersonic  Nozzle  Configuration 

CFL3D  N-S,  Mj=0.90,  a=4.0°,  Baldwin-Barth  Turbulence  Model,  Re  = 40  x 106 
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Iterations 


Convergence  Histories  for  the  Ref.  H Installed  Axisymmetric 

Transonic  Nozzle  Configuration 
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Iterations 


Upper  Surface  Pressure  Contours  for  the  Installed 

Axisymmetric  Nozzles 
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Upper  Surface  Pressure  Distributions 
Reference  H Installed  Axisymmetric  Nozzle  Configurations 

a=  4 0°  Re  = 40  x 106  CFL3D  N-S,  Baldwin-Barth  Turbulence  Model 
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Transonic  Nozzle  Configuration  Supersonic  Nozzle  Configuration 


Lower  Surface  Pressure  Contours  for  the  Installed 

Axisymmetric  Nozzles 
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Reference  H Lower  Surface  Pressure  Contours  for  the  Installed 
Axisymmetric  Nozzle  Configuration  at  Both  Supersonic  and  Transonic  Settings 

(Moo=  0.9,  a=  4.0°,  Rec=  40  x 1 06,  CFL3D  N-S,  Baldwin-Barth  Turbulence  Model) 
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Transonic  Nozzle  Configuration  Supersonic  Nozzle  Configuration 


Mach  Contours  (Horizontal  Cuts)  with  surface  Cp  Distribution 
for  the  Reference  H with  both  Transonic  and  Supersonic  Nozzle 

Configurations 
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Mach  Contours  (Horizontal  Cuts)  with  Surface  Cp  Distribution  for  the  Reference  H 
with  both  Transonic  and  Supersonic  Nozzle  Configurations 

(Moo=  0.9,  a=  4.0°,  Rec=  40  x 106,  CFL3D  N-S,  Baldwin-Barth  Turbulence  Model) 


330 


Transonic  Nozzle  Configuration  Supersonic  Nozzle  Configuration 


C Distributions  for  Reference  H Installed  Axisymmetric  Nozzle 
P Configurations 
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Cp  Distributions  for  Reference  H Installed  Axisymmetric  Nozzle  Configurations 

Moo=  0.9,  a=  4.0°,  Rec=  40  x 106,  CFL3D  N-S,  Baldwin-Barth  Turbulence  Model 
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Transonic  Nozzle  Configuration  Supersonic  Nozzle  Configuration 


Isolated  2-D/Transonic  Nozzle  Boattail 
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Isolated  Ref.  H 2D/Transonic  Nozzle  Boattail  Configuration 

Computed  Surface  Oil-Flow  Pattern 

(CFL3D  Baldwin-Barth,  M^O.9,  Rec=40x106;  12-zones,  1.5  million  grid  points 
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Isolated  Ref.  H 2D/Transonic  Nozzle  Boattail  Configuration 

Computed  Cp  Distributions 

(CFL3D  Baldwin-Barth,  Rec=40x106;  12-zones,  1.5  million  grid  points) 
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Isolated  2-D/Supersonic  Nozzle  Boattail 
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Isolated  Ref.  H 2D/Supersonic  Nozzle  Boattail  Configuration 

Computed  Surface  Oil-Flow  Pattern 

(CFL3D  Baldwin-Barth,  Rec=40x106;  12-zones,  1.5  million  grid  points) 
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Grid  Structure  for  the  Reference  H Installed  2-D  Nozzle 
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Grid  Structure  for  the  Reference  H Installed  2-D  Nozzle 
Configuration  at  the  Transonic  Setting 

( 37  blocks,  8.33  million  gridpoints,  back-view) 
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Cp  Distribution  for  the  Installed  Reference  H 2-D  Transonic 

Nozzle  Configuration 
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Cp  Distribution  for  the  Installed  Reference  H 2D  Transonic  Nozzle  Configuration 

Moo=  0.9,  a=  4.0°,  Rec=  40  x 1 06,  CFL3D  N-S,  Baldwin-Barth  Turbulence  Model 
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Zoom  View 


Surface  Cp  Distribution 
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Surface  Cp  Distribution 

Installed  Reference  H 2-D  Transonic  Nozzle  Configuration 

a=  4.0°,  Rec=  40  x 1 06,  CFL3D  N-S,  Baldwin-Barth  Turbulence  Model 
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Nozzle  Zoom  View 


Axi  and  2-D  Isolated  Drag  Comparisons  at  Mach  0.9 
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Axi  and  2-D  Isolated  Drag  at  Mach  1.10 
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Axi  and  2D  Nozzle  Drag  Comparisons 

Boattail  and  Sidewall  Pressure  Drag 
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Axi  Isolated  Drag  at  Mach  0.90  and  1.10  - Supersonic 

Geometry 
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CFD  and  IMS  Nozzle  Drag  Comparisons 

Boattail  and  Sidewall  Pressure  Drag 
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Comparison  of  CFD  Codes  for  the  Isolated  2-D  Nozzle  Drag 
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Nozzle  Drag  Comparisons  between  the  CFD  Codes 
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Nacelle  Installation  Increments,  Mach  0.90 
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Axi  and  2-D  Nacelle  Interference  Drag  Comparisons 
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Nacelle  Interference  Drag  Comparisons 


C 

To  ■£ 

■■4=  3 

w 0 © 

— rn  W 


"D 

C 

CO 

© 

-Q 

CO 

o 

> 

CO 


<0 

m 10 

O CD 

0 <D 

^ 05  a> 


0 

co 

£ 


o 

CO 


0 

3 g 
£-§ 
© 3 

0 i2 

F 13 
0 
CD 


CO 

(0 

0 

N 

N 

O 

C 


o 

CD 


0 

0 

0 


> 

0 

i.5 

N (0  ^ 

o -Q  O) 

c ° "D 
•c  O 0 

CO  _ 

X3  h— 
0 0 0 
-C  c ^ 
T c © 
o « £ 

**-  Q-  r- 

S5  w 

2 Z S 

y - o 
a>  <d  i- 

© N ^ 
§ N 


o p 

co  ^ 

2 CD 
^ 05 
CO  CD 


O -© 
C © 
^ _ CO 


vu 

® 9 ^ 

c w 


0)  10 
0 r W 


Q w " 


CVJ 

■D 

C= 

0 

X 

0 

0 


O 

X) 


O) 

0 

"D 

0 

O 

C 

0 

i_ 

0 

T= 

0 


0 

-C 

H 


c 

0 C 
0 

w 0 
CDi= 

2 0 
T3  0 
0 "O 

a 0 
c £ 

£ § 

CD  CO 

*t  a> 

© Q. 
.£  © 
© 3 

-C  to 
~ CD 
C »- 

23 

■*— * -»— * 

0 C 

? * 

.-t;  .c 

c *■' 
co  © 
co  © 

!i 

a)  co 

0 « 

0 0 

o>£ 


z o y 

0 N 


r -c= 


2 c 


0 

-C  0 
_C 

0 ■+-* 

co  c 

<D  0) 

CO  ^ 
0 a) 

<D  -D 
£ © 
5 § 

o 2 

-Q  0 

co  ^ 

CD  TJ 
■°  0 

©I 

E <D 

© 00 

© £ 

O T5 

CO  C 

c •— 

O CO 
■3  0 
0 

I? 


<D  "O 
C © 

C co 

5 8 

■D  3 

.£  i= 
^ 0 
0 £ 
£ "O 


©Cm 


© +-  ->5 
0 £ 0 
.©  jg  o 

"O  +=  'F 
.£0O 
- 7S  © 
.©0  0 
N t Q, 

N $ 3 

0^0 

r-  0 

o ;e 

■m  o ^ 

ii« 
M 1 

x © 2 
0 © «- 

-n  CD  © 
0 © £ 

0 O 
0 00  ° 

.£  3 c 
© o ° 

.c  -©  ~ 

0 *~ 

o © o 
a .Q  0 


>,§ 

0 4= 

’*  7:3 

0 c 


o « 

sj 

S .2 

o>i 

0 £ 

•D  © 

0 0 
o — * 
c 0 

0 -C 

2 i 

.£  -o 

•♦r  © 

0 *5< 

c o 

13  O 


o je  3 E 


E o 0 
- - 0 


S 

|2 

5 2 Eg 

Z;  h—  rn  X 


— 0 

2 - 

il 


o 2 

2 3 
0 0 
U L. 

5 0 
0 rE 
c 

0 
0 0 


m ° © 

©0-0 


© 


2 8 2 
I 53 

C .0T3-a 

— «»  o c c 


— 0 

■55  © 
© 0 

£ © 

0 O 
"SJ  1 
© (M 

-E*  0 

3 -C 

_ +—t  4-* 

C 

© o 


^ © 


p ? © 

r—  x: 
© o **- 

©CO 

K .©  © 

o — S? 
g © 0 

P O 0 

5 F o 


359 


Configuration  Aerodynamics  Workshop,  25  February  1 997 


Axi  and  2D  Nacelle  Interference 
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CFD  and  IMS  Jet  Effects  Drag  Comparisons, 
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CFD  and  IMS  Jet  Effects  Drag  Comparisons 
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2-D  Fixed  Chute  Nozzl 
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Comparison  of  Current  FCN  and  ATC  Nozzle  Drag 

Nozzle  Pressure  Drag  Predicted  by  IMS 


Configuration  Aerodynamics  Workshop,  25  February  1 997 


1^ 

o> 

o> 


(0 


■Q 

0)|£ 

<0 

^ LO 
D CM 

® O 

N 

N 2 
O 

Z 0) 

o c 

c .2 

o 


<0 


<0 


C c 
CC  0) 
£ CO 

T3  O 
0)  ^ 

m C 

5 0) 

(/)  t 

C 3 

o 

T3 

C 

(0 

(0 

3 

s 

55 


CD 

■D 

CO 

E 

CD 

XI 


CO 

C 

o 

03 

£ 

<D 

CO 

JD 

O 

C 

CO 

■c 

o 

Q. 

E 

CD 

0 

CD 

k_ 

0 


3 

-Q 

<D 

■4— » 
0 

Q_ 

E 

o 

o 

0 

>s 

o 

c 

CO 

>» 

“O 

3 

■4— ' 
0 

0 


374 


A npnrlimnmiAA  lA/Aifl/rhAn  OK  Pohri  ior\/  1007 


Installed  Transonic  Nozzle  Drag 
Status  and  Current  Observations  as  of  25  February  1997 
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Uncertainties  in  HSCT  Cruise  Drag  Prediction 

Shreekant  Agrawal,  Michael  G.B.  Novean,  Geojoe  Kuruvila,  and  Robert  P.  Narducci 

McDonnell  Douglas  Corporation 
Long  Beach,  California  90807-5309 


Within  the  High  Speed  Research  (HSR)  program,  NASA  and  Industry  are 
jointly  developing  various  technologies  so  that  the  U.S.  Industry  has  the 
capability  to  launch  the  High  Speed  Civil  Transport  (HSCT)  aircraft 
development  in  early  2000.  One  of  the  primary  objectives  of  the  HSR 
program  is  to  be  able  to  predict  the  cruise  aerodynamic  performance  of  the 
HSCT  configurations  with  a sufficiently  high  confidence  level  that  will  aide 
Industry  in  the  decision  to  proceed  with  the  development  of  the  aircraft  and 
guarantee  its  performance  to  its  airline  customers. 

This  paper  addresses  the  current  status  in  the  prediction  of  drag  at 
primarily  the  supersonic  cruise  Mach  number  (M  = 2.4),  however,  drag  at 
the  transonic  cruise  Mach  number  (M  = 0.9)  is  also  presented,  wherever 
appropriate.  The  thrust  of  this  paper  is  the  uncertainty  (or  the  confidence 
level)  in  drag  prediction.  Use  is  made  of  the  available  experimental,  linear 
and  nonlinear  computational,  and  empirical  database  to  the  McDonnell 
Douglas  Corporation  (MDC).  In  some  cases,  there  is  sufficient  database;  in 
some  other  cases,  there  is  very  little  database;  and  yet  in  some  other  cases, 
there  is  none  available.  However,  an  attempt  is  made  to  see  where  we  stand 
today  in  the  cruise  drag  prediction,  although  it  is  difficult  to  determine 
uncertainty  levels  in  all  the  elements  contributing  to  drag.  Please  note  that 
the  uncertainty  levels  discussed  here  are  the  views  of  the  researchers  at 
MDC  only,  and  they  may  not  represent  those  at  other  organizations. 
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Contributing  Factors  in  HSCT  Cruise  Performance  Prediction 
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Contributing  Factors  in  HSCT 
Cruise  Performance  Prediction 
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• Uncertainty  in  prediction  (primarily  at  M = 2.4)  of 
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Comparison  of  Predicted  and  Experimental  Drag  Polars 
TCA  Wing/Body  and  Wing/Body/Nacelle/Diverter  Configurations 
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Comparison  of  Predicted  and  Experimental  Drag  Polars 
TCA  Wing/Body  and  Wing/Body/Nacelle/Diverter  Configurations 

Langley  Test  1671,  UPWT-2,  Moo=2.4,  Rec=6.36x106,  Trip  Drag  Removed 
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Trip  Drag  Assessment 
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Trip  Drag  Assessment  Procedure 
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Uncertainty  in  Trip  Drag  Measurement 
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Trip  Drag  Study,  Extrapolation  to  k=0" 
1.675%  TCA  Baseline  Model,  W/B  Configuration 
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Trip  Drag  Study,  Drag  Polars 

1 .675%  TCA  Baseline  Model,  W/B  Configuration 

Langley  UPWTTest  1671,  M_=2.4,  Re=4x106/ft 
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Trip  Drag  Study,  Drag  Polars 
1.675%  TCA  Baseline  Model,  W/B  Configuration 
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Uncertainty  in  Trip  Drag  Measurement 
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Trip  Drag  Study,  Extrapolation  to  k=0" 

1 .675%  TCA  Baseline  Model,  W/B  Configuration 

Langley  UPWTTest  1671,  Mj=2A,  Re=4x106/ft 
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Trip  Drag  for  k=0.012" 
1.675%  TCA  Baseline  Model 
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Trip  Drag  for  k=0.012" 

1 .675%  TCA  Baseline  Model 

Langley  UPWTTest  1671,  M =2.4,  Re=4x106/ft 
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Lift  Coefficient,  C 


Uncertainty  in  Trip  Drag  Estimation 
1.675%  TCA  Baseline  Model,  W/B  Configuration 
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Uncertainty  in  Trip  Drag  Estimation 
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- Increasing  the  number  of  trip  heights  tested 

- Increasing  the  number  of  repeat  runs  at  a specific  k (i.e.,  reduce  ACD) 


Pressure  Distributions  for  the  TCA  Wing/Body  Configuration 
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Pressure  Distributions  for  the  TCA  Wing/Body  Configuration 
UPWT  1671,  CFL3D  Euler  and  Navier-StokesCp 
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Pressure  Distributions  for  the  TCA  Wing/Body/Nacelle/Diverter 
Configuration:  Experimental  Data,  Euler,  and  Navier-Stokes 
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and  reflected  shocks  are  nicely  predicted  as  well. 
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Pressure  Distributions  for  the  Ref.  H Wing/Body/Nacelle/Diverter 
Configuration:  Experimental  Data  and  Navier-Stokes 


<D  <D 
co  o 

22  C 
X <D 

H 43 


3 Oh 

Uh 

g.s 

^ & 

o a 

§ 3 

•e  s 

•-  Uh 

£ U 

c ^ 

r a <D 

U 2 

H03 

X) 

CD  oo 

. o 


03 

33  M 
C 

o ‘X 
X <L> 

0)  CD 
L~r 
C/5 

- g 

C3  v-i 

o 5 

00  ^ 

*c  ^ 

03  >* — i 

Cl  *“■  "Tl 
c S? 

g > £ 

o _ 

O C3  OJ 

4U  §. 

w bD 

, 30) 

C J2  C3 

•-J  C u, 

03  Q -o 
M oo  ^ 

<3  ^ 

C3  3 

CL  £ 

S 00 
<D 

O g 

O CL 


406 


0.05i 


Pressure  Drag  at  Supersonic  Condition 
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Pressure  Drag  at  Supersonic  Condition 
Navier-Stokes  (Baldwin-Lomax),  M = 2.4,  CL  =0.1 
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Skin-friction  Drag  for  TCA  W/B,  Isolated  Nacelle  and  W/B/N/D 
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Skin-friction  Drag  for  TCA  W/B,  Isolated  Nacelle  and  W/B/N/D 

M = 2.4,  C.  =0.1 

OO  7 L 
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W/B  4 Nacelles  W/B/N/D  W/B  4 Nacelles  W/B/N/D 

Rec  = 6.36x10e,  Truncated  Fuselage  Rec  = 212x106,  Full  Fuselage 


Skin-friction  Drag  for  2.2%  Ref.  H Wing/Body  Configuration 
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Prediction  of  Nacelle  Installation  Drag 
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Prediction  of  Nacelle  Installation  Drag 

CFL3D  & Experiment,  M =2.4 
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Ref.  H TCA  Ref.  H,  C =0.12  TCA,  C,  =0.10 


Prediction  of  Isolated  Nozzle  Boattail  Drag 
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Uncertainty  in  Predicting  Nozzle  Boattail  Drag 
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Axisymmetric  2-D 


Geometries  of  the  Rigid  and  Flexible  Opt  5 Model 
Corresponding  to  the  Measured  Structural  Deflections  at  oc= 


o 


Cd 


cd 


<L> 


J .5 


ao 


C/5 


•g  s 


c/5  cd 
0) 

« » 

ao  c 
c <u 
•C 

3 T3 
Q § 


r^-H 

X3 
• ^ 

.2? 

C 

<o 

dd 

+-» 

o . 

£ GhT3 

~ "O  '5b 

*c 

aj 

u- 

V 

D ’C 
£ ao 

Ui 

G 

<d 

r* 

£ 

ao  53 

• f-H 

cd 

a 

C/5 

G sz 

G 

G 

C/5 

C/5 

T3 

'C 

o ° 

/— N * 'l 

G 

CD 

Lh 

03 

X!  Q 
OOpLi 

-C  4) 
*5  JG 
cd 

<U  c/5 

T3  TD 

O cd 

s 5 


‘-c  2 


. jS 
>>  00 

X>  p 

H 

£ 

Oh 

D 

<u 

dd 


o 

Ui 

dd 

w 

C/5 

a> 

jd 

Id 

> 


CD 

T3 

O 

s 

T3 
• ^ 

&Q 


T3 

s g 

G X3 

G 

cd 


c/5 

cd 


C/5 

cd 

£ 


Uh 

&0 


418 


o 


O 

"U 

o 


II 

a 

73 

</> 


LO  O 

Q.® 

o *5) 

0)  Q 

5 75 

X T 

a)  2 

si 

■d'o 

■-  <D 

.5>5 

cc  c/> 

© 2 


o © 
</)  £ 
.2  o 

O O) 

if 

0)  c 
0O 

(/) 

0 


o 

o 


419 


wing  twist  and  bending  data  available 
at  r)=84.5  and  T)=96. 1 percent  stations 


Aeroelastic  Study  Drag  Polars 
M2.4-7A  Opt  5 1.675%  Model,  Wing/Body  Configuration 
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Aeroelastic  Study,  Drag  Polars 

M2.4-7A  Opt5  1 .675%  Model,  W/B  Configuration 

CFL3D  N-S,  Baldwin-Lomax,  M =2.4,  Re=4x106/ft 
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Drag  Coefficient,  C 


Effect  of  Reynolds  Number  on  Pressure  Drag 
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Effect  of  Reynolds  Number  on  Pressure  Drag 
CFL3D-NS  (Baldwin-Lomax),  = 2.4,  CL  =0.1 
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W/B  W/B/N/D 


Effect  of  Reynolds  Number  on  Pressure  Drag  at  Cp 


o3  0)  s-h  (D 

a x n 

m c3 

O — <U  . 

/j  ^ «i  >% 

C 3 <U 

s * £ -5 

*g  S e 

T3  2 Oh  O 
"3  X « 3 

03  v©  -O  i_i 
w cn  w 3 


J3  2 

O w 3 
c«  ^ 00 

<U  !+h 

<U  X5  c 

^ g o 
o § Cj 

75  C >> 


£ § l 
z sr| 

Q OfS  * 

CO  <D 

J J3  cS 

Cl  c ^ 

n c ^ 
U a o 


^ > <D 

_H  > +-1 

CJ 

CU  QJ  0 

C X 00 
X +->  <3-> 

03  . t-4 

X)  (o  ^ 
O t4-1  <u 

g ^ «3 

g <N  ^ 

J3  ii  *> 

3 8 O 

> *5  — < 

®P  h-j  04 

g 3 - 


u,  _Q 

6*2  6 

3 , — . 3 

£ 13  c 

;s-§  ^ 

£ £ § 
<u  <u  >L 
.c  g « 
H § d! 


o . 

£•£ 

^ a 

<D  \g 
^ C 
^ <D 


S ^ 

— c3 
a <u 
> c 


424 


Effect  of  Reynolds  Number  on  Pressure  Drag  at  CD 
CFL3D-NS  (Baldwin-Lomax),  M„  = 2.4,  C,  » 0.01 
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W/B  W/B/N/D 


Effect  of  Reynolds  Number  on  Pressure  Drag 
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Effect  of  Reynolds  Number  on  Pressure  Drag 
TLNS3D  (Baldwin-Lomax),  M = 0.9,  CL  =0.1 8 
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Reynolds  Number  Effects  on  Pressure  Distributions 
TCA  Wing/Body/Nacelle/Diverter  Configuration 
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Reynolds  Number  Effects  on  Pressure  Distributions 
TCA  Wing/Body/Nacelle/Diverter  Configuration  Cp 
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Drag  Prediction 
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Overall  Uncertainties  in  =0.90  Drag  Prediction 
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Overall  Uncertainties  in  M^O.9 
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Drag  Build-up  from  CFD  to  Flight  Condition 
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Drag  Build-up  from  Wind-Tunnel  to  Flight  Condition 
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Drag  Build-up  from  Wind-Tunnel  to  Flight  Condition 

TCA  Drag  Build  up  Components  (M=2.4,  CL=0.1) 
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CFD  and  wind-tunnel  testing  help  to  check  the  TI  data  base  and  provide  data  for 
refinement  of  assumed  drag  components  and  the  technology  projection 


Conclusions  are  summarized  in  the  next  two  slides. 
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- Aeroelastic  effects  at  flight  Re 
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Goal:  Update  S&C  assessment  with  Cycle  3 evaluation 
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Configuration  Aerodynamics  (Task  32) 
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Reference  H Cvcle  3 analytic  model: 
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Flight  Path  Stability 


Donnfa  Henderson 


Bare  Airframe  Flight  Path  Stability 

The  FCS  Requirements  document  states  that  "The  curve  of  flight  path  angle  versus  true  airspeed  for  the  unaugmented 
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atisfied  for  end-of-cruise  weight  at  all  altitudes. 
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Descent  Capability 
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Emergency  Descent  61 4K  Emergency  Descent  501 K 


Go  Around 
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Requirement  Satisfied 
+ 6 deg/sec2  pitch  acceleration  satisfied. 

+ 3 degree  pitch  change  in  1 sec.  not  satisfied 
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High  a Recovery 
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Requirement  of  -4.0  deg/sec2  NOT  satisfied 
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Diving  Pullout 


462 


Requirement  written  for  partial  failure,  but  model  does  not 
allow.  Satisfied  for  healthy  A/C. 


The  diving  pullout  requirement  states  “It  shall  be  possible  to  effect  a 1.5g  pull-out  at  the  V,-/MD  boundary  with  one 
actuator  power  supply  channel  inoperative”.  Currently,  the  Simulink  model  does  not  have  the  capability  of  partial 
hydraulic  system  failures;  therefore,  this  assessment  was  conducted  with  full  hydraulic  power  available.  Diving  pullouts 
were  simulated  for  a heavy  cruise  weight  of  614.9  klbs.  and  a 1 .5g  pull-up  was  attainable  throughout  the  envelope.  The 
figure  illustrates  time  histories  for  a diving  pullout  at  50,000  ft. 
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Roll  Performance 
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Configuration  1 : Baseline,  TE  1/8, 3/6,  and  SSDs 
- Level  1 FQ  NOT  satisfied 
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• Can  trim  at  wings  level,  but  not  y = 
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Emergency  Descent  rate  requirement  for  all  weight  cases  in  the  transonic 
and  supersonic  regimes.  In  the  subsonic  regime,  the  requirement  is  met  at 
end  of  cruise  weight,  but  is  NOT  met  at  mid-cruise  weight  and  heavier. 


An  assessment  of  the  Reference  H configuration  using  the  Cycle  3 simulation  model  has  been  conducted  and  the  results 
have  been  compared  to  the  FCS  requirements.  Primary  consideration  was  given  to  evaluating  those  requirements  that 
pertain  to  the  HSR  Configuration  Aerodynamics  task.  Assessments  currently  indicate  that  the  Reference  H aircraft: 
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Oust  Recovery 
Maneuvering  Performance 
Crosswind  Landing 
Tameness 


Summary  - cont. 
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the  "initial  guess".  Alternative  schemes  should  be  investigated. 


Drag  devices  on  the  HSCT  must  be  carefully  sized  and  configured,  as  it  was  shown  that  200  incremental  drag  counts 
the  baseline  configuration  were  required  for  emergency  descent  capabilities. 
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Outline 

• Objectives 

• Background 

- description  of  transition,  location  of  trips,  and  trip  drag 

- supersonic  vs  subsonic  issues 

• Approach 

- multiple  tests,  multiple  trip  types,  patterns,  heights 

- low  and  high  Rn,  transition  detection 

• Low  Rn  Supersonic  Testing  (UPWT) 

• High  Rn  Supersonic  Testing  (PSWT) 

• High  Rn  Subsonic  & Transonic  Testing  (NTF) 

• Concluding  Remarks 


As  outlined  above,  this  presentation  describes  the  general  objectives  of  the 
project,  followed  by  background  information  which  led  to  the  initiation  of  the 
study,  and  the  approach  taken  to  meet  the  objectives.  Next,  experimental 
studies  in  the  LaRC  Unitary  Plan  Wind  Tunnel,  the  MDA  Polysonic  Wind 
Tunnel,  and  the  National  Transonic  Facility  will  be  discussed.  Concluding 
remarks  will  close  the  presentation. 
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Objectives 

• Determine  the  best  method  (grit  vs  dot)  of  tripping  the  boundary 
layer 

• Determine  the  best  method  for  assessing  trip  effectiveness 

• Determine  the  best  method  for  quantifying  trip  drag 

• Determine  any  advantage  of  testing  at  slightly  higher  Rn«yft 
available  in  industry  blowdown  facilities 

• Determine  the  Reynolds  number  for  fully  turbulent  flow 


The  objectives  of  this  effort  were  to  determine  (if  possible)  the  best  method:  1) 
for  forcing  the  boundary  layer  to  transition,  2)  for  assessing  trip  effectiveness, 
3)  for  quantifying  trip  drag,  4)  for  testing  at  Reynolds  numbers  per  foot  from  5 
million  to  maximum  available  rather  than  1 to  5 million,  and  5)  for  boundary 
layer  state  determination. 
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Subsonic  vs  Supersonic  Issues 


Subsonic 

- Rnk=600 

- possible  to  size  for  no/minimal  trip  drag  —trip  usually  inside  b.l. 

- works  in  small  alpha  range  — sized  and  located  based  on  a given 
attachment  line  location 

Supersonic 

- Rnk  =1000 

- possible  to  size  for  minimal  trip  drag  —trip  usually  outside  b.l. 

- works  in  small  alpha  range  — sized  and  located  based  on  a given 
attachment  line  location 

Trip  Verification  Methods 

- flow  visualization  — time-consuming 

- Reynolds  number  sweeps  — less  of  impact  on  schedule,  in  general 


It  has  been  shown  (Braslow,  Hicks,  et.  al.)  that  for  subsonic  through  low 
supersonic  conditions,  the  trip  drag  can  be  considered  negligible  and  doubling 
the  effective  trip  height  (as  long  as  the  trip  does  not  stick  into  the  ffeestream) 
will  add  an  indiscernible  amount  of  drag.  However,  at  supersonic  conditions, 
this  is  not  the  case.  In  order  to  transition  the  boundary  layer,  the  trips  must  be 
sized  larger  than  the  boundary  layer  thickness.  Now,  the  portion  of  the  trips 
that  stick  out  into  the  freestream  do  produce  measurable  amounts  of  drag. 

Subsonically,  Braslow  et  al.  have  shown  that  trip  drag  (in  the  plateau  region, 
which  is  the  region  where  the  boundary  layer  has  transitioned  and  the  trips  are 
not  sticking  out  into  the  freestream)  can  be  calculated  by  testing  several  trip 
heights  and  plotting  the  associated  drag  versus  the  trip  height.  A linear  curve 
can  be  drawn  through  the  data  in  the  plateau  region  and  extrapolated  to  k=0 
(no  height).  The  trip  drag  then  is  equal  to  the  slope  of  the  curve  multiplied  by 
the  specific  trip  height  in  question.  Braslow  et  al.  also  have  shown  that 
supersonically,  the  best  method  is  to  use  drag  versus  k2.  A slope  is  then  found 
in  a similar  fashion,  however,  due  to  the  fact  that  drag  on  an  object  in  a 
freestream  is  not  linear  with  height  (it  is  closer  to  quadratic)  a much  lower  drag 
is  estimated. 

Finally,  sublimation  and  other  flow  visualization  techniques  take  a lot  of  time, 
whereas  running  Rn  sweeps  takes  a relatively  short  time.  By  accurately 
reducing  the  data,  a very  good  understanding  of  the  boundary  layer  state  can  be 
made. 
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CD  vs  Rn 


The  ideal  method  of  insuring  that  turbulent  boundary  layer  conditions  are  met  is  shown  above 
(open  circles).  The  technique  would  require  that  the  test  be  run  at  Rn  conditions  that  would 
insure  that  one  would  have  fully  laminar  b.l.  conditions  at  the  low  end  of  the  Rn  capability  to 
fully  turbulent  b.l.  conditions  at  the  high  end.  This  technique  can  be  done  in  some  facilities 
(NTF,  20-Inch  Mach  6,  PSWT,  etc.).  However,  for  most  tunnels,  the  operating  Rn  range  lies 
completely  in  the  transition  region  (shown  above  as  the  location  that  the  data  does  not  lie  on 
the  fully  laminar  nor  the  fully  turbulent  curves).  Therefore,  one  must  trip  the  boundary  layer. 

The  figure  above  also  illustrates  the  effect  of  a transition  trip  on  the  values  of  CD  >miD.  In  the 
same  operating  range,  fully  turbulent  boundary-layer  conditions  can  be  achieved  and  the  value 
of  trip  drag  is  equal  to  the  difference  in  the  forced-transitioned  and  the  free-transitioned,  fully- 
turbulent  b.l.  CD,min.  However,  if  one  can  not  ever  reach  free- transitioned,  fully-turbulent  b.l. 
conditions  in  the  tunnel,  how  can  one  determine  the  trip  drag?  Note,  it  is  always  a good 
practice  to  run  through  a Rn  sweep  in  the  facility  at  Co.min  conditions,  just  to  determine  that  the 
transition  trips  are  working  (and  this  usually  can  be  done  with  minimal  impact  on  the  test 
schedule).  Why  is  this  important  and  what  useful  data  can  you  get?  First,  as  shown  by  the 
solid  circles  above,  CD.min  increases  until  fully  turbulent  conditions  are  reached  as  Rn  is 
increased.  As  Rn  is  increased  further,  the  CD,min  values  will  decrease  (and  the  data  will  fall  on 
a curve  that  is  parallel  with  the  fully  turbulent  b.l.  curve).  Thus,  while  CDmin  is  still  increasing 
or  has  reached  a plateau,  transitioning  b.l.  conditions  are  occurring.  If  the  CD  ^ values 
decrease  at  a fairly  "constant"  rate  with  Rn  (same  local  slope  as  that  for  fully  turbulent  flow 
found  for  free  transition),  then  fully  turbulent  b.l.  conditions  have  been  reached. 
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Determination  of  Trip  Drag 

• Trip  Drag  may  be  determined  by  calculating  the  slope  of  the 
curve  for  drag  (at  a specific  lift  condition)  versus  k or  k2. 


ki  ki  kj  k > 6 kjl  k>J 

K,  trip  tvaight  kJ,  tnp  height  »qu»r«l 


There  are  several  methods  of  determining  trip  drag.  When  a large  Rn 
capability  is  achievable  in  the  facility,  the  technique  described  previously  is 
preferred.  However,  the  reduced  capability  test  facility  is  usually  the  norm.  In 
such  a facility,  one  can  test  the  model  with  different  sized  trips.  Whether  the 
trips  are  smaller  or  larger  than  the  boundary  layer  determine  how  the  trip  drag 
is  calculated. 

In  subsonic  to  low  supersonic  flow  conditions,  the  grit  particles  when  sized 
and  located  properly  on  the  wing  for  the  freestream  conditions,  can  allow  a 
doubling  of  the  size  of  the  trip  to  have  very  little  effect  on  the  drag. 

Note,  once  the  boundary  layer  is  fully  turbulent  and  as  long  as  the  trip  height, 
k,  is  not  greater  than  the  boundary-layer  thickness,  8,  then  the  values  of 
CD,min  are  not  overly  affected  by  the  trip.  Therefore,  a trip  drag  can  be 
calculated  by  extrapolating  a linear  curve  fit  (of  the  CD,min  data  for  the  k<  8 
conditions)  back  to  k=0.  This  will  define  a CD,min  that  should  be  equivalent 
to  fully  turbulent  flow  with  no  trip  drag.  Then  the  trip  drag  for  any  k value  can 
be  calculated  by  subtracting  the  CD,min  value  for  k=0  from  the  value  of 
CD,min  for  a given  k.  Also,  note  that  when  k>  8 , the  CD,min  values  for  those 
k values  increase  rapidly.  Thus,  for  conditions  that  require  the  size  of  the  trip, 
k,  to  be  greater  than  the  boundary  layer,  a similar  method  can  be  employed, 
however,  one  must  use  k2  instead  of  k to  get  an  accurate  extrapolated  curve  to 
k or  k2  = 0 and  thus,  CD,min  conditions  for  fully  turbulent  b.l.  with  no  trip 
drag. 
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The  figure  above  shows  how  the  actual  distance,  xstreamiine,  that  the  flow  covers 
from  the  leading  edge  to  the  trip  dot  (by  following  the  streamline  rather  than 
the  axial  or  "streamwise"  direction  behind  the  leading  edge)  can  change  the 
height,  k,  required  to  transition  the  boundary  layer.  For  highly  swept  wings, 
Xstreamiine  can  be  much  less  than  x and  thus,  the  Rnk  cr  is  never  reached  and  the 
trip  dot  does  not  have  any  effect  on  the  boundary  layer  state. 
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Transition  Detection  Methods 


Trip  Dot 


Trip  Dot 

imbedded  in  ftourine  layer 


Trip  Dot 

Glued  on  top  of  TSP 


The  figure  above  shows  the  effective  height,  k^r,  after  the  sublimating 
chemical  is  applied  to  the  surface.  Since  keff  is  smaller  than  k,  the  value  of 
Rnk,cris  never  reached  and  the  trip  dot  does  not  have  any  (or  very  little)  effect 
on  the  boundary-layer  state.  Thus,  a method  of  applying  sublimating 
chemicals  that  does  not  cover  the  leading  edge  up  to  the  trips  may  be  required. 


Shown  to  the  right  (above)  is  a sketch  of  the  same  trip  dot  glued  onto  a surface 
with  TSP  (Temperature  Sensitive  Paint).  Note,  here  the  effective  height  of  the 
trip  dot  is  the  same  as  that  specified.  In  addition,  the  paint  surface  is  hard,  thin, 
smooth,  and  consistent  throughout  a series  of  runs.  Note  that  paint  does  not 
have  to  be  reapplied  for  each  angle-of-dttack. 
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Approach 

• The  M2.4-7A  Arrow  Wing  was  tested  in  UPWT  and  generated  a 
data  base  at  M~=2.4  and  Rru,  = 1 to  5 million/foot  with  free 
(natural)  and  fixed  (grit  & dot  of  various  size)  transition 

• Test  the  M2.4-7A  Arrow  Wing  in  the  MDA-E  Polysonic  Wind 
Tunnel  (PSWT)  at  Rn~  greater  than  5 million/foot  foot  with  free 
(natural)  and  fixed  (grit  & dot  of  various  size)  transition 

• Obtain  force  and  image  data  for  free  (natural)  and  fixed 
transition 


The  M2.4-7A  Arrow  wing  configuration  was  chosen  because  it  was  a clean 
model,  had  no  “planned”  entries  in  the  near  term,  and  needed  a trip  drag 
analysis  at  the  cruise  condition. 

The  model  is  slightly  larger  than  ideal  for  testing  in  the  PSWT  due  to  its  large 
projected  area  relative  to  start-up  and  ending  conditions.  The  model  is  of 
typical  size  for  supersonic  testing  in  the  HSR  program,  and  results  in  the 
PSWT  were  very  promising. 
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BOUNDARY  LAYER  HEIGHT  AT  TWO  X-LOCATIONS 
LAMINAR  FLAT  PLATE  CALCULATION 


M = 2.40 


5,  inches  o.ow 
(B.L.  Height) 


- **0.42'  ttrMmww*  (0.2*  normal) 
(Outboard  panal:  a»«i  S*) 

> (0.2*  normal) 
: A-71*) 


—0.01  *7*  Tnp  Dot 
—0.0155*  Tnp  Dot 
—0.0121*  Tnp  Dot 

♦ 0.01 00*  Tnp  Dot 
—0.0101*  Trip  Dot 

0.0151*  Sand  Grt 
-Q.01S2*  Sand  Gnt 
-0.012i"  Sand  Qnl 

* 0.O10r  Sand  Grit 


3t 


2.0  3.0 

Re  (x10*/ft) 


Shown  above  is  a plot  of  the  predicted  laminar,  flat-plate,  boundary-layer 
height  at  x = 0.42”  and  0.61”  from  the  leading  edge  in  the  streamwise  direction 
(outboard  panel  trip  location  and  inboard  panel  trip  location,  respectively). 
Since  the  free  stream  Mach  number  is  2.4,  the  required  trip  height  would  have 
to  be  greater  than  the  boundary  layer  thickness  to  be  effective.  Thus,  the  plot 
also  includes  the  sand  grit  sizes  and  trip  dot  heights  that  were  tested  on  the 
M2.4-7A  configuration  in  the  NASA  LaRC  UPWT. 
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TEST  1667  SAND  GRIT  DRAG  DATA 
Baseline  M2.4-7A  Configuration 
M = 2.4,  C ~ 
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2.0  3.0 

Re  (x  10*/ft) 


4.0  5.0 


The  plot  above  shows  the  Co.min  values  plotted  versus  Reynolds  number  for  the 
three  sand  grit  sizes  tested  as  well  as  the  free  transition  values.  Also  plotted 
are  the  flat  plate  fully  turbulent  skin  friction  drag  estimates  corrected  to  the 
Reynolds  number  4 million/foot  condition. 

The  results  show  that  for  Reynolds  numbers  less  than  2 million/foot,  a 
transitional  stage  of  the  boundary  layer  is  occurring. 
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The  plot  above  shows  the  CD,min  values  plotted  versus  Reynolds  number  for  the 
five  trip  dot  heights  tested  as  well  as  the  free  transition  values.  Also  plotted 
are  the  flat  plate  fully  turbulent  skin  friction  drag  estimates  corrected  to  the 
Reynolds  number  4 million/foot  condition. 

The  results  show  that  for  Reynolds  numbers  less  than  2 million/foot,  a 
transitional  stage  of  the  boundary  layer  is  occurring  as  was  seen  with  the  sand 
grit  runs. 
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Baseline  M2.4-7A  Configuration 
(Free  Transition) 

M = 2.40,  Re  = 1 to  5 x10s/ft 


-3.000  -2.000  -1.000  0.000 


The  data  above  shows  the  curve  fit  values  of  the  Axial  Force  Coefficient  (Ca) 
versus  angle  of  attack,  a,  for  the  baseline  M2.4-7A  configuration  at  a Mach 
number  of  2.40  and  Reynolds  numbers  varying  from  1 to  5 million  per  foot  for 
the  free  transition  condition.  The  CA  data  is  very  valuable  in  determining  if  a 
“laminar  bubble”  occurs  at  any  conditions  in  the  polar  as  can  be  clearly  seen 
for  the  lowest  Reynolds  numbers  tested.  As  seen  above  for  all  Reynolds 
number  conditions,  the  curves  are  not  continuous.  By  examining  the  local 
slope  of  the  curves  (dCA/da),  this  effect  can  more  readily  be  observed. 
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The  above  plot  illustrates  how  the  derivative  of  Ca  (dCA/da)  varies  with  angle 
of  attack  for  the  free  transitioning  boundary  layer  over  the  Reynolds  numbers 
tested.  The  data  clearly  shows  changes  in  the  derivative  for  the  lowest 
Reynolds  number  conditions  and  fairly  large  variations  for  the  higher  Reynolds 
numbers  as  well. 
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The  data  above  shows  the  curve  fit  values  of  the  Axial  Force  Coefficient  (Ca) 
versus  angle  of  attack,  a,  for  the  baseline  M2.4-7A  configuration  at  a Mach 
number  of  2.40  and  Reynolds  numbers  varying  from  1 to  5 million  per  foot  for 
forced  b.l.  transition  using  sand  grit  of  nominal  size  0.0181”.  The  Ca  data 
above  show  that  at  a Reynolds  number  of  1 million/foot  the  curve  does  not 
look  smooth  (implying  that  the  boundary  layer  has  not  fully  transitioned). 
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The  above  plot  illustrates  that  the  derivative  of  CA  (dCA/da)  for  the  Reynolds 
number  of  1 million/foot  appears  wavy  implying  transitioning  is  occurring. 
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Baseline  M2.4-7A  Configuration 
(Trip  Dots  = 0.01 47") 

M = 2.40,  Re  = 1 to  5 x10*/ft 
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The  data  above  shows  the  curve  fit  values  of  the  Axial  Force  Coefficient  (CA) 
versus  angle  of  attack,  a,  for  the  baseline  M2.4-7A  configuration  at  a Mach 
number  of  2.40  and  Reynolds  numbers  varying  from  1 to  5 million  per  foot  for 
forced  b.l.  transition  using  trip  dots  of  nominal  height  0.0147”.  The  CA  data 
above  show  that  at  a Reynolds  number  of  1 million/foot  the  curve  does  not 
look  smooth  (implying  that  the  boundary  layer  has  not  fully  transitioned). 
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The  above  plot  illustrates  that  the  derivative  of  CA  (dCA/da)  for  the  Reynolds 
number  of  1 million/foot  appears  wavy  implying  transitioning  is  occurring. 
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Baseline  M2.4-7A  Configuration 
(Sand  Grit  = 0.0107") 

M_  = 2.40,  Re_  = 1 to  5 x10*/ft 


The  data  above  shows  the  curve  fit  values  of  the  Axial  Force  Coefficient  (Ca) 
versus  angle  of  attack,  a,  for  the  baseline  M2.4-7A  configuration  at  a Mach 
number  of  2.40  and  Reynolds  numbers  varying  from  1 to  5 million  per  foot  for 
forced  b.l.  transition  using  sand  grit  of  nominal  size  0.0107”.  The  Ca  data 
above  show  that  at  a Reynolds  number  of  1 million/foot  the  curve  does  not 
look  smooth  (implying  that  the  boundary  layer  has  not  fully  transitioned). 
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The  above  plot  illustrates  that  the  derivative  of  Ca  (dCA/da)  for  the  Reynolds 
numbers  of  1 and  1.5  million/foot  appear  wavy  implying  transitioning  is 
occurring. 
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The  data  above  shows  the  curve  fit  values  of  the  Axial  Force  Coefficient  (Ca) 
versus  angle  of  attack,  a,  for  the  baseline  M2.4-7A  configuration  at  a Mach 
number  of  2.40  and  Reynolds  numbers  varying  from  1 to  5 million  per  foot  for 
forced  b.l.  transition  using  trip  dots  of  nominal  height  0.0101”.  The  Ca  data 
above  show  that  at  Reynolds  numbers  of  1 and  1 .5  million/foot  the  curves  do 
not  look  smooth  (implying  that  the  boundary  layer  has  not  fully  transitioned). 
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The  above  plot  illustrates  that  the  derivative  of  Ca  (dCA/da)  for  the  Reynolds 
numbers  of  1 and  1 .5  million/foot  appear  wavy  implying  transitioning  of  the 
boundary  layer  is  occurring. 
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The  above  plot  utilizes  the  two  methods  of  calculating  trip  drag  described 
earlier  (either  calculating  a slope  for  Cd  at  a given  CL  versus  k or  k2)-  The 
estimated  drag  using  the  “k2”  method  tends  to  provide  values  of  CD  t|ipthat  are 
about  half  those  found  using  the  “k”  method.  Both  methods  predict  rather 
large  variations  with  differing  lifting  conditions. 
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High  Rn  Supersonic  Testing  I 

• Objectives 

- extension  of  HSR  experimental  database  to  higher  Rn  at 
supersonic  Mach 

- fully  turbulent,  free  transition  data  for  improved  trip  drag 
assessment  at  low  Rn 

• Background 

- existing  HSR  experimental  database  has  Rn.mac  max  of  approx.  7- 
8 million;  cruise  flight  Rn.mac  - 150-200  million 

- disagreement  in  HSR  community  as  to  best  trip  drag  assessment 
methodology;  partially  due  to  questions  of  boundary  layer  state 


In  order  to  further  understand  forced  transition  issues,  it  is  valuable  to  acquire 
test  data  with  naturally  occurring,  fully  turbulent  flow.  The  importance  lies  in 
the  ability  to  anchor  data  on  fully  turbulent  skin  friction  drag  predictions  which 
can  then  be  extrapolated  to  flight  conditions  without  the  uncertainty  of  trip 
drag.  On  the  other  side,  fully  turbulent  trip  free  data  provides  a target  for  low 
Rn  data  with  trip  drag  corrections. 

To  date,  the  existing  HSR  supersonic  experimental  database  has  an  Rn,mac 
maximum  in  the  range  of  7 to  8 million;  thecruise  flight  Rn,mac  is 
approximately  150  to  200  million.  The  high  Rn  supersonic  testing  effort  set  out 
to  extend  the  HSR  supersonic  data  base  to  high  Rn,  and  address  low  Rn  forced 
transition  methods  and  trip  drag  calculation  methods. 
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High  Rn  Supersonic  Testing  II 

• Approach 

- Test  free  & fixed  transition  configurations  from  low  to  high  Rn  in 
highest  Rn  supersonic  facility  available 

• free  transition  from  low  Rn  to  highest  Rn  attainable 

• fixed  transition  at  low  Rn  (variable  trip  heights,  types,  patterns) 

• apply  transition  detection  methods  to  assess  boundary  layer  state 

• MDA-E  Polysonic  Wind  Tunnel  has  highest  Rn  capability  in  the 
US 

- calibrated  at  M = 2.48  (close  to  design  cruise  Mach  for  HSCT) 

- 4 FT  test  section  (same  as  LaRC  UPWT) 

- blowdown  facility  (starting/ending  loads  are  an  issue) 


The  approach  taken  was  to  test  a representative  HSCT  configuration  with  both 
free  and  fixed  transition  from  low  Rn  to  the  highest  Rn  available  in  the  US.  It 
was  planned  to  evaluate  several  trip  types,  patterns,  and  locations,  and  several 
transition  detection  approaches. 

The  facility  in  the  US  with  the  highest  supersonic  Rn  capability  is  the  MDA-E 
Polysonic  Wind  Tunnel  (PSWT)  located  in  St.  Louis,  MO.  The  PSWT  has  a 4 
ft  test  section  allowing  existing  models  sized  for  the  LaRC  UPWT  to  be  tested. 
Start-up  & ending  loads  are  however  larger  than  experienced  at  the  UPWT  due 
to  the  fact  that  the  PSWT  is  a blowdown  facility.  At  present,  the  PSWT  is  not 
calibrated  for  Mach  = 2.4  (HSCT  nominal  cruise),  but  is  calibrated  at  Mach  = 
2.48.  For  the  purposes  of  this  investigation,  this  difference  is  not  significant. 
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High  Rn  Supersonic  Testing  III 

• Status 

- Test  1 (PSWT  689)  completed  November  25  1996 

- Model  Tested:  1 .675%  scale  M2.4-7A  baseline  wing/body 

- Rn.mac  from  8.3  to  20.7  million,  at  Mach  = 2.48 

- Good  shakedown  test 

• experience  gained  in  the  facility,  largest  model  tested  in  facility 

• HSR  cruise  data  quality  issues  addressed 

- questions  remain  about  flow  dew  point 

• good  balance/sting  combination;  currently  in  post-test  inspection 

• temperature  sensitive  paint  used  for  transition  detection 

• Plans 

- Complete  data  reduction  & analysis  from  test  1 

- 2nd  test  currently  scheduled  for  3QFY97 


The  model  chosen  for  testing  in  the  PSWT  was  the  M2.4-7A  baseline 
configuration  (same  model  as  used  in  UPWT  1667).  The  initial  test  in  the 
PSWT  was  conducted  in  November,  1996,  and  was  designated  PSWT689. 
Data  were  obtained  from  Rn,mac  = 8.3  to  20.7  million  (5.3  to  13.2  million  per 
foot)  at  Mach  = 2.48  over  a nominal  angle  of  attack  range  from  -3  to  + 3 deg. 

PSWT689  was  a good  shakedown  test.  Experience  was  gained  relative  to 
testing  in  the  facility  with  an  HSR  model,  and  data  quality  requirements.  The 
1.675%  scale  model  tested  was  the  largest  ever  tested  in  the  facility,  and 
accordingly,  significant  time  and  attention  was  given  to  monitoring  start-up 
and  ending  loads.  In  fact,  half  way  through  the  test  the  balance/sting 
combination  was  changed  to  a stronger  system.  The  second  system  consisted 
of  LaRC  balance  756  and  a 2000  lb.  normal  force  sting.  This  combination 
worked  well,  and  is  currently  in  post  test  inspection. 

TSP  was  used  for  transition  detection.  At  this  point,  the  results  are 
inconclusive.  Analysis  of  the  images  continues,  as  does  development  of  the 
system  for  application  in  a blowdown  facility. 

Force  data  from  the  test  is  to  be  reduced  again  to  correct  several  errors 
discovered  in  post  test  analysis.  It  should  be  noted  that  cavity  pressures  lagged 
during  the  1st  run  of  most  “blows”,  and  that  the  dew  point  was  not  known. 
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High  Rn  Supersonic  Testing  IV 

• Drag  plot  vs  Rn  using  available  data  & accounting  for  the  cavity 
pressure  lag 

• Forced  T ransition  runs  appear  to  follow  flat  plate  predicted 

*renC*'  M2.4-7A  Baseline  - M = 2.48 


The  above  chart  shows  minimum  drag  data  with  and  without  forced  transition 
as  compared  to  an  equivalent  flat  plate  prediction.  Experimental  data  is  that 
available  prior  to  a final  data  reduction  (as  discussed  previously),  but  has  been 
adjusted  manually  to  account  for  the  cavity  pressure  lag. 

Forced  transition  data  appear  to  follow  the  predicted  trend,  as  do  the  free 
transition  data  above  Rn/ft  - 8.5  million.  Note  also  that  the  trip  increment 
above  Rn/ft  ~ 8.5  million  is  on  the  order  of  2-3  drag  counts  which  is  similar  to 
that  determined  in  the  UPWT  tests  using  the  k2  trip  drag  assessment  method. 

This  analysis  will  be  updated  upon  final  data  reduction,  and  it  is  expected  to  be 
repeated  and  extended  in  a follow  on  test.  Balance  and  sting  limitations  should 
allow  extension  of  data  to  approximately  15  to  16  million  per  foot.  The 
limiting  factor  at  this  point  is  sting  divergence. 
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High  Rn  Subsonic^" ransonic  Testing  I 

• Objectives 

- develop  the  Temperature  Sensitive  Paint  (TSP)  test  technique 

- detection  of  tree  & fixed  transition  on  the  2.2  % HSR  Ref.  H wing 

- grit  drag  assessment  at  subsonic  & transonic  conditions  with  & 
without  flaps 

• Background 

- TSP  technique  identified  as  most  promising  approach  for  “global” 
transition  detection  in  a cryogenic  environment 

- joint  NASA/Industry/University  team  formed  to  develop  system 

- risk  reduction  experiments  executed  in  the  0.3m  TCT  prior  to  NTF 
test  to  evaluate  issues  such  as  paint  chemistry  & camera 
requirements 


Considerable  testing  of  the  HSR  Ref.  H model  has  been  conducted  in  the  NTF. 
The  primary  purpose  of  these  tests,  and  of  all  tests  in  the  NTF,  is  the 
assessment  of  Rn  effects.  Key  to  understanding  Rn  effects  is  the  understanding 
of  the  boundary  layer  state  at  all  conditions,  both  with  free  and  forced 
transition.  Several  methods  for  transition  detection  in  the  cryogenic 
environment  of  the  NTF  have  been  identified;  TSP  has  been  identified  as  the 
most  likely  candidate  to  be  successful  in  providing  a productive,  high-quality, 
global  assessment  of  boundary  layer  state. 

A joint  NASA/Industry/University  team  was  formed  to  develop  all  aspects  of 
such  a system.  Risk  reduction  experiments  were  conducted,  including  one  in 
the  0.3m  TCT  at  LaRC,  prior  to  application  in  the  NTF. 
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High  Rn  Subsonic/Transonic  Testing  II 

• Approach 

- Apply  the  TSP  system  in  a test  of  the  2.2%  Ref.  H model  in  the 
NTF 

- Evaluate  free  transition  boundary  layer  state 

- Confirm  boundary-layer  fixing  with  grit  trip 

- Assess  trip  drag 

- Realize  that  the  TSP  test  was  high  risk  with  potentially  high  pay-off 


Once  confidence  in  the  technique  was  gained  through  the  risk  reduction 
activities,  the  system  was  ready  for  the  NTF.  The  initial  test  was  still 
considered  to  be  high  risk  in  the  sense  that  successful  integration  of  a complex 
system  with  a complex  tunnel  would  be  difficult  on  the  first  try.  However, 
success  would  bring  high  pay-offs  in  terms  of  flow  field  understanding. 

The  first  part  focussed  on  applying  the  TSP  system  to  acquire  images  showing 
the  free  transition  boundary  layer  state.  In  addition,  grit  was  then  applied  and 
images  acquired  to  assess  the  trip  effectiveness.  Gaps  were  intentionally  placed 
in  the  trip  to  aid  in  the  assessment.  Several  results  are  shown  on  the  following 
charts.  Additional  analysis  is  presented  in  a paper  entitled:  “Use  of  Boundary 
Layer  Transition  Detection  to  Validate  full-scale  Flight  Performance 
Predictions,”  by  Hamner,  Owens,  Wahls,  and  Yeh  (presented  in  High  Lift 
session  of  this  workshop). 

The  second  phase  of  the  test  addressed  subsonic  and  low  transonic  trip  drag. 
Details  are  not  presented  here,  but  as  expected,  conventional  trip  location  and 
sizing  forced  transition  in  the  usual  manner  without  significant  trip  drag. 
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High  Rn  Subsonic  /Transonic  Testing  III 


The  chart  above  shows  several  typical  images  for  Mach  = 0.9  and  undeflected 
flaps  near  the  minimum  drag  condition.  Images  were  acquired  across  a wide 
temperature  range. 

As  noted  on  the  chart,  the  amount  of  laminar  flow  decreased  with  increasing 
Rn.  Although,  not  shown  here,  laminar  flow  did  exist  on  the  transonic  cruise 
configuration  at  Mach  = 0.9  and  the  cruise  angle  of  attack. 
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High  Rn  Subsonic/Transonic  Testing  IV 

• 10/3  outboard  flaps  w/  grit  strip 

• Mach  = 0.9,  Rn.mac  = 10.3  million,  a = 2 deg,  Tt  = 120°F 


upper  surface 


The  chart  above  demonstrates  the  effectiveness  of  the  conventional  grit  trip  at 
forcing  transition.  Laminar  flow  is  clearly  shown  progressing  through  the 
intentional  gaps  in  the  trip. 

The  TSP  transition  detection  technique  was  considered  very  successful.  Many 
lessons  relative  to  its  application  were  learned  and  work  on  the  technique  will 
continue. 
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Concluding  Remarks 

• Trip  dots  and  grit  performed  equally  well.  (Existing  sand  grit  at 
UPWT  was  mislabeled  and  not  well  “screened”  to  ensure 
consistent  size.) 

• Reynolds  number  sweeps  backed  up  with  flow  visualization 
techniques  worked  well. 

• Trip  drag  appears  to  be  best  calculated  by  Cd  (at  a given  CL) 
versus  k2  (However,  the  HSR  community  has  not  come  to 
consensus  on  this!) 

• Testing  in  high(er)  Reynolds  number  facilities  should  be  limited 
due  to  the  restrictions  caused  by  higher  loads,  shortness  of  run 
times,  cost,  etc. 

• Free-transitioning,  fully  turbulent  flow  appears  to  have  been 
reached  in  the  Polysonic  facility. 

• Temperature  Sensitive  Paint  for  transition  detection  is  viable 
and  worth  further  development. 


The  tests  showed  that  the  boundary  layer  could  be  forced  to  transition  equally 
well  using  either  grit  or  trip  dots.  The  trip  dots  still  appear  to  send  small 
vortices  downstream  through  the  boundary  layer  whereas  the  vortices 
generated  on  the  grit  seem  to  combine  and  not  travel  all  the  way  to  the  trailing 
edge  without  mixing  with  the  neighboring  disturbances). 

Very  large  differences  in  the  value  of  trip  drag  estimates  are  obtained  whether 
using  Cd  vs  k or  k2.  Follow  on  tests  in  the  Polysonic  facility  should  help 
understand  which  method  is  more  acceptable.  The  “k-method”  does  appear  to 
predict  levels  of  trip  drag  that  can  cause  fully  turbulent  drag  levels  that  would 
fall  below  free  transition  values  (found  experimentally)  implying  that  free 
transition  has  higher  drag  than  fully  turbulent  flow  which  is  clearly  incorrect. 

Future  testing  in  the  Polysonic  facility  should  produce  results  that  will  benefit 
all  supersonic  testing  done  in  the  HSR  program. 

Finally,  initial  application  of  temperature  sensitive  paint  for  transition 
detection  has  shown  it  to  be  a viable  test  technique,  and  worth  the  support  for 
further  system  development  to  improve  its  routine  use. 
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Pressure-Sensitive  Paint  and 
Video  Model  Deformation  Systems 
at  the  NASA  Langley 
Unitary  Plan  Wind  Tunnel 


G.  E.  Erickson 
A.W.Burner 
R.  DeLoach 


25  February  1997 
NASA  Langley  Research  Center 
Hampton,  Virginia 
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Pressure-Sensitive  Paint  and  Video  Model 
Deformation  Systems  at  the  NASA  Langley 
Unitary  Plan  Wind  Tunnel 

G.  E.  Erickson 
A.  W.  Burner 
R Deloach 

NASA  Langley  Research  Center 
Hampton,  Virginia  23681-0001 

Pressure-sensitive  paint  (PSP)  and  video  model  deformation  (VMD)  systems 
have  been  installed  in  the  Unitary  Plan  Wind  Tunnel  at  the  NASA  Langley 
Research  Center  to  support  the  supersonic  wind  tunnel  testing  requirements 
of  the  High  Speed  Research  (HSR)  program.  The  PSP  and  VMD  systems 
have  been  operational  since  early  1996  and  provide  the  capabilities  of 
measuring  global  surface  static  pressures  and  wing  local  twist  angles  and 
deflections  (bending).  These  techniques  have  been  successfully  applied  to 
several  HSR  wind  tunnel  models  for  wide  ranges  of  the  Mach  number, 
Reynolds  number,  and  angle  of  attack.  A review  of  the  UPWT  PSP  and  VMD 
systems  is  provided,  and  representative  results  obtained  on  selected  HSR 
models  are  shown.  A promising  technique  to  streamline  the  wind  tunnel 
testing  process,  Modern  Experimental  Design,  is  also  discussed  in  conjunction 
with  recently-completed  wing  deformation  measurements  at  UPWT. 
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• Review  of  UPWT  PSP  system  and  selected 
results  obtained  on  HSR  models 

• Review  of  UPWT  VMD  system  and  selected 
results  from  recently-concluded  HSR  testing 

• Discussion  of  a Modern  Experimental  Design 
method  for  improved  wind  tunnel  productivity 
--  recently  applied  at  UPWT  in  conjunction 

with  the  VMD  system  to  predict  HSR  model 
deformation 


A comprehensive  facility  enhancement  program  is  underway  at  the 
NASA  Langley  Research  Center  Unitary  Plan  Wind  Tunnel  to  provide 
state-of-the-art  test  techniques  to  support  the  supersonic  testing 
requirements  of  the  High-Speed  Research  program.  This  paper 
provides  a review  of  the  UPWT  pressure-sensitive  paint  system  for 
global  surface  static  pressure  measurements  and  the  UPWT  video 
model  deformation  system  to  measure  wing  local  twist  and  deflections. 
In  an  effort  to  improve  wind  tunnel  productivity,  a Modern 
Experimental  Design  technique  was  used  in  parallel  with  the  model 
deformation  system  to  determine  its  effectiveness  in  predicting  wing 
twist  at  supersonic  speeds.  Representative  results  obtained  with  the 
PSP,  VMD,  and  Modern  Experimental  Design  techniques  in  recent 
HSR  wind  tunnel  model  testing  at  UPWT  are  presented. 
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PSP  and  VMD  Systems  at  UPWT 


• PSP  and  VMD  systems  established  at  UPWT  in  1996  to  support 
HSR  experimental  programs 

• Strong  cooperative  efforts  involving  Aero-Gas  Dynamics  Division 
and  Experimental  Testing  Technology  Division  (ETTD) 

• PSP  system  evolved  from  previous  installation  in  the  NASA  Langley 
8-Foot  Transonic  Pressure  Tunnel  in  1994  in  cooperation  with  NASA 
Ames  and  ETTD 

• VMD  systems  established  at  the  National  Transonic  Facility  and 
Transonic  Dynamics  Tunnel  in  1994  provided  foundation  for  UPWT 
installation 

• PSP  and  VMD  systems  are  now  operated  by  UPWT  personnel  after 
extensive  training  with  ETTD  (setup,  operation,  image  acquisition 
and  processing,  and  data  analysis  and  plotting) 


Cooperative  efforts  involving  personnel  from  the  NASA  Langley  Aero-Gas 
Dynamics  Division  and  the  Experimental  Testing  Technology  Division  have 
led  to  the  establishment  of  PSP  and  VMD  deformation  systems  at  UPWT.  A 
PSP  system  previously  installed  in  the  NASA  Langley  8-Foot  Transonic 
Pressure  Tunnel  was  upgraded  and  installed  in  UPWT  in  early  1996.  This 
system  was  modeled  after  a similar  setup  currently  in  use  at  the  NASA  Ames 
Research  Center.  The  VMD  systems  established  by  ETTD  at  the  National 
Transonic  Facility  and  the  Transonic  Dynamics  Facility  served  as  models  for 
the  UPWT  installation,  which  was  initiated  in  1996.  Experience  gained  from 
several  HSR-sponsored  tests  has  provided  resident  expertise  at  UPWT  in  all 
aspects  of  the  PSP  and  VMD  systems. 
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UPWT  Description 

• Closed-circuit,  continuous-flow,  variable-density  tunnel 

• Two  4-ft  by  4-ft  test  sections 

• “Low  Mach”  test  section  has  a design  Mach  number 
range  from  1.5  to  2.9 

• “High  Mach”  test  section  has  a design  Mach  number 
range  from  2.3  to  4.6 

• Both  test  sections  use  asymmetric  sliding-block  nozzles 
that  allow  continous  variation  in  Mach  number 

• Maximum  Reynolds  number/foot  varies  from  6xl06to 
llxlO6,  depending  on  Mach  number 


The  NASA  Langley  UPWT  is  being  extensively  used  in  the  HSR  program  to 
provide  aerodynamic  performance  and  stability  and  control  characteristics  at 
supersonic  speeds.  The  ranges  of  Mach  number  and  Reynolds  number,  the 
control  of  the  dewpoint  and  stagnation  temperature,  the  optical  access  to  the 
test  section,  and  the  benign  environment  for  the  installation  of  digital  and 
video  imaging  equipment  are  factors  that  render  UPWT  well-suited  to  the 
application  of  PSP  and  VMD  techniques.  The  PSP  and  VMD  systems  are 
portable  between  the  “low  Mach”  and  “high  Mach”  number  test  sections.  In 
addition,  the  present  facility  scheduling  features  one  “active”  test  section 
while  the  “idle”  test  section  is  available  for  PSP  and  VMD  system  set-up, 
check-out,  and  test  technique  enhancements. 
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UPWT  PSP  System 


* Method  is  based  on  oxygen  sensitivity  of  photoluminescent  material 
in  the  form  of  a “paint”  (University  of  Washington  formulation) 

* PSP  sprayed  onto  model  surface  after  application  of  white  undercoat 

* Purge  air  is  applied  through  model  surface  pressure  orifices  during 
painting  process  via  electronic  pressure  scanner  modules  installed  in 
the  wind  tunnel  model 

* Registration  marks  applied  to  model  using  overlay  template 

9 PSP  excitation  source  is  2S0-W  lamps  that  emit  ultraviolet  light  in  a 
broadband  centered  around  360  nm  and  a cutoff  filter  to  block 
emission  in  the  visible  wavelengths 

9 All  other  illumination  sources  are  eliminated  by  installing  “light- 
tight” enclosures  on  both  sides  of  the  wind  tunnel  test  section 

The  UPWT  PSP  technique  is  based  on  a system  developed  by  the  NASA  Ames 
Research  Center  and  the  University  of  Washington.  The  photoluminescent 
paint  chemistries  developed  by  the  University  of  Washington  have  been  used  in 
all  of  the  NASA  Langley  UPWT  PSP  tests  and  have  proved  to  be  very  robust  for 
the  supersonic  experiments.  Approximately  one  shift  is  required  to  the 
application  and  curing  of  the  white  base  and  PSP  coatings.  Specially  modified 
ESP  scanners,  when  operated  in  a purge  mode,  route  air  through  the  wing 
surface  pressure  lines  to  prevent  the  paint  from  clogging  the  orifices. 

Registration  marks  are  applied  to  the  painted  surfaces  using  an  overlay 
template  and  a black  marking  pen.  The  marks  are  tyically  0..125  to  0.188  inches 
in  diameter  and  are  positioned  along  the  edges  of  the  model  and  at  selected 
other  locations  on  the  wing.  UV  lamps  are  mounted  to  the  webbing  of  the  the 
test  section  side  wall  to  provide  a continuous  illumination  source.  Manual 
shutters  are  used  to  block  the  UV  light  between  runs.  The  large  image  areas 
that  are  typically  mapped  on  HSR  models  requires  at  least  two  UV  lamps. 
Photodegradation  of  the  PSP  is  a concern  because  of  the  proximity  of  the  UV 
light  source  to  the  model  (approximately  2 feet),  so  double  and  triple  filters  are 
applied  to  the  lamps  to  reduce  the  UV  intensity  at  the  model.  A hand-held 
digital  radiometer  is  used  to  measure  the  effectiveness  of  the  UV  lamp 
positioning  and  filtering  arrangements.  Wooden  enclosures  have  been  built 
that  are  bolted/clamped  to  the  test  section  side  walls  to  eliminate  all  extraneous 
light  sources. 
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UPWT  PSP  System  (continued) 

• PSP  imaging  is  conducted  using  2 scientific-grade,  cooled  CCD 
digital  cameras  (12-bit  and  14-bit  intensity  resolution,  512x512 
and  1024x1024  spatial  resolution)  installed  in  the  “webbing” 
of  the  test  section  side  wall 

• Optical  filters  are  installed  on  the  camera  lenses  to  permit  the 
passage  of  the  luminescence  emission  wavelengths  while 
blocking  other  wavelengths 

• Model  is  rolled  90°  for  best  optical  access  and  the  model  pitch 
angle  is  varied  using  the  support  system  yaw  mechanism 

• Camera  integration  time  and  image  acquisition  are  controlled 
by  host  computers  located  in  the  UPWT  Data  Room,  about 
125  feet  from  the  camera  positions 


The  high  CCD  performance,  low  noise,  linear  response,  and  good  signal- 
to-noise  ratio  of  scientific-grade  digital  cameras  provide  high-precision, 
quantitative  light  measurements.  Two  digital  cameras  are  available  at 
UPWT  for  PSP  image  acquisition.  In  a typical  installation,  both  cameras 
are  mounted  inside  the  webbing  of  the  test  section  side  wall  with  lenses 
that  are  selected  to  provide  a detailed  view  of  an  area  of  particular 
interest  on  the  model  and  a global  view  of  the  wing  surface.  Special 
optical  filters  are  mounted  to  the  front  of  each  lens  so  that  the  camera 
detects  only  the  luminescence  emission  spectra.  Optical  access  to  the  test 
section  is  available  through  the  schlieren  windows  in  the  side  walls,  so  the 
model  is  rolled  90°  to  be  roughly  orthogonal  to  the  cameras.  Variation  in 
the  model  pitch  angle  is  obtained  using  a mechanized  yaw  mechanism. 
The  camera  exposure,  or  integration,  time  and  image  acquisition  are 
remotely  controlled  via  a Windows  95-based  PC  and  a UNIX-based 
workstation  that  are  located  in  the  UPWT  Data  Room.  The  integration 
times  are  selected  to  provide  high  image  intensity  while  avoiding  local 
saturation.  Typical  integration  times  are  1 to  1.5  seconds.  PSP  imaging 
has  not  been  compromised  by  the  relatively  minor  model  dynamics  that 
are  encountered  at  the  supersonic  speeds. 
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This  photograph  shows  a 1.675%-scale  HSR  arrow  wing  model  installed  in 
UPWT  test  section  1.  The  model  is  rolled  -90°  for  this  photograph,  although  the 
model  was  rolled  in  the  opposite  direction  for  the  testing.  The  right  (upper) 
wing  surface  is  coated  with  pressure-sensitive  paint,  while  the  left  (lower)  wing 
features  an  application  of  temperature-sensitive  paint. 
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The  photograph  shows  a close-up  view  of  a PSP  hardware  installation 
in  the  webbing  of  UPWT  test  section  2.  Two  scientific-grade  CCD 
digital  cameras,  one  standard  video  CCD  camera,  and  three  250W  UV 
lamps  are  installed  using  articulated  mounting  arms  or  C-clamps.  The 
webbing  provides  a stable,  virtually  vibration-free  mounting  surface 
for  the  PSP  imagers  and  light  sources. 
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The  PSP  digital  camera  electronic  control  units  and  chiller  units  are  shown 
in  this  figure.  Excess  length  of  a 200-foot  fiber  optic  cable  is  shown  at  the 
bottom  of  the  mobile  cart,  which  is  positioned  adjacent  to  the  test  section. 


518 


UPWT  PSP  System  (concluded) 

• Personal  computer  controls  the  14>bit  camera  via  a proprietary 
interface  card  and  electronics  cable,  while  a UNIX-based 
workstation  controls  the  12-bit  camera  using  a separate  interface 
and  a fiber  optic-based  SCSI  link 

* Video  camera  with  optical  filter  provides  real-time  PSP  response 

9 Wind-off  and  wind-on  images  are  processed  on  the  UNIX  machine 
using  the  NASA  Ames-developed  “paintcp”  software  package, 
which  performs  the  image  ratio,  registration,  and  paint  calibration 
operations 

9 PSP  is  calibrated  via  an  “in-situ”  method  using  surface  static 
pressures  obtained  from  discrete  orifices  on  the  model  surface 
connected  to  internally-mounted  ESP  scanners 

9 Image  processing,  analysis,  and  plotting  are  performed  on-site 
and  results  posted  on  WWW  site  established  for  each  UPWT  test 

The  14-bit  digital  camera  with  512x512  pixel  array  is  controlled  by  a personal 
computer  and  camera  interface  card.  A camera  electronics  cable  extends  from 
the  interface  card  to  an  electronics  control  unit  and  camera  chiller  unit  assembly 
located  adjacent  to  the  test  section.  The  12-bit  digital  camera  with  1024x1024 
pixel  array  is  controlled  by  a high-end  workstation  that  features  a fiber  optic- 
based  SCSI  bus  extender  from  the  workstation  to  the  camera  control  unit/chiller 
unit  assembly,  also  positioned  in  proximity  to  the  test  section.  Electronics  cables 
and  fiber  optics  cables  are  permanently  routed  from  the  Data  Room  to  both  test 
sections.  A separate  video  CCD  camera  with  optical  filter  is  mounted  to  the  test 
section  webbing  to  provide  real-time  display  and  recording  of  the  paint  response, 
which  can  include  the  footprints,  or  signatures,  of  shock  waves  and  vortices.  An 
extensive  disk  array  has  been  assembled  to  accommodate  the  image  storage 
requirements  of  PSP  testing.  All  images  are  transferred  to  the  workstation, 
where  the  image  ratioing,  image  registration,  and  PSP  calibrations  are 
performed  using  a software  package  developed  by  NASA  Ames  Research  Center. 
An  “in-situ”  calibration  is  performed  whereby  the  paint  is  calibrated  using  the 
static  pressures  measured  at  discrete  locations  with  internally-mounted  pressure 
scanners.  All  image  processing  operations  and  data  analysis  and  plotting  are 
done  on-site.  World  Wide  Web  sites  are  typically  established  for  each  UPWT  test 
to  allow  posting  of  the  processed  PSP  images. 
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The  host  computers  that  control  the  two  PSP  digital  cameras  are  shown  in  the 
photograph  above.  The  PC  (Windows  95  OS)  and  UNIX  workstation  are  situated 
side-by-side  along  with  high-capacity  disk  drives  and  color  postscript  printer. 
Additional  magneto-optical  hard  drives  and  recordable  CD  drives  have  recently 
been  acquired  to  augment  the  UPWT  PSP  system. 
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UPWT  PSP  Applications  to  Date 


• 1.675%-scale  HSR  arrow  wing  model 

• 1.675%-scale  HSR  TCA  2a  model 

• Test  sections  1 and  2 

• M=1.6  to  2.7,  Re/ft=3,  4 million 

• a = -2°  to  8° 

• Attached  flow,  separated  (vortex)  flow,  shock  waves 


The  UPWT  PSP  system  has  been  applied  to  several  HSR  models,  including  a 
1.675%-scale  HSR  arrow  wing  model  in  test  sections  1 and  2 and  a 1.675%- 
scale  model  of  the  HSR  TCA2a  model  in  test  section  2.  PSP  results  have 
been  obtained  for  a wide  range  of  Mach  number  and  angle  of  attack  that 
encompass  flow  regimes  dominated  by  attached  flow,  vortices,  and  shock 
waves.  Time  required  to  set  up  the  PSP  system  and  acquire  all  wind-on  and 
wind-off  images  is  approximately  2 shifts.  Additional  time  is  required  at  the 
outset  of  the  wind  tunnel  entry  to  acquire  flow  angle  corrections  (upflow  and 
sideflow)  to  provide  accurate  determination  of  the  model  angle  of  attack. 
Runs  are  also  made  of  the  unpainted  model  at  the  same  test  conditions  to 
quantify  any  obtrusive  effects  of  the  paint  thickness  on  the  wing  surfaces. 
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Arrow  Wing  PSP  Image  at  M=1.65,  a=8° 

UPWT  Test  1836  Re/ft  = 3(10*)  June  1996 


The  photograph  shown  above  is  a false-colored,  ratioed  and  registered  image 
of  the  wing  upper  surface  pressure  field  on  a 1.675%-scale  model  of  an  HSR 
arrow  wing  configuration.  The  test  conditions  correspond  to  a free-stream 
Mach  number  of  1.65,  Reynolds  number  per  foot  of  3(106),  and  angle  of 
attack  of  8°.  Free-stream  stagnation  temperature  is  125  °F.  The  PSP  image 
clearly  shows  the  signatures  of  leading-edge  vortices  that  develop  from  the 
inboard  and  outboard  wing  regions.  The  inboard  wing  vortex  passes  over  the 
outboard  nacelles,  and  the  effect  of  the  vortex  passage  can  be  correlated  with 
the  nacelle  base  pressure  measurements.  The  paint  was  calibrated  using 
pressure  measurements  obtained  at  discrete  ports  with  two  ESP  modules. 
This  test  was  conducted  in  June  1996  and  was  the  first  application  of  the 
UPWT  PSP  system  in  UPWT  test  section  1. 
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A processed  PSP  image  of  the  wing  lower  surface  on  a 1.675%-scale  model 
of  the  HSR  TCA  2a  in  UPWT  test  section  2 is  shown  above.  The  test 
conditions  correspond  to  a Mach  number  of  2.4,  Re/ft  = 4(106),  Tstag=125°, 
and  a = 3.5°.  The  inboard  and  outboard  nacelles  were  installed  for  this 
run,  but  were  painted  flat  black  to  eliminate  the  effect  of  reflected  light 
from  the  sides  of  the  nacelles  on  to  the  wing  surface.  Of  particular  interest 
in  this  application  was  the  character  of  the  interacting  shock  waves 
developed  by  the  nacelle  diverters.  Reflected  shocks  from  the  diverters  are 
also  discernible  in  the  original  image  and  in  the  PSP-derived  static  pressure 
distributions.  Several  streamwise  rows  of  wing  lower  surface  static 
pressure  orifices  were  plumbed  to  an  ESP  module  without  the  purge  air 
capability.  As  a result,  thin  strips  of  masking  tape  were  applied  to  these 
rows  during  the  painting  process.  These  unpainted  strips  are  visible  in  the 
image  above.  The  in-situ  paint  calibration  required  the  selection  of  pixel 
locations  outside  of  these  regions. 
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Pressure-  Sensitive  Paint  Results  on  HSR  TCA  2 Wing  Lower  Surface 
M=2.40  (X=3.5°  Re/ft=4(106)  UPWT  Test  1671  (PXL  Camera) 


This  composite  plot  compares  the  streamwise  distributions  of  the  wing  lower 
surface  static  pressure  coefficient  obtained  with  the  PSP  technique  and  the 
electronically-scanned  pressure  modules  at  Mach=2.4,  Re/ft=4(106),  and  a = 
3.5°  (same  case  as  on  previous  page).  The  ratioed  and  registered  PSP  image 
and  a model  installation  image  are  also  included.  The  PSP  and  ESP  pressure 
data  compare  very  well,  and  the  maximum  difference  in  the  coefficients 
obtained  with  the  two  methods  is  within  approximately  5%  of  full-scale  range 
of  the  ESP  transducers  . The  abrupt  pressure  rise  across  the  shocks  is 
apparent  in  the  first  four  pressure  distribution  plots.  Note  that  the  PSP  data 
plots  are  restricted  to  values  obtained  at  a single  pixel  in  proximity  to  each 
pressure  orifice.  The  advantage  of  the  PSP  method  is  that  all  image  pixels 
are  “pressure  tap”  locations,  and  the  corresponding  hundreds  of  thousands  of 
pixels  (depending  on  the  camera  resolution)  can  provide  much  higher 
resolution  of  the  CP  distribution,  particularly  across  the  shock  waves. 


UPWT  VMD  System 

• VMD  technique  is  based  on  a single  video  camera  photogram  metric 
determination  of  two-dimensional  coordinates  of  wing  targets  with 
a known  fixed  third  dimensional  coordinate  (spanwise  location) 

• Primary  application  of  UPWT  VMD  system  is  to  determine  local 
wing  twist,  while  secondary  applications  include  wing  deflection 
(bending)  and  model  angle  of  attack  measurements 

• Retroreflective  dots  with  adhesive  backing  are  applied  in  several 
chordwise  rows  from  the  wing  root  to  the  wing  tip  to  provide  high- 
contrast  targets 

• Images  are  acquired  using  a standard  RS-170  CCD  video  camera 
with  752  horizontal  by  240  vertical  pixel  resolution 

• Illumination  source  is  a fiber  optic-based  ring  light  mounted  to  the 
front  of  the  camera’s  10  to  100mm  focal  length  remote  zoom  lens 

A unique  aspect  of  the  video  model  deformation  system  developed  by  NASA 
Langley  ETTD  is  the  photogrammetric  determination  of  two-dimensional  wing 
targets  using  a single  video  camera.  A requirement  is  that  the  third  dimensional 
coordinate  be  known  and  fixed,  namely,  the  spanwise  location  of  the  targets.  The 
primary  application  of  the  VMD  system  is  to  measure  the  wing  local  twist  angle, 
although  the  wing  deflection  (bending)  and  model  angle  of  attack  measurements 
may  be  equally  significant  depending  on  the  experimental  objectives.  Targets  in  the 
form  of  retroreflective  dots  with  an  adhesive  backing  are  applied  at  precisely 
known  locations  in  chordwise  rows  at  several  wing  span  stations.  The  inboard  row 
of  targets  is  placed  in  a region  of  the  wing  that  may  be  considered  rigid.  This  row 
serves  as  a reference  to  all  other  target  rows  and  provides  an  “onboard”  angle  of 
attack  measurement.  The  dots  provide  extremely  high-contrast  targets  for  image 
acquisition,  and  any  glints  or  other  undesired  sources  of  high  contrast  on  the  model 
surface  are  eliminated  by  applying  a thin  coat  of  flat  black  paint  to  these  regions 
(Note:  This  is  not  a standard  practice  at  all  facilities.)  The  thickness  of  the  targets 
is  somewhat  instrusive  and  may  cause  drag  coefficient  increments  of  a few  counts  at 
the  supersonic  speeds;  as  a result,  the  VMD  measurements  are  generally  made  in  a 
separate  run  series  in  a manner  similar  to  the  PSP  technique.  A standard  CCD 
video  camera  with  characteristics  that  have  been  well-documented  by  ETTD  is  used 
to  acquire  images  of  the  targeted  region.  Uniform  illumination  of  the  model  is 
provided  by  a fiber  optic-based  ring  light  that  easily  attaches  to  the  front  of  the 
camera’s  remote  zoom  lens. 
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The  1.5%-scale  HSR  TCA  20  model  installed  in  UPWT  test  section  1 is  shown 
in  the  above  photograph,  taken  at  the  conclusion  of  a recent  video  model 
deformation  experiment.  The  5 chordwise  rows  of  retroreflective  targets  are 
visible  on  the  right  wing  upper  surface. 
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The  figure  above  shows  a close-up  view  of  the  right-hand  wing  upper  surface  on 
the  1.5  /o-scale  HSR  TCA  model  20  installed  in  UPWT  test  section  1.  Five  rows 
of  retroreflective  targets  are  visible;  the  first  four  chordwise  rows  (starting  from 
the  wing  root  region)  feature  four  0.188-inch  diameter  targets  while  the  fifth 
row  (at  the  wing  tip)  has  three  0.125-inch  diameter  targets.  The  wing  twist  and 
deflection  measurements  that  are  presented  in  following  figures  correspond  to 
the  row  near  the  wing  tip.  Note  that  several  smaller  screw  holes  in  the  wing 
surface  are  filled  with  dental  plaster,  and  these  holes  can  be  miscontrued  in  the 
photograph  as  VMD  targets.  These  holes  are  not  visible  during  the  image 
acquisition  process. 
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A close-up  view  of  the  VMD  system  camera  installed  in  the  webbing  of  the 
test  section  is  shown  in  this  photograph.  The  standard  video  CCD  camera  is 
mounted  to  the  remote  zoom  lens  which,  in  turn  is  bolted  to  an  angle  plate 
that  is  C-clamped  to  the  webbing.  The  fiber  optic  link  to  the  lens-mounted 
ring  light  is  also  discernible  in  the  figure.  Considerable  care  is  required  to 
ensure  that  that  the  focal  length  and  camera  position  are  not  changed  once 
the  camera  calibration  is  completed.  The  video  signal  from  the  camera  is 
routed  to  the  Data  Room  via  an  RG-59  cable  to  a video  distribution  amp  and 
to  the  video  framegrabber  board  in  the  PC  image  acquisition  system.  Set-up 
of  the  VMD  camera  equipment  is  very  straightforward  and  requires  less  time 
than  the  PSP  hardware  installation. 
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UPWT  VMD  System  (concluded) 


• The  video  signal  is  routed  to  a frame  grabber  controlled  by  a 
120  MHz  Pentium  PC  in  the  UPWT  Data  Room 

• Detailed  in-tunnel  static  calibrations  are  performed  using  a target 
plate  to  determine  the  camera  position  and  pointing  angles 

• Wind-off  pitch  sweeps  of  the  model  (in  the  upright  position)  and 
retroreflective  targets  installed  are  then  conducted 

• Automated  system  analyzes  several  digitized  video  images  at  each 
angle  of  attack  during  the  wind-off  and  wind-on  pitch  sweeps  and 
displays  “raw”  values  of  the  wing  local  angle  of  attack  and  vertical 
(“z”)  coordinates 

• Commercially-available  numeric  computation  and  visualization 
software  package  is  used  to  compute  and  plot  final  wing  twist  and 
wing  bending  results 


Images  are  acquired  from  the  video  CCD  camera  using  a frame  grabber 
board  installed  in  a Windows  95-based  PC.  Acquisition  of  digitized  video 
images  is  triggered  by  a “pickle  switch”  or  a keyboard  command,  and  the 
automated  system  identifies  the  model  targets  and  analyzes  several  video 
images  at  each  angle  of  attack.  Tunnel  test  condition  information  is  also 
acquired  at  this  time  from  the  wind  tunnel  data  acquisition  system  via  an  RS- 
232  interface.  The  tunnel  test  conditions,  test  point  information,  and  the 
values  of  the  uncorrected  local  angle  of  attack  and  vertical  displacement  at  all 
target  rows  are  then  displayed  on-screen,  at  which  point  the  system  is  ready 
for  the  next  data  point.  Target  plate  calibrations  and  wind-off  and  wind-on 
data  are  acquired  in  the  same  manner.  The  calibration  of  the  camera  is  a 
detailed  procedure  which  uses  a target  plate  rig  constructed  for  the  UPWT 
system  and  yields  the  camera  location,  pointing  angles,  and  effective  focal 
length.  Wind-off  pitch  sweeps  of  the  model  in  the  upright  position  are 
conducted  to  verify  the  camera  calibration  and  to  provide  static  “tares”  that 
are  subtracted  from  the  wind-on  data.  Post-run  processing  of  the  VMD  data 
is  done  using  a commercially-available  software  package  that  computes  and 
plots  the  corrected  wing  twist  and  deflection  results. 
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This  is  a close-up  view  of  the  1.5%-scale  HSR  TCA  20  model  with  a VMD 
system  calibration  target  plate  placed  on  the  right  wing  upper  surface.  The 
target  plate  consists  of  49  targets  with  precisely  known  x-  and  y-coordinates 
measured  on  a NASA  Quality  Assurance  validator.  In  practice,  the  target 
plate  is  mounted  to  a platform  that  has  precise  control  of  the  y-  and  z-position 
of  the  plate  relative  to  the  model  centerline.  The  target  plate  is  a critical 
element  in  the  determination  of  the  VMD  system  camera  position  and 
pointing  angles. 
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The  calibration  rig  that  is  used  to  determine  the  camera  constants  (location, 
pointing  angles,  effective  focal  length)  is  shown  in  this  photograph.  The  target 
plate  is  situated  atop  the  rig  and  slides  over  the  top  of  the  wing  surface.  The  x- 
position  of  the  calibration  assembly  is  set  to  provide  a satisfactory  image  from 
the  video  CCD  camera,  and  the  y-  and  z-positions  of  the  plate  are  varied  using 
optical  rail  and  lab  jack  arrangements.  The  y-  and  z-displacements  of  the 
target  plate  are  measured  using  dial  gauges  (only  one  gauge  is  shown  installed 
in  the  present  photograph).  As  the  calibration  progresses,  the  calibration  rig 
displacements  measured  with  the  dial  gauges  are  compared  to  similar 
measurements  obtained  with  the  VMD  system. 
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The  host  computer  for  the  VMD  system  is  illustrated  in  the  above 
photograph.  The  mini-tower  case  contains  the  video  frame  grabber  board 
that  acquires,  stores,  and  analyzes  the  digitized  images.  All  image 
acquisition,  processing,  analysis,  and  plotting  of  the  VMD  results  can  be 
performed  from  this  site. 
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VMD  Applications  to  Date 


• 1.675%-scale  HSR  Reference  H model 

• 1.675%-scale  HSR  TCA  2a  model 

• 1.675%-scale  HSR  TCA  20  model 

• Test  sections  1 and  2 

• M=1.6  to  2.7,  Re/ft=l-5  million 

• a = -4°  to  12° 


An  early  proof-of-concept  test  of  the  VMD  system  applied  to  a 1.675%- 
scale  HSR  Reference  H model  focused  on  the  measurement  of  the  wing 
twist  at  supersonic  speeds.  Each  subsequent  test  led  to  enhancements  of 
the  UPWT  VMD  system.  Primary  improvements  include  additional 
target  rows  across  the  wing  span  to  provide  wing  twist,  deflection,  and 
secondary  model  pitch  angle  measurements;  improved  method  of 
“spatially  mapping”  the  wind-off  and  wind-on  results;  and  development 
of  an  effective  target  rig  that  significantly  streamlined  the  camera 
calibration  process.  UPWT  provides  an  excellent  environment  for  this 
test  technique,  which  has  been  successfully  applied  in  both  test  sections 
over  wide  ranges  of  Mach  number,  Reynolds  number,  and  angle  of 
attack. 
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Mach  Number  Effect  on  Wing  Twist 
HSR  TCA  20,  Re/ft=3(106),  UPWT  Test  1844 


a,  deg 

FEDD  R««triel*d  • P itlminary  Dan  Modal  dalomtation  data.  UPWT  Ta«l  1844. 


The  effect  of  the  free-stream  Mach  number  on  the  wing  twist  at  the 
y/(b/2)=0.989  span  station  is  presented  above.  The  Mach  number  varies  from 
1.60  to  2.70  at  a constant  Reynolds  number  per  foot  of  3(106).  Increasing  the 
Mach  number  decreases  the  wing  twist  (washout)  at  a given  angle  of  attack. 
This  effect  is  caused  by  the  reduced  wing  lift  as  the  Mach  number  is  increased. 
A maximum  twist  angle  of  approximatly  -2.9°  was  obtained  at  Mach=1.60  and 
a = 12°;  at  the  same  angle  of  attack,  the  twist  angle  was  about  -1.2°  at 
Mach=2.70. 
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Mach  Number  Effect  on  Wing  Deflection 
HSR  TCA  20,  Re/ft=3(106),  UPWT  Test  1844 


0 Ot,  deg 

FE  D D R«uacii£)  P ralim  inary  Data  M o4el  dtform ation  data.  U P WT  Ttct  1844 


0 

0 


The  Mach  number  effect  on  the  wing  deflection  (bending)  at  the  y/(b/2)=0.989  span 
station  is  shown  in  this  figure.  Increasing  the  Mach  number  decreases  the 
deflection  in  the  z (vertical)  axis  at  a given  angle  of  attack  (less  upward  bending  at 
the  wing  tip).  The  maximum  deflection  of  approximately  031  inches  was  obtained 
at  Mach =1.60  and  a = 12°;  the  z displacement  was  0.15  inches  at  Mach=2.7  and  the 
same  angle  of  attack. 
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Reynolds  Number  Effect  on  Wing  Twist 
HSRTCA20,  Mach=2.10,  UPWT  Test  1844 
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The  Reynolds  number  effect  on  wing  twist  at  the  y/(b/2)  = 0.989  span  station 
and  a constant  Mach  number  of  2.10  is  shown  in  the  data  plot  above.  The 
Reynolds  number  per  foot  varies  from  1.0  (106)  to  5.0(10*)  in  increments  of 
1(106).  The  trend  in  the  data  plot  is  more  of  a free-stream  dynamic  pressure 
(“q”)  effect  than  Reynolds  number,  since  q varied  from  approximately  221 
psf  at  Re/ft  = 1(106)  to  1100  psf  at  Re/ft  = 5(10*).  The  twist  angle  is 
approximately  a linear  function  of  the  Reynolds  number  (“q” );  for  example, 
a five-fold  increase  in  the  Reynolds  causes  a corresponding  increase  in  the 
twist  angle  near  the  wing  tip.  For  the  range  of  angle  of  attack  tested,  a 
maximum  twist  of -2.75°  occurs  at  a Reynolds  number  of  5(106)/ft  and  a = 
12°,  while  the  corresponding  twist  at  Re/ft=  1(10*)  is  -0.55°. 
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Reynolds  Number  Effect  on  Wing  Deflection 
HSRTCA20,  Mach=2.10,  UPWT  Test  1844 
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Wing  deflection  measurements  obtained  during  a Reynolds  number  “sweep” 
at  constant  Mach  number  (Mach=2.10)  are  shown  in  this  data  plot.  The 
Reynolds  number  per  foot  varies  from  1.0  (106)  to  5.0(10*)  in  increments  of 
1(106).  The  z-displacement  is  approximately  a linear  function  of  the 
Reynolds  number.  Similar  to  the  results  shown  in  the  previous  figure,  the 
primary  factor  affecting  the  wing  displacement  is  the  free-stream  dynamic 
pressure.  For  the  a -range  in  the  present  test,  a maximum  deflection  of  0.36 
inches  occurs  at  a Reynolds  number  of  5 (106)/ft  and  a = 12®,  while  the 
corresponding  displacement  at  Re/ft=  1(106)  is  0.10  inches. 


FE  0 0 Rtftritrfcd  • P reim  inary  Data 
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The  term  "Classical  Design"  is  used  to  describe  an  approach  to  experimentation 
in  which  one  variable  is  changed  at  a time  while  all  other  variables  are  held 
constant.  Classical  designs  have  been  used  in  wind  tunnel  research  at  Langley 
since  the  earliest  days  of  flight,  and  are  widely  used  in  wind  tunnel  testing 
elsewhere  as  well. 

Today,  important  aircraft  design  decisions  can  turn  on  fractional  drag  count 
results,  and  practitioners  of  an  alternative  experiment  design  philosophy  called 
"Modern  Design"  recognize  the  futility  of  "holding  everything  constant"  which 
might  affect  results  at  this  level.  Instead,  they  exploit  their  knowledge  of  the 
stochastic  nature  of  experimental  variables  to  control  error  through  balance  and 
randomization.  Modern  and  classical  design  philosophies  also  differ  in  their 
approach  to  productivity  enhancement.  Classical  designers  attempt  to  maximize  • 
the  data  volume  for  a given  resource  budget  while  modem  designers  attempt  to 
achieve  a specific  technical  result  at  a prescribed  level  of  confidence  with  the 
smallest  expenditure  of  resources  possible. 
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Classical  designs  divide  a given  inference  space  into  a "grid"  or  "matrix"  of  test 
conditions  at  which  response  variables  of  interest  are  measured  (forces, 
moments,  etc.)  The  extent  of  this  grid  and  the  size  of  the  cells  which  comprise  it 
are  influenced  by  the  amount  of  resources  available  for  a given  test. 

Modern  design  practitioners  use  the  concept  of  a "response  surface"  to  guide 
their  design  efforts.  A response  surface  is  a logical  extension  of  the  simple  one- 
variable  line  graph  in  which  the  dependence  of  the  response  variable  on  all 
relevant  independent  variables  is  simultaneously  considered  in  a small  region  of 
interest  in  the  inference  space.  The  extent  of  this  region  is  purposely  limited  to 
that  in  which  the  response  variable  can  be  approximated  adequately  by  a low- 
order  Taylor  series.  Methods  such  as  regression  and  contrast  analysis  are  used 
to  elucidate  the  response  surface.  Various  curvature  tests  and  optimization 
procedures  are  then  used  to  quickly  identify  regions  in  the  inference  space 
which  are  the  most  interesting  (peaks,  ridge  systems,  etc.),  which  reduces 
resources  that  would  otherwise  be  spent  in  less  profitable  regions. 
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The  aerodynamically-induced  increase  in  wing  twist  for  an  HSR  stability  and 
control  model  has  been  measured  for  a range  of  angles  of  attack,  Mach 
numbers,  and  Reynolds  numbers  as  described  elsewhere  in  this  paper.  A 
classical  design  requiring  330  data  points  was  initially  conducted,  followed  by  a 
modern  design  to  likewise  quantify  the  wing  twist  change  for  the  same  model. 
The  modem  design  required  only  20  data  points  to  define  wing  twist  as  a 
second-order  response  function  in  3 variables  with  a design-goal  precision  of 
0.05°  at  a prescribed  95%  confidence  level,  given  the  0.04°  standard  deviation 
in  measured  wing  twist  that  was  anticipated.  The  figure  above  compares  95% 
confidence  intervals  for  the  classical  and  modem  designs.  Both  methods 
generated  results  with  a precision  that  met  the  0.05°  design  goal. 
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The  two  solid  curves  in  the  above  figure  mark  the  upper  and  lower  limits  of  the 
95%  prediction  interval  for  the  modern  design  results  at  a given  Mach  number, 
Reynolds  number,  and  normalized  semi-span  location.  This  modern  design 
prediction  was  confirmed  by  plotting  the  33  data  points  acquired  on  a different 
day  at  the  same  conditions,  using  the  classical  pitch-sweep  method.  Similar 
results  were  obtained  at  other  combinations  of  Mach  and  Reynolds  number. 

Note  that  the  above  combination  of  Mach  number  and  Reynolds  number  was 
never  actually  run  in  the  modern  design.  This  figure  simply  represents  a slice 
through  the  modern  design  response  surface  in  a direction  parallel  with  the 
"angle  of  attack  axis"  at  the  specified  values  of  Mach  number  and  Reynolds 
number.  This  illustrates  the  fact  that  modem  design  response  surface  methods, 
once  the  response  surface  is  adequately  defined  the  response  can  be  quantified 
for  other  combinations  of  the  independent  variables  besides  those  measured 
directly. 
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The  modem  design  method  only  requires  enough  data  to  fit  a low-order  (typically 
first  or  second  order)  function  of  the  independent  variables  in  regions  of  interest, 
plus  sufficient  additional  data  points  to  insure  that  design  precision  goals  are 
met  with  a prescribed  level  of  confidence.  A Central  Composite  Rotatable 
Design  (CCRD)  was  employed  in  this  test  which  could  accomplish  these 
objectives  with  only  20  data  points.  This  resulted  in  considerably  fewer  wind-on 
minutes  than  the  classical  design  (approximately  one  third  in  this  test.) 

Additional  comparison  tests  involving  other  response  variables,  other 
independent  variables  and  different  ranges  of  variables,  and  other  facilities, 
must  be  conducted  before  a body  of  practical  tunnel-testing  experience  will  have 
been  accumulated  which  is  sufficient  to  warrant  a general  implementation 
recommendation.  However,  modern  design  methods  have  been  shown  in  this 
test  to  have  some  promising  potential  for  wind  tunnel  research  in  an  era  in 
which  external  pressures  continue  to  dictate  that  more  be  accomplished  with 
less. 
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Summary 


• PSP  and  VMD  systems  are  installed  and  operational  at  UPWT 

• Test  techniques  provide  global  surface  pressure  mapping, 
qualitative  surface  flow  visualization,  and  model  deformation 
measurements  (twist  and  bending)  at  supersonic  speeds 

• Expertise  has  been  developed  at  UPWT  that  allows  autonomous 
operation  of  both  systems 

• PSP  and  VMD  systems  are  “works  in  progress”  that  will  be 
subject  to  continued  enhancements 

• Modern  experimental  design  technique  was  effective  in  capturing 
wing  twist  characteristics  and  may  provide  a means  of 
streamlining  the  wind  tunnel  test  process 

Pressure-sensitive  paint  and  video  model  deformation  systems  are 
installed  in  the  NASA  Langley  Unitary  Plan  Wind  Tunnel  and  have  been 
operational  since  early  1996.  The  PSP  and  VMD  systems  has  been 
effectively  used  in  support  of  HSR  supersonic  wind  tunnel  testing  to 
provide  global  surface  pressure  mapping,  qualitative  surface  pressure 
field  response  to  shock  waves  and  vortex  flows,  and  measurements  of  the 
wing  local  twist  angle  and  deflections  (bending).  Time  to  set  up  and 
calibrate  the  PSP  and  VMD  systems  is  one  shift  (each),  while  one  shift  for 
each  technique  is  necessary  to  acquire  a typical  set  of  wind-on  runs. 
Simultaneous  installation  of  the  PSP  and  VMD  systems  has  been  done, 
although  the  images  from  each  system  were  acquired  in  a concurrent, 
rather  than  simultaneous,  manner  because  of  system  conflicts.  Future 
enhancements  to  these  systems  that  may  lead  to  a “turn-key”  operation 
include  the  ability  to  remotely  control  all  illumination  sources,  including 
mechanized  shutters  for  the  UV  lamps  to  reduce  the  effects  of 
photodegradation,  and  full  automation  of  the  image  acquisition  process. 
The  experience  gained  from  PSP  and  VMD  testing  in  cooperation  with 
ETTD  has  resulted  in  resident  expertise  at  UPWT  regarding  virtually  all 
aspects  of  the  system  operations  (application  of  the  PSP  coating  continues 
to  be  performed  by  ETTD).  A Modern  Experimental  Design  technique 
was  used  during  a recent  VMD  test  where  all  critical  test  parameters  were 
varied  at  each  of  20  data  points.  The  resultant  response  surface  proved 
effective  in  predicting  the  wing  twist  over  ranges  of  Mach  number, 
Reynolds  number,  and  angle  of  attack. 
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Analysis  and  Multipoint  Design 
of  the  TCA  Concept 


Steven  E.  Krist 
Steven  X.  S.  Bauer 


Analysis  and  Multipoint  Design  of  the  TCA  Concept 

The  goal  in  this  effort  is  to  analyze  the  baseline  TCA  concept  at  transonic  and  supersonic  cruise,  then  apply  the 
natural  flow  wing  design  concept  to  obtain  multipoint  performance  improvemnts.  Analyses  are  conducted  with 
OVERFLOW,  a Navier-Stokes  code  for  overset  grids,  using  PEGSUS  to  compute  the  interpolations  between 
the  overset  grids. 
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Grid  Generation  and  Codes 


The  initial  20  block  overset  grid  for  the  TCA  baseline 
configruation,  Mach  2.4  cruise  condition,  was  developed  by 
Boeing 

Slight  topology  changes  were  made  to  the  nacelles  to  permit  the 
use  of  PEGSUS41_45  (7  times  faster  than  PEGSUS41_36) 

Grids  for  transonic  cruise  at  Mach  0.95 

- 1/3  smaller  wall-normal  spacing 

- 10%  more  grid  points 

OVERFLOW  version  1 .7t 


Grid  Generation  and  Codes 

The  initial  overset  grids  utilized  in  this  effort  were  developed  by  Steve  Chaney  and  Steve  Ogg  at  Boeing.  The 
grids  were  sized  for  Mach  2.4  supersonic  cruise  conditions,  maintaining  a y+  value  on  the  order  of  one.  The 
grids  were  modified  for  the  Mach  0.95  transonic  cruise  condition  by  decreasing  the  wall-normal  spacing  to  a 
third  of  its  original  value,  thereby  maintaining  a y+  value  of  one.  A slight  change  to  the  topology  was  also  made 
to  the  nacelles,  using  two  grids  for  the  external  and  internal  surfaces  rather  than  three.  This  change  permits  the 
user  of  PEGSUS  version  41_45,  which  is  roughly  seven  times  faster  than  version  41_36,  which  was  required  for 
the  successful  interpolation  of  the  initial  set  of  overset  grids. 

The  version  of  OVERFLOW  utilized  is  overl  ,7t.  This  version  contains  upgraded  treatments  of  the  multigrid  and 
mesh  sequencing  boundary  conditions. 


545 


• exp.  upper  surf 

■ exp.  lower  surf 

Spalart-Almaris 

Baldwin-Barth 


Wing/Body  Pressure  Coefficient  at  Mach  0.95 

Results  using  the  Spalart-Almaris  and  Baldwin-Barth  turbulence  models  for  the  baseline  wing/body 
configuration  at  Mach  0.95,  a = 3.5,  are  compared  against  the  experimental  results.  While  the  pressure 
distributions  are  virtually  identical  inboard,  by  53.9%  span  the  Spalart-Almaris  solution  shows  a slower  recovery 
from  the  acceleration  over  the  wing  leading  edge.  This  result  is  consistent  with  OVERFLOW  results  for 
transonic  transport  configurations,  where  the  transonic  shock  is  typically  further  aft  of  the  data  for  the 
Spallart-Almaris  model  than  for  Baldwin-Barth.  Forces  from  the  two  solutions  are  compared  against  data  in  the 
following  table: 


Cl 

Cd 

Cm 

experiment 

.12712 

.01184 

.00370 

Splart-Almaris 

.13106 

.01256 

.01007 

Baldwin-Barth 

.13161 

.01202 

.00071 

In  order  to  investigate  reports  suggesting  that  a y+  of  roughly  0.5  is  required  for  the  Spalart-Almaris  turbulence 
model  in  some  cases,  the  wall-normal  spacing  was  halved  and  the  grid  refined  by  10%.  No  significant 
difference  was  seen  in  either  solution. 

Since  discrepancies  with  the  data  are  larger  for  Spalart-Almaris  than  for  Baldwin-Barth,  all  subsequent 
transonic  solutions  utilize  the  Baldwin-Barth  turbulence  model. 
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• exp.  upper  surf 
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Wing/Body  Pressure  Coefficient  at  Mach  2.40 

Results  using  the  Spalart-Almaris  and  Baldwin-Barth  turbulence  models  for  the  baseline  wing/body 
configuration  at  Mach  2.40,  a = 3.0,  are  compared  against  the  experimental  results.  In  this  case  the  pressure 
distributions  are  virtually  identical  accept  at  the  leading  edge.  Once  again,  Spalart-Almaris  recovers  somewhat 
more  slowly  than  Baldwin-Barth,  but  in  this  case  the  Spalart-Almaris  solution  shows  better  agreement  with  the 
data.  Forces  from  the  two  solutions  are  compared  against  data  in  the  following  table: 


Cl 

Cd 

Cm 

experiment 

.07452 

.01080 

.00457 

Splart-Almaris 

.07414 

.01047 

.00421 

Baldwin-Barth 

.07394 

.01037 

.00420 

Since  discrepancies  with  the  data  are  larger  for  Baldwin-Barth  than  Spalart-Almaris,  all  subsequent  supersonic 
solutions  utilize  the  Spalart-Almaris  turbulence  model. 

While  the  solutions  generally  agree  with  experiment  quite  well,  a noticable  discrepency  occurs  on  the  lower 
wing  surface  at  32.6%  and  41 .3  % span. 
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Forces  and  Pitching  Moment 

tCA  Wing/Body 
M = 0.95  Re  = 4.0x10® 


• 1 6ft  TT:  runs  52, 63, 73 

- Baldwin-Barth 
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Wing/Body  Forces  and  Moments  at  Mach  0.95 

Solutions  with  the  Baldwin-Barth  turbulence  model  were  generated  for  the  baseline  wing/body  configuration  at 
Mach  0.95  , a = 2.96, 3.5,  and  3.92.  The  results  are  compared  against  experiment  in  the  figures.  The  lift  vs. 
angle  of  attack  curve  shows  the  slight  difference  in  slope  which  is  typical  of  aeroelastic  deformation  effects.  The 
lift  vs.  drag  polars  also  have  sligtly  different  trends  , with  the  discrepency  in  drag  being  around  4 counts  at  the 
lower  Cl’s  and  from  1 to  2 counts  at  higher  Cl’s. 
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Forces  and  Pitching  Moment 

TCA  Wing/Body 
M = 2.40  Re  = 4.0x1 06 


UPWT 

Spalart-Almaris 


Wing/Body  Forces  and  Moments  at  Mach  2.40 

Solutions  with  the  Spalart-Almaris  turbulence  model  were  generated  for  the  baseline  wing/body  configuration  at 
Mach  2.40 , a = 3.0, 3.26,  and  3.47.  The  results  are  compared  against  experiment  in  the  figures.  Once  again 
the  lift  vs.  angle  of  attack  curve  shows  the  slight  difference  in  slope  which  is  typical  of  aeroelastic  deformation 
effects.  The  lift  vs.  drag  polars  show  the  same  trend,  but  the  discrepency  in  drag  is  around  4 counts. 
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The  Effect  of  Trailing  Edge  Spacing 


The  initial  grid  has  a trailing  edge  spacing  on  the  order  of  0.5%  of 

chord.  The  sufficiency  of  this  was  checked  by  decreasing  the 

spacing  to  0.05%  of  chord  and  refining  the  streamwise  resolution 

by  1 0 /o.  Results  for  transonic  and  supersonic  flow  are  tablulated 
below. 


Mach 

T.E.  Spacing 

a 

Cl 

Cd 

0.95 

0.95 

0.5% 

0.05% 

3.5 

3.5 

.13205 

.13106 

.01204 

.01256 

2.40 

2.40 

0.5% 

0.05% 

3.0 

3.0 

.07414 
.0741 5 

.01047 

.01051 

Effect  of  Trailing  Edge  Spacing 

In  transonic  transport  problems,  solutions  from  OVERFLOW  have  been  found  to  be  sensitive  to  trailing  edge 
spacing,  particularly  in  terms  of  shock  location.  To  determine  the  sensitivity  for  HSR  configurations,  the 
streamwise  spacing  at  the  trailing  edge  was  decreased  from  0.5%  of  chord  to  0.05%  of  chord  and  the  grid  was 
refined  with  10%  more  grid  points.  The  results  at  both  supersonic  and  transonic  conditions  indicate  that  drag 
increases  by  about  half  a count  with  the  grid  refinement.  However,  no  significant  differences  in  pressure 
distributions  are  seen. 
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TCA 

Fuselage  & Wing  Cp’s 


M = 0.95 
Re  = 4.0x1 06 
a = 3.50° 
exp.  a = 3.66° 

WBNDF 


• exp.  upper  surf 

■ exp.  lower  surf 

Baldwin-Barth 


x/c  x/c 


Wing/Body/Nacelle/Diverter  Pressure  Coefficient  at  Mach  0.95 

Results  using  the  Baldwin-Barth  turbulence  model  for  the  baseline  wing/body/nacelle/diverter  configuration, 
including  the  fairing  where  the  nacelle  protrudes  through  the  upper  surface  of  the  wing,  at  Mach  0.95,  a = 3.5, 
are  compared  against  the  experimental  results.  While  the  computed  and  experimental  distributions  are 
comparable  inboard  and  immediately  outboard  of  the  diverters,  a large  discrepency  is  seen  in  the  location  of  the 
compression  between  the  diverters,  at  both  32.6%  and  41 .3%  span.  Further  outboard,at  68%  span,  the  effects 
of  aerolastic  deformation  are  seen  with  the  unloading  of  the  leading  edge. 
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• exp.  upper  surf 

■ exp.  lower  surf 

Spalart-Almaris 


Wing/Body/Nacelle/Diverter  Pressure  Coefficient  at  Mach  2.40 

Results  using  the  Spalart-Almaris  turbulence  model  for  the  baseline  wing/body/nacelle/diverter  configuration, 
including  the  fairing  where  the  nacelle  protrudes  through  the  upper  surface  of  the  wing,  at  Mach  2.40,  a = 3.0, 
are  compared  against  the  experimental  results.  Once  again,  the  computed  and  experimental  distributions  are 
comparable  inboard  and  outboard  of  the  diverters,  but  a large  discrepancy  is  seen  in  the  shock  location  at 
32.6%  span. 
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Forces  and  Pitching  Moment 

TCAWBND  * 16ft  TT:  runs  90,100,107 

M = 0.95  Re  = 4.0x1 06  -■ Baldwin-Barth 


Wing/Body/Nacelle/Diverter  Forces  and  Moments  at  Mach  0.95 

Solutions  with  the  Baldwin-Barth  turbulence  model  were  generated  for  the  baseline  wing/body/nacelle/diverter 
configuration  at  Mach  0.95  , a = 3.2, 3.5,  and  4.2.  The  results  exhibit  the  same  trend  the  wing/body  solutions, 
with  the  discrepancy  in  drag  varying  from  2 to  4 counts. 
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Forces  and  Pitching  Moment 

TCAWBNDF 
M = 2.40  Re  = 4.0x10® 


UPWT 

Spalart-Almaris 


Wing/Body/Nacelle/Diverter  Forces  and  Moments  at  Mach  2.40 

Solutions  with  the  Spalart-Almaris  turbulence  model  were  generated  for  the  baseline  wing/body  configuration  at 
Mach  2.40  , a = 3.0, 3.26,  and  3.5.  The  results  exhibit  the  same  trend  as  for  the  wing/body,  with  the  drag 
descrepancy  being  on  the  order  of  8 counts. 
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Natural  Flow  Wing  Design 


Based  on  Wing/Body  Euler  Analyses 

- Maintain  leading  edge  thickness  constraint 

- Blunt  leading  edge  outboard  of  leading  edge  break 

- Design  upper  surface  for  strong  conical  flow 

- Use  lower  surface  to  satisfy  spar  constraints 

Regridding  for  NFW  Designs 

- Given  fuselage  and  wing  definitions 

- Wing/Body  Script 

- Automatic  installation  script  to  place  nacelle/diverter  grids 
from  baseline  configuration  on  redesigned  wing. 


Natural  Flow  Wing  Design 

Designs  are  based  on  the  Natural  Flow  Wing  Design  philosophy  developed  by  Rick  Woods  and  Steve  Bauer,  as 
reported  in  the  1996  HSR  workshop  proceedings  under  the  title  "Application  of  the  Natural  Flow  Wing  Design 
Philosophy  to  the  HSR  Arrow  Wing  Configuration".  An  additional  aspect  of  the  current  design  is  to  maintain  a 
blunt  leading  edge  outboard  of  the  leading  edge  break  for  transonic  performance. 

Wing  body  grids  for  the  redesigned  configruation  are  generated  using  a modified  form  of  the  wing/body  script 
developed  in  the  AST  program.  Nacelle/diverter  grids  for  the  installed  configuration  are  generated  through  use 
of  a script  file  to  translate,  rotate  and  project  the  diverter/nacelle  component  grids  from  the  baseline 
configuration  onto  the  new  fuselage/wing  surface  definition.  Constraints  imposed  on  the  installation  are  to  keep 
the  same  inboard  diverter/wing  trailing  edge  intersection  point,  intersect  the  same  point  on  the  outboard  diverter 
with  the  wing  trailing  edge,  and  maintain  a constant  distance  between  the  nacelle  lip  and  wing  surface. 
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Upper  Surface  Pressure  Coefficient  on  the  Baseline  and  Design  NFW1  at  Mach  2.40 

The  first  redesign  of  the  TCA  baseline  configuration  using  the  Natural  Flow  Wing  design  philosophy  is 
designated  as  NFW1 . Pressure  distributions  on  the  wing  upper  surface  for  the  two  configurations  at  Mach  2.4 
are  illustrated  in  the  figure.  In  the  NFW  design,  a strong  conical  expansion  is  enforced  over  the  leadiing  edge. 
In  addition  the  compression  seen  on  the  outboard  section  of  the  baseline  is  eliminated.  The  compressions  seen 
at  the  front  of  the  nacelle  fairings  are  due  in  part  to  the  fact  that  the  nacelles  protrude  through  the  upper  surface 
to  a greater  extent  than  the  baseline.  However,  it  also  indicates  an  inadequate  regridding  process  in  the 
automatic  installation  of  the  nacelles. 
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Pressure  Coefficient 
Wing  Lower  Surface 


Lower  Surface  Pressure  Coefficient  on  the  Baseline  and  Design  NFW1  at  Mach  2.40 

Pressure  distributions  on  the  wing  lower  surface  for  the  two  configurations  at  Mach  2.4  are  illustrated  in  the 
figure.  The  NFW  design  exhibits  a much  stronger  shock  interaction  between  the  diverters,  as  well  as  a larger 
expansion  over  the  outboard  leading  edge  due  to  the  bluntness. 
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Pressure  Coefficient 
Wing  Upper  Surface 


Lower  Surface  Pressure  Coefficient  on  the  Baseline  and  Design  NFW1  at  Mach  0.95 

Pressure  distributions  on  the  wing  lower  surface  for  the  two  configurations  at  Mach  0.95  are  illustrated  in  the 
figure.  The  NFW  design  exhibits  a much  stronger  shock  at  the  nacelle  and  diverter  leading  edges,  as  well  as  a 
stronger  expansion  between  the  diverters. 
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Summary 

NFW1  has  poor  performance 

- Drag  is  3 counts  higher  than  the  baseline  at  Mach  2.4 

- Drag  is  23  counts  higher  than  the  baseline  at  Mach  0.95 

The  wing/body  script  is  fast  and  robust 

The  nacelle  rerigging  script  is  fast  but  incomplete 

- 2 hour  turnaround  on  grid  generation  for  NFW2 

- a sophisticated  nacelle  fairing  treatment  is  required 

Further  investigation  of  the  Spalart-Almaris  turbulence  model  is 
required  at  transonic  conditions 


Summary 

The  first  redesign  of  the  TCA  Baseline  Configuration  using  the  natural  flow  wing  design  has  poor  performance, 
with  the  drag  3 counts  higher  than  the  baseline  at  supersonic  cruise  and  23  counts  higher  at  transonic  cruise.  A 

portion  of  the  drag  increment  can  be  attributed  to  the  poor  installation  of  the  fairing  in  the  automatic  installation 
of  the  nacelles  and  diverters. 

The  wing  body  and  nacelle  rerigging  scripts  are  fast  and  robust  Turnaround  time  for  generating  grids  for  the 
second  design  was  on  the  order  of  2 hours  (this  does  not  include  the  time  to  get  PEGS  US  through  the  que). 
However,  the  nacelle/wing  fairing  treatment  is  incomplete,  requiring  a more  sophisticated  treatment  to  ensure 
captruing  the  nacelle  outer  boundary. 


While  Baldwin-Barth  is  the  preferred  method  for  transonic  OVERFLOW  solutions  at  this  time,  further 
investigation  of  the  Spalart-Almaris  turbulence  model  is  warranted. 
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TLNS3D/CDISC  MULTIPOINT  DESIGN 
OF  THE  TCA  CONCEPT 


Richard  L.  Campbell  and  Michael  J.  Mann 


1997  HSR  Aerodynamic  Performance  Workshop 
NASA  Langley  Research  Center 
February  25-28, 1997 


This  paper  presents  the  work  done  to  date  by  the  authors  on  developing  an  efficient  approach  to 
multipoint  design  and  applying  it  to  the  design  of  the  HSR  TCA  configuration.  While  the  title 
indicates  that  this  exploratory  study  has  been  performed  using  the  TLNS3DMB  flow  solver  and 
the  CDISC  design  method,  the  CDISC  method  could  have  been  used  with  any  flow  solver,  and  the 
multipoint  design  approach  does  not  require  the  use  of  CDISC.  The  goal  of  the  study  was  to 
develop  a multipoint  design  method  that  could  achieve  a design  in  about  the  same  time  as  10  anal- 
ysis runs. 
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OUTLINE  OF  PRESENTATION 


• Review  of  CDISC  design  method 

• Single-point  designs 

• Multipoint  design  approaches 

- weighted  average  of  geometries  (WAG) 

- transonic  flap 

• Concluding  remarks 


This  paper  will  begin  with  a review  of  the  Constrained  Direct  Iterative  Surface  Curvature 
(CDISC)  design  method,  then  look  at  its  application  to  design  of  the  TCA  configuration  at  a 
supersonic  and  a transonic  cruise  point.  Two  approaches  to  the  multipoint  design  problem  will 
then  be  considered:  a new  method  that  uses  a weighted  average  of  the  geometries  from  the  previ- 
ous point  designs  (referred  to  as  the  WAG  method),  and  a second  approach  that  involves  the  use 
of  a flap  on  the  supersonic  point  design  geometry  to  improve  the  transonic  performance.  The  con- 
cluding remarks  will  summarize  the  lessons  learned  so  far  and  present  the  future  plans  for  appli- 
cation of  the  multipoint  design  methods. 
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CDISC  DESIGN  METHOD 


• Efficient  constrained  design  with  Navier-Stokes 
codes  (design  run  time  * analysis  run  time) 

• Target  pressures  automatically  generated  based 
on  flow  constraints 

• Geometry  constraints  allow  design  to  be  impacted 
by  requirements  from  other  disciplines 

• “Optimization”  available  through  constraint  specifi- 
cation and  relaxation 

• Modular  coupling  of  CDISC  and  flow  solvers 


This  chart  gives  a description  of  some  of  the  key  features  of  the  CDISC  design  method.  It  is  a 
knowledge-based  design  approach  that  uses  rules  and  guidance  obtained  from  analytical,  experi- 
mental, and  computation  sources  to  allow  new  designs  to  be  obtained  in  about  the  same  time  as  a 
single  flow  analysis.  The  target  pressure  distributions  required  by  the  basic  DISC  method  are 
automatically  generated  based  on  flow  constraints,  and  a suite  of  geometry  constraints  are  avail- 
able that  allow  requirements  from  other  disciplines  such  as  structures  and  manufacturing  to  be 
included  in  the  design  process.  In  addition  to  the  TLNS3DMB  Navier-Stokes  code  already  men- 
tioned, the  CDISC  design  module  has  also  been  coupled  with  the  CFL3D,  OVERFLOW,  and 
USM3D  flow  solvers  to  allow  a variety  of  options  for  viscous  design  of  complex  configurations. 
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FLOWCHART  OF  CDISC  DESIGN  SYSTEM 


FLOW  SOLVER 


1 

■■H 

1 

8 

HUH 

This  flow  chart  shows  the  components  of  the  CDISC  design  module  and  how  it  is  coupled  with  a 
flow  solver.  The  design  process  begins  by  obtaining  a partially-converged  flow  solution  for  the 
initial  configuration.  The  CDISC  method  then  extracts  the  surface  geometry  and  pressure  infor- 
mation that  it  needs  from  the  grid  and  restart  files.  Initial  target  pressure  distributions  are  defined 
from  the  current  analysis  pressures,  then  modified  as  required  to  meet  the  flow  constraints.  After 
the  basic  DISC  method  is  used  to  alter  the  surface  geometry  based  on  these  target  pressures,  the 
geometry  constraints  are  applied  and  a grid  perturbation  scheme  is  used  to  modify  the  volume 
grid  to  accommodate  the  new  surface  shape.  This  new  grid  is  then  returned  to  the  flow  solver  for 
further  analysis. 
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SINGLE  POINT  DESIGN  USING  CDISC 

• Flow  solver:  TLNS3DMB  (Euler) 

• Grid:  117x25x25  C-H 

• Objective:  reduce  drag 

• Constraints: 

- twist  and  camber  changes  on  wing  only 

- maintain  original  lift  and  spanload  distribution 

• Design  variable:  chordwise  loading  parameter 


Since  the  primary  focus  of  this  study  was  the  development  and  evaluation  of  a multipoint  design 
method,  a simplified  approach  has  been  used  for  the  single  point  designs.  A coarse  wing/body 
grid  has  been  used  with  TLNS3DMB  run  in  an  Euler  mode  to  allow  rapid  flow  analysis  and  the 
design  was  limited  to  twist  and  camber  changes  on  the  wing  only.  The  objective  was  to  reduce  the 
drag  at  the  design  point  while  maintaining  the  original  lift  and  spanload  distribution.  For  these  ini- 
tial designs,  the  only  design  variable  was  the  chordwise  loading  parameter. 
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CDISC  DESIGN  STATIONS  ON  TCA6  PLANFORM 


Nine  wing  stations  are  used  for  design  in  all  cases,  with  the  locations  indicated  by  the  dashed  lines 
in  the  figure,  along  with  the  root  and  tip  stations.  The  root  station  is  fixed  in  order  to  maintain  the 
fuselage  geometry,  but  is  used  to  interpolate  changes  from  the  second  station  onto  grid  lines 
located  between  the  first  two  design  stations.  The  changes  at  the  tip  are  aliased  to  the  changes  at 
the  changes  at  the  station  just  inboard  of  it.  The  arrows  indicate  design  stations  at  which  sectional 
pressures  and  airfoils  will  be  shown  in  later  figures. 
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SUPERSONIC  DESIGN  USING  CDISC 


• Conditions: 

- M = 2.40 

- CL  = 0.0896 

• Results: 

- ACd  = -.00002 


The  supersonic  point  design  was  performed  at  a Mach  number  of  2.40  and  a lift  coefficient  of 
0.0896  (this  value  corresponds  to  lift  of  the  baseline  TCA  configuration  at  an  angle  of  attack  of 
3.5  degrees).  A number  of  combinations  of  chordwise  loading  parameter  were  tried,  but  little  drag 
reduction  was  obtained.  This  was  consistent  with  previous  studies  that  indicated  that  most  of  the 
potential  for  drag  improvement  involves  wing  root  design  and  thickness  changes. 
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SUPERSONIC  DESIGN  PRESSURES  AT  ti=0.43 
M=2.40  Cl=0.0896 


BASELINE 

DESIGN 


The  resulting  pressure  distributions  for  the  supersonic  design  are  compared  with  the  pressures  for 
the  baseline  configuration  at  a station  on  the  inboard  portion  of  the  wing.  The  leading  edge  expan- 
sion is  similar  for  both  configurations,  but  the  CDISC  design  pressure  distribution  is  smoother 
and  somewhat  more  aft-loaded.  It  should  be  noted  that  the  loading  could  be  shifted  slightly  more 
forward  without  a significant  drag  increase  if  pitching  moment  is  a constraint. 
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SUPERSONIC  DESIGN  AIRFOIL  AT  r|=0.43 
M=2.40  Cl=0.0896 


0 .1  .2  .3  .4  .5  .6  .7  .8  .9  1.0 


X/C 


A comparison  of  the  design  and  baseline  airfoils  at  the  same  inboard  station  shows  that  the  angle 
of  attack  is  nearly  the  same,  but  the  design  airfoil  has  more  camber,  which  may  be  helpful  at  the 
transonic  design  point. 
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SUPERSONIC  DESIGN  PRESSURES  AT  11=0.89 
M=2.40  Cl=0.0896 


x/c 


BASELINE 

DESIGN 


The  pressure  distributions  for  the  baseline  and  design  airfoils  at  a station  on  the  outboard  portion 
of  the  wing  are  very  similar.  Both  have  the  nearly  uniform  chordwise  loading  that  appears  to  be 
optimal  for  airfoils  on  a wing  with  a supersonic  leading  edge. 
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SUPERSONIC  DESIGN  AIRFOIL  AT  ti=0.89 
M=2.40  Cl=0.0896 


As  with  the  inboard  station,  the  baseline  and  design  airfoils  are  at  a similar  angle  of  attack,  with 
the  design  airfoil  having  slightly  more  camber. 
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TRANSONIC  DESIGN  USING  CDISC 

• Conditions: 

- M = 0.95 
-CL  = 0.1438 

• Results: 

- ACd  = -.00137 


The  transonic  design  was  performed  at  a Mach  number  of  0.95  and  a lift  coefficient  of  0.1438 
which,  as  with  the  supersonic  design,  was  the  lift  of  the  baseline  configuration  at  an  angle  of 
attack  of  3.5  degrees.  Again,  several  values  of  the  chordwise  loading  parameter  were  tried,  with 
the  best  case  producing  a drag  reduction  of  almost  14  counts  relative  to  the  baseline. 
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TRANSONIC  DESIGN  PRESSURES  AT  r|=0.43 
M=0.95  Cl=0.1438 


x/c 


BASELINE 

DESIGN 


The  baseline  pressure  distribution  at  the  inboard  station  has  a moderate  leading-edge  peak  and  an 
aft  re-expansion  that  creates  a fairly  steep  gradient  as  the  trailing  edge  is  approached.  The  design 
pressures  soften  both  of  these  features.  At  this  Mach  number,  the  sonic  pressure  coefficient  is 
about  -0.1.  Although  the  leading  edge  peak  as  well  as  the  aft  expansion  reach  supersonic  speeds, 
there  is  probably  not  much  wave  drag  associated  with  them  because  of  the  sweep  of  the  isobars 
(not  shown). 
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TRANSONIC  DESIGN  AIRFOIL  AT  ti=0.43 
M=0.95  Cl=0.1438 


x/c 


The  airfoils  that  correspond  to  the  pressures  in  the  previous  figure  are  shown  in  this  figure.  The 
angle  of  attack  is  nearly  the  same  for  the  two  configurations,  but  the  design  airfoil  has  quite  a bit 
more  camber,  especially  near  the  leading  edge.  This  suggests  that  a leading-edge  flap  could  be 
used  to  obtain  much  of  the  benefits  of  the  transonic  camber  design  at  this  station. 


574 


TRANSONIC  DESIGN  PRESSURES  AT  11=0.89 
M=0.95  Cl=0.1438 


BASELINE 

DESIGN 


At  the  outboard  design  station  near  the  tip,  the  baseline  pressures  indicate  that  a strong  leading- 
edge  shock  is  present.  The  design  has  eliminated  this  shock  while  recovering  the  lift  through 
increased  aft  loading,  although  the  aft  recovery  gradient  is  still  mild.  The  leading-edge  shock  was 
probably  the  source  of  most  of  the  wave  drag  for  the  baseline  configuration  at  these  conditions, 
and  its  elimination  is  the  likely  cause  of  the  drag  reduction  for  this  case. 
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TRANSONIC  DESIGN  AIRFOIL  AT  r|=0.89 
M=0.95  Cl=0.1438 


The  comparison  of  the  baseline  and  design  airfoils  at  the  outboard  station  indicates  that  the  reduc- 
tion in  the  leading-edge  pressure  peak  was  achieved  by  a combination  of  twist  and  camber,  with 
again  a noticeable  increase  in  camber  near  the  leading  edge.  This  suggests  that  only  some  of  the 
benefit  of  the  transonic  design  could  be  achieved  by  simply  deflecting  a leading-edge  flap. 
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MULTIPOINT  DESIGN  APPROACHES 


• Weighted  Average  of  Geometries  (WAG) 

- single  point  grids  used  as  design  variables 

- drag  estimation  for  improved  convergence 

with  combined-drag  objective  function 

- drag  constraint  option  also  available 

- procedure  is  fully  automated 

• Transonic  flaps 

- I.e.  flap  deflected  for  Mpeak  constraint 

- t.e.  flap  deflected  for  c(  constraint  (not  used) 


Having  completed  the  initial  single-point  designs  at  supersonic  and  transonic  cruise  conditions, 
the  problem  of  obtaining  good  performance  for  a mission  that  includes  flight  at  both  design  points 
was  then  addressed.  Two  approaches  to  this  multipoint  design  problem  were  considered:  a new 
method  that  employed  a weighted  average  of  the  point  design  geometries  (WAG  method),  and  a 
second  approach  that  used  a leading-edge  flap  to  improve  the  transonic  performance. 

The  WAG  method  automatically  determines  a grid  weighting  factor  that  will  be  used  to  blend  the 
grids  from  the  two  point  designs  in  a manner  that  will  minimize  a combined-drag  objective  func- 
tion (this  process  will  be  described  in  detail  in  the  next  figure).  A drag  estimation  procedure  is 
used  to  help  improve  the  convergence  of  this  process.  For  this  study,  the  supersonic  and  transonic 
drag  coefficients  are  arbitrarily  assigned  factors  of  0.8  and  0.2,  respectively,  for  computing  a com- 
bined-drag coefficient.  The  WAG  method  also  has  the  option  of  minimizing  the  drag  at  one  point 
subject  to  a drag  constraint  at  the  other  point. 

For  the  transonic  flap  approach,  a leading-edge  flap  is  simulated  on  the  supersonic  point  design 
geometry  to  reduce  the  drag  at  transonic  speeds.  In  this  case,  a maximum  Mach  number  constraint 
is  imposed,  and  the  flap  is  automatically  deflected  in  an  attempt  to  meet  the  constraint. 
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FLOW  CHART  FOR  WAG  MULTIPOINT  DESIGN  METHOD 


A flow  chart  of  the  WAG  multipoint  design  method  is  shown  in  this  figure.  The  procedure  begins 
with  a program  that  estimates  a grid  weighting  factor  for  use  in  the  grid  blending  program.  For  the 
first  two  passes  through  the  procedure,  weighting  factors  of  1 .00  and  0.00  are  specified  in  order  to 
obtain  flow  analyses  at  the  two  design  points  for  each  of  the  initial  point  designs.  For  the  third 
pass,  the  drag  is  assumed  to  vary  quadratically  between  configurations  at  a given  design  point 
(e.g,  at  the  supersonic  cruise  conditions,  a 50/50  blend  of  the  supersonic  and  transonic  point 
design  grids  will  yield  a drag  that  is  higher  than  the  drag  of  the  supersonic  design  by  25  per  cent 
of  the  drag  difference  between  the  transonic  and  supersonic  designs).  On  subsequent  passes,  a 
weighting  factor  for  the  drag  minimum  is  estimated  using  a curve  fit  of  the  previous  drag  values. 
Currently,  convergence  is  assumed  when  the  new  weighting  factor  is  within  0.05  of  a previous 
value. 
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MULTIPOINT  DESIGN  RESULTS  FOR  A GENERIC 
FIGHTER  AT  MILD  MANEUVER  CONDITIONS 


CDISC-WAG 

OPTIMIZATION 

BASELINE 


(DRA) 


Although  the  WAG  approach  is  a very  simple  one,  it  was  shown  in  one  study  to  produce  drag 
reductions  comparable  to  those  achieved  from  multipoint  numerical  optimization,  but  at  several 
orders  of  magnitude  less  cost.  In  this  study  done  cooperatively  with  the  Defense  Research  Agency 
(DRA),  a generic  fighter  wing/body  configuration  was  redesigned  to  minimize  the  combined-drag 
coefficient  for  a number  of  mission  (i.e.,  relative  drag)  weighting  factors.  The  two  design  points 
were  transonic  and  supersonic  mild  maneuver  conditions.  Since  different  grids  and  flow  solvers 
were  used  for  the  WAG  design  and  the  numerical  optimization,  the  design  results  have  been  nor- 
malized by  the  drag  of  the  baseline  configuration  at  each  design  point  as  predicted  by  the  flow 
solver  used  by  each  design  method  (the  circle  symbol  on  the  plot).  The  lines  represent  a series  of 
multipoint  optimizations  performed  by  each  method,  where  the  left  end  is  a purely  supersonic 
mission  and  the  right  end  is  a purely  transonic  mission.  Both  methods  were  effective  at  reducing 
drag  relative  to  the  baseline,  with  the  CDISC-WAG  method  showing  slightly  greater  drag  reduc- 
tions for  the  missions  dominated  by  the  supersonic  drag  point.  Each  multipoint  numerical  optimi- 
zation required  several  hundred  flow  analyses,  while  the  WAG  method  required  6-8  for  the  first 
mission  and  typically  only  two  additional  analyses  for  each  new  mission  after  that. 
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INITIAL  MULTIPOINT  DESIGN  RESULTS 
USING  WAG  METHOD 


C0 
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This  figure  shows  the  initial  results  from  applying  the  WAG  design  method  to  the  TCA  configura- 
tion. A grid  weighting  factor  (GWF)  of  0.0  corresponds  to  the  initial  transonic  point  design  while 
GWF=1.0  corresponds  to  the  initial  supersonic  point  design.  The  first  two  cycles  generated  the 
four  data  points  at  the  ends  of  the  dashed  lines,  with  the  dashed  lines  generated  by  the  quadratic 
drag  increase  assumption  mentioned  earlier.  Using  these  assumed  drag  distributions,  the  variation 
of  the  combined-drag  objective  function  is  computed  (solid  line)  and  a minimum  is  found  to  occur 
at  GWF=0.72.  The  symbols  at  this  value  of  GWF  indicate  the  computed  supersonic  and  transonic 
drag  values  for  this  new  blended  grid.  The  estimated  supersonic  drag  is  nearly  identical  to  the 
computed  value,  while  the  estimated  transonic  drag  is  slightly  below  the  computed  value. 
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FINAL  MULTIPOINT  DESIGN  RESULTS 
USING  WAG  METHOD 


Using  curve  fits  through  the  computed  values  of  drag,  a new  estimate  of  the  minimum  of  the 
objective  function  is  obtained  (GWF=0.80).  A new  blended  grid  is  generated  using  this  value  of 
GWF  and  analyzed  at  the  two  design  points.  The  resulting  drag  values  are  nearly  identical  to  the 
estimated  values.  Updating  the  drag  curve  fits  with  the  new  values  changes  the  minimum  of 
objective  function  by  less  than  0. 1 count  of  drag.  Thus  the  WAG  method  required  eight  analysis 
runs  to  achieved  a converged  multipoint  design. 
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TRANSONIC  FLAP  DESIGN 


• Conditions:  M = 0.95,  CL  = 0.1438 

• Geometry:  t\  = 0.66-1 .00,  hinge  line  x/c  = 0.15 

• Constraints:  Mpeak  < 1.3,  Mh(  < Mpeak 

• Results:  ACD  = -.00077 


As  an  alternative  to  a single  geometry  that  compromises  the  performance  at  the  two  design  points, 
the  design  at  the  dominant  design  point  (supersonic  in  this  case)  is  maintained  and  a flap  is  used  to 
improve  the  performance  at  the  other  design  point.  Since  most  of  the  wave  drag  at  the  transonic 
design  point  appeared  to  result  from  the  leading  edge  shock  on  the  outboard  portion  of  the  wing,  a 
flap  was  defined  from  the  leading-edge  planform  break  to  the  tip,  with  a hinge  line  at  x/c=0.15. 
(This  hinge  line  location  was  selected  based  on  a brief  study  of  several  locations).  The  CDISC 
flap  constraint  was  used  with  a peak  Mach  number  limit  of  1.3  which,  with  the  shock  sweep  con- 
sidered, gives  a normal  Mach  number  of  slightly  over  1 .0.  Also  included  in  the  flap  constraint  is  a 
requirement  that  the  hinge-line  Mach  number  not  exceed  the  leading-edge  peak  Mach  number. 
The  flap  deflection  quickly  converged  to  a value  of  5.84  degrees,  with  a resulting  drag  reduction 
of  about  8 counts  relative  to  the  baseline. 
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TRANSONIC  DESIGN  PRESSURES  AT  t|=0.89 
M=0.95  Cl=0.1438 


The  resulting  pressure  distribution  is  compared  with  the  pressures  from  the  transonic  single  point 
design  (solid  line)  and  the  WAG  multipoint  design  (dashed  line)  at  the  outboard  wing  station 
shown  before.  The  peak  Mach  number  limit  corresponds  to  a pressure  coefficient  of  about  -0.56. 
The  flap  design  did  not  quite  obtain  this  value  because  of  the  limitation  also  imposed  on  the 
hinge-line  Mach  number.  While  the  flap  design  did  not  eliminate  the  shock  as  was  the  case  with 
the  single  point  design,  the  shock  strength  was  greatly  reduced  relative  to  the  WAG  multipoint 
case.  (It  should  be  noted  that  the  WAG  case  did,  however,  reduce  the  shock  strength  relative  to  the 
baseline). 
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TRANSONIC  DESIGN  AIRFOILS  AT  ri=0.89 
M=0.95  Cl=0.1438 


x/c 


The  resulting  flap  geometry  at  the  outboard  station  is  compared  with  the  airfoils  from  the  single 
point  and  multipoint  designs  in  this  figure.  While  there  is  an  obvious  twist  difference  between  the 
single  point  design  and  the  flap  design,  the  camber  lines  in  the  region  of  the  flap  are  very  similar. 
The  large  surface  curvature  that  can  be  seen  at  the  hinge  line  is  the  source  of  the  secondary  Mach 
number  peak  shown  in  the  previous  figure. 
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SUMMARY  OF  MULTIPOINT  DESIGN 
DRAG  COEFFICIENTS 


Baseline 

WAG 

Flap 

Supersonic 

0.00624 

0.00625 

0.00622 

Transonic 

0.00672 

0.00637 

0.00595 

Combined 

0.00634 

0.00627 

0.00617 

A summary  of  the  drag  results  for  the  two  multipoint  design  approaches  is  given  in  the  table 
above.  At  the  supersonic  design  point,  the  WAG  design  is  only  slightly  worse  than  the  baseline, 
while  the  flap  design  retains  the  small  improvement  of  the  single  point  supersonic  design.  At  the 
transonic  design  point,  both  designs  have  improved  performance  relative  to  the  baseline,  with  the 
flap  approach  providing  about  twice  as  much  drag  reduction  as  the  fixed  geometry  WAG  design. 
The  variation  in  the  combined  drag  coefficient  is  small,  reflecting  the  dominance  of  the  supersonic 
portion  of  the  mission,  but  the  WAG  and  flap  multipoint  design  approaches  do  provide  about  1 
and  2 counts,  respectively,  of  overall  drag  improvement.  On  this  basis,  the  flap  approach  would  be 
preferred,  but  other  systems  issues  such  as  weight  or  safety  may  affect  this  conclusion. 
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CONCLUDING  REMARKS 


• An  efficient  automated  multipoint  design 
approach  (WAG)  has  been  developed 

• CDISC  camber/twist  design  at  supersonic  design 
point  produced  no  significant  drag  benefit 

• CDISC  camber/twist  design  was  effective  at 
reducing  drag  at  the  transonic  design  point 

• Flap  design  was  less  effective  than  camber/twist 
design  at  reducing  drag  at  the  transonic  design 
point,  but  produced  a slightly  lower  combined 
drag  value  than  the  WAG  approach 


In  conclusion,  an  efficient,  fully-automated  multipoint  design  approach,  referred  to  as  the  WAG 
method,  has  been  developed.  It  combines  single  point  design  geometries  (developed  in  this  study 
using  the  CDISC  design  method)  in  a systematic  fashion  to  reduce  a combined-drag  objective 
function.  For  the  limited  camber/twist  designs  performed  in  this  initial  study,  CDISC  did  not  pro- 
duce an  appreciable  drag  reduction  at  the  supersonic  design  point,  but  did  significantly  reduce  the 
drag  at  the  transonic  design  point.  A transonic  flap  approach  to  multipoint  design  was  also  evalu- 
ated in  this  study.  While  the  flap  was  not  as  effective  as  the  CDISC  point  design  at  reducing  tran- 
sonic drag,  it  did  produce  a lower  combined  drag  value  than  the  WAG  method. 
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FUTURE  PLANS 


• Include  thickness  and  fuselage/wing-root 
design  changes  (within  TCA  constraints)  in 
CDISC  point  designs 

• Apply  WAG  and/or  transonic  flap  multipoint 
design  approaches  to  promising  designs  on 
finer  viscous  grids,  including  full-configuration 
cases 


Based  on  the  results  of  this  initial  study,  it  is  recommended  that  future  design  work  with  CDISC 
allow  for  changes  in  wing  thickness  (within  the  TCA  constraints)  and  also  include  wing  root  and 
fuselage  design  changes,  especially  for  supersonic  point  designs.  While  this  relaxed  set  of  design 
constraints  should  initially  be  evaluated  using  the  coarser  grids  from  this  study,  the  efficiency  of 
the  CDISC  and  WAG  methods  should  allow  the  methods  to  be  practically  used  for  viscous  design 
on  fine  grids  for  even  full-configuration  cases. 
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Prediction  and  Assessment  of  Reynolds 
Number  Sensitivities  Associated  with  Wing 
Leading-Edge  Radius  Variations 

Richard  A.  Wahls 
Melissa  B.  Rivers 
Lewis  R.  Owens,  Jr. 

NASA  Langley  Research  Center 
February  25,  1997 


The  primary  objectives  of  this  study  were  to  expand  the  data  base 
showing  the  effects  of  LE  radius  distribution  and  corresponding 
sensitivity  to  Rn  at  subsonic  and  transonic  conditions,  and  to  assess  the 
predictive  capability  of  CFD  for  these  effects.  Several  key  elements  led 
to  the  initiation  of  this  project:  1 ) the  necessity  of  meeting  multipoint 
design  requirements  to  enable  a viable  HSCT,  2)  the  demonstration 
that  blunt  supersonic  leading-edges  can  be  associated  with 
performance  gain  at  supersonic  speeds  , and  3)  limited  data.  A test  of  a 
modified  Reference  H model  with  the  TCA  planform  and  2 LE  radius 
distributions  was  performed  in  the  NTF,  in  addtion  to  Navier-Stokes 
analysis  for  an  additional  3 LE  radius  distributions.  Results  indicate  that 
there  is  a tremendous  potential  to  improve  high-lift  performance  through 
the  use  of  a blunt  LE  across  the  span  given  an  integrated,  fully 
optimized  design,  and  that  low  Rn  data  alone  is  probably  not  sufficient 
to  demonstrate  the  benefit. 
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1997  HSR  Aerodynamic  Performance  Workshop 

Outline 

• Objectives 

• Background 

• Approach 

• Model  Geometry 

• NTF  Results 

• Navier-Stokes  Results 

• Concluding  Remarks 


As  outlined  above,  this  presentation  will  begin  with  a statement  of  the 
general  objectives  of  the  project,  followed  by  background  information 
which  led  to  the  initiation  of  the  study,  and  the  approach  taken  to  meet 
the  objectives.  Next,  the  wind  tunnel  model  is  described  including  its 
relationship  to  both  the  Reference  H and  Technology  Concept  Airplane 
(TCA)  geometries.  Finally,  preliminary  analysis  of  results  from  both  the 
experimental  and  computational  portions  of  the  study  will  be  discussed. 
Concluding  remarks  will  close  the  presentation. 
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Objectives 

• Expand  subsonic/transonic  data  base  of  Rn  sensitivities 
associated  with  LE  radius  variations,  including  the  supersonic 
LE  of  an  outboard  wing  panel 

• Obtain  data  for  CFD  code  validation 


The  general  objectives  of  the  project  are  shown  above.  The  primary 
goal  was  the  expansion  of  the  data  base  showing  the  effects  of  LE 
radius  distribution  and  corresponding  sensitivity  to  Rn  at  subsonic  and 
transonic  conditions.  Particular  emphasis  was  placed  on  the  under 
exploited  supersonic  LE  of  the  outboard  wing  panel.  The  experimental 
data  generated  meets  the  goal  of  data  for  CFD  validation. 

Additional  objectives  addressed  in  the  course  of  this  study,  but  not 
presented  herein,  included  preliminary  assessments  of  the  Rn  effects 
associated  with  the  planform  change  from  the  Reference  H to  the  TCA 
and  of  the  corresponding  change  to  the  high-lift,  inboard  LE  flap 
configuration.  These  topics  were  addressed  in  the  experimental  portion 
of  the  study,  and  results  are  described  in  a separate  paper  in  this 
workshop  (High-Lift  Session)  entitled  “Testing  a 2.2%  HSR  Ref.  H 
Model  with  Modified  Wing  Planform  in  the  NTF,”  by  Owens,  Wahls,  and 
Hamner. 
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Background  I 

• Multipoint  Design  Requirements 

- Take-off,  Approach,  Transonic  Cruise,  & Supersonic  Cruise 

- High-lift  needs  quieter  engines  and/or  more  L/D  to  be  viable 

• Supersonic  Issues 

- wave  drag  increase  with  LE  bluntness  for  supersonic  LE? 

- reduced  drag  supersonic  airfoil  design  w/  blunt  LE  reported  by 
Wilby 

- NFW  success  on  M2.4-7A  configuration 

- improved  multipoint  design  incorporating  a blunt  supersonic  LE  is 
possible 

• Existing  Data 

- AST  210  configuration  tested  in  the  NTF 

• Mach  = 0.3r  LE  radius  variation  inboard  only:  “sharp”  vs  “blunf 

- Generic  65  deg  delta  wing  tested  in  the  NTF 

• Mach  = 0.4  -»  0.9,  4 interchangeable  LE’s  w/  various  radius 
distributions 


Several  key  elements  led  to  the  initiation  of  this  project:  1 ) the  necessity 
of  meeting  multipoint  design  requirements,  2)  the  demonstration  that 
blunt  supersonic  leading-edges  can  be  associated  with  performance 
gain  at  supersonic  speeds  , and  3)  limited  data. 

In  addition  to  supersonic  cruise,  the  mission  of  an  HSCT  also  includes 
transonic  cruise,  take-off  and  landing  requirements.  Currently,  high-lift 
success  requires  quieter  engines  and/or  more  L/D.  Blunt  LE  geometry 
can  enhance  performance  at  these  conditions,  particularly  if  the  less 
swept,  outboard  portion  of  the  wing  is  allowed  to  contribute. 

Can  improved  subsonic  performance  be  realized  without  adversely 
affecting  supersonic  performance?  Yes.  Wilby  (Aeronautical  Research 
Council,  CP-921 ) shows  a reduced  drag  airfoil  at  supersonic  speeds 
incorporating  a relatively  blunt  LE.  The  Natural  Flow  Wing  (NFW) 
design  philosophy  as  applied  in  the  redesign  of  the  M2.4-7A 
configuration  (NASA/CP-1999-209690,  Bauer  and  Krist)  demon- 
strates the  possibilities  for  an  HSCT.  The  key  to  success  is  an 
integrated  design;  that  is,  the  LE  geometry  cannot  be  changed 
independent  of  overall  geometry. 

Data  addressing  LE  radius  effects  over  a large  Rn  range  are  limited. 

Two  sources  are  NTF  tests  of  the  AST  210  (73/60  arrow  wing)  with 
sharp  and  blunt  inboard  LE  radi  (outboard  LE  remained  sharp) (NASA/ 
TP-1999-209695,  Williams,  et  al.),  and  of  a generic  65  deg  delta  wing  with 
4 interchangeable  LE  sections  (NASA  TM  4645,  Chu  and  Luckring). 
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Background  II 

• HSR  Program 

- Baseline  configuration  changed  from  Ref  H to  TCA 

- Timing  issues 

• affect  next  downselect  to  Tech.  Configuration 

• material  availability,  new  vs.  modified  model,  NTF  shutdown 

• want  to  demonstrate  that  a blunt  supersonic  LE  design  is  worth 
pursuing  when  considering  multipoint  design  requirements 

TCA  lints  frozen  NTF089  Data  to  T1  Tach.  Config.  Downselsrf 


NTF  Shutdown  NTF  Shutdown 
Starts  Ends 


The  HSR  program  is  currently  in  a 3 year  phase  centered  around  the 
evaluation  and  redesign  of  the  TCA  configuration.  It  was  desired  to 
generate  Rn  effects  data  on  the  TCA  planform,  examine  the  high-lift  LE 
flap  configuration,  and  demonstrate  that  a blunt  supersonic  LE  design  is 
worth  pursuing  in  time  to  provide  input  to  the  definition  of  the  follow-on 
baseline  configuration.  Given  the  NTF  schedule  and  major  shutdown 
for  upgrade,  model  material  availability,  and  insufficient  funds/support 
for  a new  model,  the  decision  was  made  to  target  a test  window  in  the 
NTF  in  the  1st  quarter  of  FY97  prior  to  the  NTF  shutdown. 
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Approach 

• Modify  2.2%  Ref  H model  to  TCA  planform 

- include  alternate  LE  radius  distribution  & high-lift  flaps 

• Execute  NTF  test 

- Rn  effects  assessment  & CFD  validation  data 

• Execute  transonic  Navier-Stokes  analysis 

- prediction  assessment  for  Rn  and  LE  radius  effects 


The  approach  to  meet  the  objectives  within  the  program  and  facility 
availability  constraints  was  as  follows.  First,  modify  an  existing  model 
suitable  for  the  NTF  test  environment.  The  obvious  choice  was  the 
2.2%  HSR  Reference  H model.  Second,  execute  a test  in  the  NTF  at 
high-lift  and  transonic  conditions  to  provide  a wide  range  of  Rn 
conditions  to  allow  experimentally  based  assessments  and  provide  data 
for  CFD  validation.  Finally,  and  concurrent  with  the  experimental 
testing,  execute  a complementary  CFD  study  to  assess  predictive 
capability  and  to  expand  the  study  to  geometries  not  tested 
experimentally. 
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Model  Geometry  I 

• Comparison  of  Modified  Ref  H and  Ref  H Models 

• Geometric  Constants  at  2.2%  scale 

S.ref  mac  span  AR  LE  sweep 

fr  (gross)  in.  in.  deg 

Ref.  H 3.674  22.71  34.23  2.21  76/68.5/48 

Modified  Ref.  H.  4.114  25.07  34.65  2.03  71/52 
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The  first  step  was  the  modification  of  the  existing  2.2%  HSR  Ref.  H 
model  to  represent  the  TCA  wing  as  closely  as  possible.  Geometric 
constants  are  shown  above;  the  Modified  Ref.  H values  are  identical  to 
the  TCA.  Note,  that  the  reference  area  for  the  Ref.  H is  the  gross  wing 
area  (rather  than  the  wimpress  area  used  during  Ref.  H testing)  to  be 
consistent  with  the  TCA  definition.  The  Ref.  H (truncated)  body  and 
inboard  wing  center  section  and  TE  (indicated  by  the  dotted  lines)  were 
maintained,  while  the  LE  and  outboard  wing  panels  (indicated  by  the 
dashed  lines)  were  not.  New  LE  and  outboard  wing  panels  were 
designed  and  fabricated  to  provide  the  TCA  planform  while  not 
restricting  a return  to  the  Ref.  H geometry. 
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Model  Geometry  II 


• Airfoil  modification  process  was  as  follows: 

- align  TCA  & Ref.  H TE  (inboard  sections:  existing  model) 

- rotate  TCA  section  around  TE  to  align  to  with  existing  Ref.  H parts 

- blend  overlap  section  between  TCA  LE  and  existing  Ref.  H parts 

- spanwise  blending  outboard  of  existing  Ref.  H parts 
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y a 5.7646  in.  (2.2%  SCALE) 
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The  modification  process,  or  more  specifically  the  blending  process,  is 
demonstrated  above  for  a typical  inboard  airfoil  section.  First,  the  TCA 
section  at  a given  span  location  is  translated  to  match  the  TE  of  the 
existing  Ref.  H model  hardware.  Next,  the  TCA  section  is  rotated 
around  the  TE  to  align  with  the  existing  model  parts  with  emphasis  on 
the  upper  surface  to  avoid  unwanted  surface  inflections.  Finally, 
blending  occurs  over  a small  region  forward  of  the  existing  hardware  in 
to  the  TCA  LE  region.  This  sequence  was  repeated  for  several  airfoils 
over  the  span  of  the  existing  wing  center  section/TE  hardware; 
outboard  of  this  point,  a small  blending  region  existed  in  the  spanwise 
direction  until  the  TCA  outboard  airfoil  definitions  could  be  maintained. 
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The  resulting  geometry  had  the  characteristics  shown  above.  Note  that 
wing  LE  radius  distribution  of  the  modified  Ref.  H is  identical  to  that  of 
the  TCA,  and  that  both  the  TCA  and  the  Ref.  H have  a sharp  LE  on  the 
outboard  wing  panel.  Existing  Ref.  H model  hardware  inboard  drives 
the  differences  in  wing  twist,  maximum  thickness,  and  the  location  of 
the  maximum  thickness.  Outboard  of  the  pre-existing  hardware,  the 
modified  Ref.  H and  TCA  geometries  more  closely  match. 

The  resulting  geometry  was  smooth  and  sufficient  to  address  the 
objectives  of  the  study.  However,  in  no  way  should  this  geometry  be 
considered  optimized  aerodynamically. 


596 


1997  HSR  Aerodynamic  Performance  Workshop 

Model  Geometry  IV 

• Comparison  of  Modified  Ref.  H Airfoil  Sections 

- Baseline  & “Full  Blunt”  LE’s  tested  in  the  NTF 

- shown  in  2.2%  model  dimensions 

• y = 13  inches  — »r|=  75% 


• BasaMn*  (TCA  LE) 

- Alternate  (FuB  Bkint  LE) 
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Once  the  baseline,  modified  Ref.  H geometry  was  established,  several 
alternative  LE  radius  distributions  were  quickly  assessed  using  the 
linear  theory  code,  AER02S,  by  Carlson,  et.  al.  Time  constraints 
permitted  the  selection  of  one  alternative  radius  distribution  for 
fabrication  and  testing.  The  chosen  alternative  is  referred  to  as  the  “full 
blunt”  LE,  which  is  characterized  by  an  inboard  LE  radius  identical  to 
the  baseline,  but  a constant  LE  r/c  outboard  of  the  crank  and  matching 
that  at  the  crank.  Shown  above  is  a comparison  of  airfoils  at  the  75% 
semispan  station  highlighting  the  increased  LE  bluntness  of  the 
alternative  distribution.  The  alternative  LE  blends  into  the  baseline 
airfoil  forward  of  the  maximum  thickness  location. 

Three  other  alternative  LE  radius  distributions  evaluated 
computationally  will  be  described  later. 
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NTF  Test  Variables 


Mach  = 0.30 

Rn,mac  = 9.4  -4  1 00  x 1 06 
a = -3°  ->  24° 
nacelles  on 

0/0  & 30/10  full-span  flaps 
baseline  & alt.  LE  radius 


• Mach  = 0.90 

• Rn,mac=  11  -*89  x 106 

• a = -2°  1 2° 

• nacelles  off 

• 0/0  flaps 

• baseline  & alt.  LE  radius 


The  range  of  test  conditions  in  the  NTF  test  (designated  NTF089) 
pertinent  to  this  study  are  shown  above.  All  data  shown  herein  were 
obtained  with  natural  transition  on  the  wing.  A complete  set  of  low  Rn 
data  with  fixed  transition  was  planned  but  not  obtained  due  to 
significant  facility  downtime  associated  with  a pitch  system  failure. 

Force  and  moment  data  were  obtained.  Limited  pressure  data  on  the 
existing  Ref.  H wing  center  section  were  also  obtained;  LE  and 
outboard  wing  panel  pressures  were  not  obtained  due  to  limited  funding 
and  design/fabrication  time  constraints. 
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NTF  Results  la 


• Comparison  of  alternate  & baseline  LE  effect  on  CD 

- Mach  = 0.9,  0/0  flaps 

- Generally  a small  but  favorable  effect  due  to  the  alternate  LE 
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The  above  chart,  and  the  following  chart,  show  data  as  a function  of  LE 
radius,  Rn,  and  CL  at  Mach  = 0.90  with  undeflected  flaps.  In  each  chart, 
the  left  hand  plot  shows  data  for  each  LE  as  a function  of  CL  at  low  Rn. 
The  right  hand  plot  shows  the  increment  due  to  the  alternate  LE  radius 
for  each  Rn  as  a function  of  CL.  Data  has  not  been  corrected  for 
aeroelastic  effects  (wing  twist  and  bending  differences  due  to  two 
significantly  different  dynamic  pressures  needed  to  span  the  Rn  range). 

The  incremental  drag  data  above  generally  show  a drag  reduction  due 
to  the  alternate  LE,  particularly  below  CL  = 0.3.  The  large,  adverse 
increment  in  the  22e6  and  55e6  data  at  high  CL  are  currently  attributed 
to  curve  fitting  uncertainty  due  to  sparse  data  in  this  range.  This 
explanation  will  be  further  investigated. 
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NTF  Results  lb 


• Comparison  of  alternate  & baseline  LE  effect  on  L7D 

- Mach  = 0.9,  0/0  flaps 

- Generally  a small  but  favorable  effect  due  to  the  alternate  LE 
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The  L/D  incremental  data  are  consistent  with  the  CD  data  of  the 
previous  chart,  indicating  generally  small,  favorable  effects  due  to  the 
alternate  LE.  Results  are  particularly  favorable  at  CL’s  near  UDmax. 


600 


1997  HSR  Aerodynamic  Performance  Workshop 


NTF  Results  lla 


• Comparison  of  alternate  & baseline  LE  effect  on  CD 


- Mach  = 0.3,  0/0  flaps 

- Generally  a significant  favorable  effect  due  to  the  alternate  LE 


below  CL  = 0.5 
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The  above  chart,  and  the  following  chart,  show  data  as  a function  of  LE 
radius,  Rn,  and  CL  at  Mach  = 0.30  with  undeflected  flaps.  In  each  chart, 
the  left  hand  plot  shows  data  for  each  LE  as  a function  of  CL  at  low  Rn. 
The  right  hand  plot  shows  the  increment  due  to  the  alternate  LE  radius 
for  each  Rn  as  a function  of  CL.  Data  has  not  been  corrected  for 
aeroelastic  effects  (wing  twist  and  bending  differences  due  to  two 
significantly  different  dynamic  pressures  needed  to  span  the  Rn  range). 

The  incremental  drag  data  above  generally  show  a significant  drag 
reduction  due  to  the  alternate  LE,  particularly  below  CL  = 0.5.  The 
large,  adverse  increment  in  the  55e6  data  at  high  CL  are  currently 
attributed  to  curve  fitting  uncertainty  due  to  sparse  data  in  this  range. 
This  explanation  will  be  further  investigated. 
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NTF  Results  lib 


• Comparison  of  alternate  & baseline  LE  effect  on  L/D 


- Mach  = 0.3,  0/0  flaps 

- Generally  significant  favorable  effect  due  to  the  alternate  LE,  but 
decreasing  with  increasing  CL  above  L/Dmax 
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The  L/D  incremental  data  are  consistent  with  the  CD  data  of  the 
previous  chart,  indicating  significant,  favorable  effects  due  to  the 
alternate  LE.  Results  are  particularly  favorable  at  CL’s  near  UDmax, 
but  decrease  with  increasing  CL. 
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NTF  Results  Ilia 

• Comparison  of  alternate  & baseline  LE  effect  on  CD 

- Mach  = 0.3,  30/10  flaps 

- Significant  drag  reduction  due  to  the  alternate  LE  in  the  vicinity  of 
CL=  0.5  (-design) 

Rn.mac  = 9.4e€  <*•!«»  « (alternate  - baseline)  LE 


The  above  chart,  and  the  following  chart,  show  data  as  a function  of  LE 
radius,  Rn,  and  CL  at  Mach  = 0.30  with  leading  and  trailing  edge  flaps 
deflected  30  and  10  degrees,  respectively.  In  each  chart,  the  left  hand 
plot  shows  data  for  each  LE  as  a function  of  CL  at  low  Rn.  The  right 
hand  plot  shows  the  increment  due  to  the  alternate  LE  radius  for  each 
Rn  as  a function  of  CL.  Data  has  not  been  corrected  for  aeroelastic 
effects  (wing  twist  and  bending  differences  due  to  two  significantly 
different  dynamic  pressures  needed  to  span  the  Rn  range). 

The  incremental  drag  data  above  generally  show  a significant  drag 
reduction  due  to  the  alternate  LE,  particularly  near  CL  = 0.5. 
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NTF  Results  lllb 

• Comparison  of  alternate  & baseline  LE  effect  on  L/D 

- Mach  = 0.3,  30/10  flaps 

- Generally  a favorable  effect  due  to  the  alternate  LE 


The  L7D  incremental  data  indicate  favorable  effects  due  to  the  alternate 
LE.  Results  are  particularly  favorable  in  the  vicinity  of  CL  = 0.5 
(-design). 
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NTF  Results  IV 

• Comparison  of  CD  & LVD  increments  as  a function  of  Rn 

- positive  benefit  due  to  blunt  supersonic  LE  across  the  Rn  range 

- low  Rn  data  alone  is  not  sufficient  to  demonstrate  benefit 


The  above  chart  shows  CD  and  L/D  incremental  data  due  to  the 
alternate  LE  at  Mach  = 0.3  and  0.9  and  CL  = 0.2  (near  L/Dmax)  with 
undeflected  flaps  cross  plotted  against  Rn.  Jumps  in  increment  values 
at  constant  Rn  (~33e6)  are  due  to  uncorrected  aeroelastic  effects. 
Generally,  the  alternate  LE  provides  a positive  benefit  across  the  Rn 
range  as  exhibited  by  reduced  CD  and  increased  UD.  The  benefits  are 
considerably  more  significant  at  the  low  Mach  number,  indicating 
tremendous  potential  to  enhance  high-lift  performance  with  an 
integrated  design.  Nonmonotonic  trends  with  Rn  are  not  understood  at 
this  point,  but  remain  under  investigation.  Another  key  element 
demonstrated  above  is  the  fact  that  low  Rn  data  alone  is  not  sufficient 
to  demonstrate  the  potential  benefit  of  the  alternate  LE  geometry. 
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CFD  Approach 


• Code  & Grid  Test  Matrix 

- CFL3D  v4.1  - Mod.  Ref.  H w/  5 LE  radius  distrib. 

- BUDS  turb.  model  - Mach  = 0.90 

- C-0  grid,  single  block  - Rn.mac  = 10,  30,  80e6 

- 141  x257x65  - a = 1, 3, 5,  7, 10  deg 


The  test  matrix  for  the  computational  study  included  5 geometries 
(baseline  Modified  Ref.  H and  4 alternate  LE  radius  distributions),  5 
angles-of-attack  (1 ,3,5,7,  and  10  deg),  3 Rn  (10e6,  30e6,  and  80e6 
based  on  the  mac),  and  one  Mach  number  (0.90).Navier-Stokes 
predictions  were  made  using  CFL3D,  version  4.1 , and  primarily  the 
Baldwin-Lomax  with  Degani-Schiff  turbulence  model.  The  Baldwin- 
Barth,  Spalart-Allmaras,  and  Mentei^s  SST  turbulence  models  were 
also  used  on  the  baseline  geometry  at  an  angle-of-attack  of  5 deg.  The 
single  block  C-0  grid  has  141  points  streamwise,  257  spanwise,  and  65 
normal  to  the  surface  for  total  approaching  2.4  million  grid  points. 
Normal  spacing  near  the  surface  was  held  constant  for  the  Rn  = 1 0e6 
and  30e6  cases,  and  modified  for  the  Rn=80e6  case;  in  all  cases,  the 
average  y-plus  value  is  no  larger  than  1 . C90  times  for  converged 
solutions  ranged  from  2 to  1 1 hours,  with  angle-of-attack  being  the 
dominant  factor,  followed  by  the  Rn,  and  finally  the  LE  geometry. 

Timing  for  the  code  was  approximately  5 microseconds/grid  point/ 
iteration. 
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The  LE  radius  detail  for  the  baseline  and  4 alternates  studied  are 
shown  here.  Recall  that  the  baseline  LE  is  identical  to  that  of  the  TCA. 
The  “full  blunt”  LE  matches  the  inboard  LE  of  the  baseline,  then 
maintains  the  blunt  r/c  ratio  at  crank  across  the  entire  outboard  panel. 
The  “full  sharp”  LE  is  identical  to  the  baseline  on  the  outboard  panel,  in 
addition  to  being  sharp  inboard  (r=0).  The  “blunt  to  sharp”  LE  begins 
with  baseline  radius  at  the  side  of  body,  then  “r”  linearly  decreases  to 
zero  at  the  tip.  Finally,  the  “sharp  to  blunt”  LE  begins  as  sharp  (r=0)  at 
the  side  of  body,  then  linearly  increases  to  that  of  the  “full  blunt”  LE  at 
the  tip.  Recall  that  the  baseline  and  “full  blunt”  geometries  were  tested 
in  the  NTF. 
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CFD  Results  - Validation 


• Comparison  of  predicted  and  experimental  drag  increments 

- Mach  = 0.9,  0/0  flaps,  no  aeroelastic  correction  to  expt.  data 

• Prediction  levels  seem  optimistic,  trends  up  to  Rn.mac  - 55e6 
are  similar 


The  chart  above  shows  a comparison  of  drag  increments  due  to  the  “full 
blunt”  LE  at  a = 5 deg  (CL-0.2)  between  experiment,  CFL3D,  and 
AER02S  (AER02S  was  used  for  a quick  parametric  study  to  choose 
the  primary  alternate  LE).  Experimental  data  has  not  been  corrected  for 
aeroelastic  effects;  the  vertical  data  shift  at  Rn.mac  ~ 33e6  represents 
the  difference  between  test  dynamic  pressures  of  approximately  1000 
and  1800  psf.  The  prediction  levels  seem  optimistic,  but  the  trends  with 
Rn  are  similar,  particularly  up  to  a Rn.mac  ~ 55e6.  Although  not  done 
as  yet,  correction  of  the  experimental  data  to  that  of  the  rigid  geometry 
shape  (q  = 0 p$f)  used  in  the  predictions  will  shift  the  experimental 
results  toward  the  predicted  results. 

An  additional  factor  not  accounted  for  in  this  comparison  is  the  fact  that 
the  “sharp”  LE  of  the  experimental  model  is  in  reality  a finite  thickness 
on  the  order  of  0.0075  inches  while  the  predictions  shown  are  for  a 
perfectly  sharp  LE.  A quick  order  of  magnitude  analysis  for  this 
difference  was  made  using  AER02S  with  a “sharp”  LE  radius  defined 
as  half  the  finite  thickness.  The  results  indicate  an  effect  on  the  order  of 
5-6  drag  counts  in  the  direction  of  reducing  the  disagreement  shown 
above. 
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CFD  Results  - Assessment  I 

• Effect  of  alternate  LE  radius  distributions 
- Mach  = 0.90,  0/0  flaps 

• LE  bluntness  & removal  of  discontinuity  in  radius  distribution 
favorable 

Q on.  delta  = (’ alternate ' - baseline)  LE 


The  chart  above  compares  the  effects  of  the  4 alternate  LE  geometries 
investigated.  The  missing  point  at  Rn,mac  = 80e6  for  the  “full  sharp”  LE 
is  due  to  the  lack  of  a converged  solution.  This  occurred  at  angles-of- 
attack  of  1 , 5,  and  1 0 deg;  angles-of-attack  of  3 and  7 deg  were  not 
attempted. 

The  two  LE’s  with  the  most  “sharpness”  show  degraded  performance, 
while  the  two  LE’s  with  the  most  “bluntness”  show  improved 
performance.  The  best  LE  is  the  “blunt  to  sharp”  which,  relative  to  the 
baseline,  has  slightly  reduced  bluntness  just  inboard  of  the  crank,  but 
significantly  more  outboard.  The  other  characteristic  of  this  LE  is  the 
smooth,  continuous  variation  of  “f  across  the  crank;  this  may  be  why  it 
out  performs  the  “full  blunt”  LE,  which  is  blunter  but  still  has  some 
geometric  discontinuity  at  the  crank. 
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1997  HSR  Aerodynamic  Performance  Workshop 

CFD  Results  - Assessment  lla 

• Upper  surface  pressure  difference  due  to  the  “full  blunt”  LE 
- Rn,mac  = 10,  30,  and  80  million,  Mach  = 0.9,  a = 5 deg 

• Relatively  lower  pressure  in  the  LE  region 

• Significant  differences  emanating  from  the  crank 


The  chart  above  shows  upper  surface  pressure  differences  due  to  the 
‘lull  blunt”  LE  on  the  outboard  wing  panel.  Although  somewhat  difficult 
to  see,  the  LE  region  has  increased  suction  levels  due  to  the  bluntness 
The  more  easily  seen  differences  emanate  from  the  crank.  The  “full 
blunt”  LE  has  a smoother  geometric  transition  across  the  crank  region 
than  does  the  baseline,  and  in  turn  affects  the  vortex  formation 
emanating  from  the  crank.  The  comparisons  also  show  significant 
differences  with  increasing  Rn  as  well. 
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19^7  HSR  Aerodynamic  Performance  Workshop 

X&SfmdSf—*  . . 1 " '''  

Concluding  Remarks 

• The  Modified  Ref.  H model  tested  is  in  no  way  aerodynamically 
optimized 

• Smooth,  continuous  transition  from  inboard  to  outboard  (at  the 
crank)  is  important 

• Computational  modelling  of  “sharp”  LE’s  vs.  physical  model 
geometry  may  be  important 

• Positive  benefit  due  to  blunt  supersonic  LE  across  the  Rn  range 

• Low  Rn  data  alone  is  not  sufficient  to  demonstrate  benefit 

• Multipoint  designs  incorporating  blunt  supersonic  LE’s  are 
definitely  worth  pursuing 

— supersonic  improvement  (NFW) 

- tremendous  potential  for  high-lift  conditions 


In  conclusion,  the  following  points  are  reiterated.  The  modified 
Reference  H model  with  the  TCA  planform  as  tested  should  in  no  way 
be  considered  an  aerodynamically  optimized  configuration.  Rather,  it  is 
a test  bed  for  examining  sensitivities  to  localized  geometric  changes. 
Results  indicate  the  benefits  of  a continuous  and  smooth  LE  radius 
distribution  across  the  entire  span.  Accurate  modelling  of  LE 
geometries,  in  particular  that  of  sharp  LE’s,  may  significantly  affect 
computational  predictions. 

Most  importantly,  the  current  study  demonstrated  positive  benefit  due  to 
the  use  of  a blunt  supersonic  LE  across  the  Rn  range,  and  that  low  Rn 
data  alone  is  not  sufficient  to  demonstrate  this  benefit.  Multipoint 
designs  incorporating  blunt  supersonic  LE’s  are  definitely  worth 
pursuing,  both  at  fundamental  and  configuration  levels  of  study.  The 
NFW  wing  design  philosophy  incorporates  this  idea  by  design  and  has 
been  successful  at  reducing  drag  even  at  supersonic  cruise  conditions. 
The  tremendous  potential  for  gain  at  high-lift  conditions  alone  warrants 
further  study  be  given  to  LE  radius  effects,  and  these  studies  should 
include  a wide  range  of  Rn. 
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Lateral  Characteristics 
Summary 
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1 .5%  TCA  Controls  Model  (#20) 
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.5%  TCA  Controls  Model  (#20) 


Extended  Aftbody 


Reference  H / TCA  Model  Comparison 
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Ref  H Vertical  Tail 
Planform  with  Rudder 


Model  Highlights 
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Trip  dot  height(T.s.#2):  0.012”  0.0085 

TCA  Pressures:  (2)  Inb.  and  (2)  Outb.  Nacelle  Base 

(2)  Model  Chamber 


Test  Highlights 
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Pitch  Polars:  -4°  < aF<  1 2°  for  p = 0°,  ±3° 

Yaw  Polars:  -6°  < p < 6°  for  aF  = 0°,  4°,  8 


TCA  Configurations  Tested 
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Longitudinal  Control  (W/B/N/D/H/V): 

- iH  = ±6°,  ±15°;  8 =±20° 
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Tested  at  M=1.8  and  2.1  only 
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Long-Term  Repeatability 
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From  these  data,  it  was  concluded  that  the  long-term  repeatability  is  goo 
Mach  number. 
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Longitudinal  Characteristics 
Component  Build-up 
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Longitudinal  Characteristics 
Component  Build-up 
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Longitudinal  Characteristics 
Component  Build-up 
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Configuration 


Longitudinal  Characteristics 
Component  Build-up 
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Lateral  Characteristics 
Stability  Comparison 
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Lateral  Characteristics 
Control  Effectiveness 
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in  either  direction  is  approximately  the  same  in  magnitude. 
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Lateral  Characteristics 
Control  Comparison 
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Summary 
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(4)  The  directional  stability  of  the  TCA  is  improved  relative  to  Reference  H at 
higher  angles-of  attack.  The  chine  is  effective  for  improving  directional  stability. 
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Increased  at  higher  aF’s  with  and  without  chines 


Summary  (Cont’d) 
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edge  flap  deflections 


HSR  - 1997  Aerodynamic  Performance  Workshop 


Effect  of  Boattail  and  Sidewall  Curvature 
on  Nozzle  Drag  Characteristics 

Francis  J.  Capone,  Karen  A.  Deere, 

Linda  S.  Bangert  and  Paul  S.  Pao 


NASA  Langley  Research  Center 
February  25, 1997 
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HSR  - 1997  Aerodynamic  Performance  Workshop 

Program  Objectives 


• Objectives: 

- For  a representative  HSR  2D  nozzle,  determine  the 
effects  of 

• Nozzle  external  flap  curvature  and  length 

• Sidewall  boattail  angle  and  curvature 

- Develop  an  experimental  data  base  for  2D  nozzles 
with  long  divergent  flaps  and  low  boattail  angles 

- Provide  validation  data  for  isolated  transonic  nozzle 
boattail  drag  program  CFD  prediction  studies 


The  NASA-industry  team  has  sponsored  several  studies  in  the 
last  two  years  to  address  the  installed  nozzle  boattail  drag  issues. 
Some  early  studies  suggested  that  nozzle  boattail  drag  could  be  as 
much  as  25  to  40  percent  of  the  subsonic  cruise.  As  part  of  this  study 
tests  have  been  conducted  at  NASA-Langley  to  determine  the 
uninstalled  drag  characteristics  of  a proposed  nozzle.  The  overall 
objective  was  to  determine  the  effects  of  nozzle  external  flap  curvature 
and  sidewall  boattail  variations.  This  test  would  also  provide  data  for 
validating  CFD  predictions  of  nozzle  boattail  drag. 
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Full-scale  geometry  of  the  “Best  DSM”  nozzle  chosen  for  the 
installed  nozzle  drag  program  is  shown  in  the  upper  left  of  this  figure. 
This  nozzle,  which  comprises  over  40-percent  of  the  total  nacelle  is 
78.90  by  74.88  inches  with  an  aspect  ratio  of  1 .040.  In  order  to  provide 
data  in  a timely  fashion,  tests  were  to  be  conducted  on  an  existing 
propulsion  simulation  system  used  in  the  16-Foot  Transonic  Tunnel.  Of 
the  three  simulators,  the  on  with  the  smallest  aspect  ratio  was  chosen 
for  tests.  This  simulator  is  6.80  by  6.20  inches  with  an  aspect  ratio  of 
1 .096.  From  this,  it  was  decided  to  use  the  height  of  the  model  for 
scaling  purposes.  This  dimension  was  chosen  because  it  gives  the 
best  representation  of  nozzle  boattail  closure.  The  resulting  scale  was 
8.17%.  From  this  scale,  an  appropriate  reference  area  can  be  obtained 
for  subsequent  use  in  nondimensionalizing  drag  in  order  to  produce  a 
meaningful  drag  coefficient  in  terms  of  airplane  drag  counts. 
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HSR  - 1997  Aerodynamic  Performance  Workshop 

Nozzle  Flap  Definition 
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This  chart  summarizes  the  various  parameters  used  to  define  the 
nozzle  external  shape.  The  term  Lf  represents  the  length  of  the 
external  flap  of  the  nozzle.  The  baseline  nozzle  has  a flap  length  of 
8.4233  inches.  Two  nozzles  with  shorter  flap  lengths  were  also  tested. 
Note  that  the  overall  length  of  the  nozzle  was  not  changed  when  nozzle 
flap  length  was  varied.  The  boattail  curvature  parameter  is  rc/rc,max- 
Nozzles  with  curvatures  from  o to  100%  were  tested.  A nozzle  with  no 
curvature  would  probably  be  the  simplest  to  build  since  this  flap  would 
have  a simple  hinge  joint.  All  the  nozzles  had  the  same  internal 
contour,  throat  area  and  exit  area.  Thus  changes  in  performance 
should  only  be  attributed  to  external  flow  effects  over  the  nozzle  flaps. 

The  geometry  of  the  sidewalls  was  defined  in  a similar  manner. 
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Configuration  Matrix 


Nozzle 

Nozzle 

Config 

Flap 

■5BH 

W3SH3B 

K'fl'l 

Height 

m 

F-3 

1.2 

0% 

■HUM 

$-i 

4 

0% 

FuB 

N-2 

F-2 

1.2 

10% 

■IT  " ■ 

S-1 

4 

0% 

Fui 

N-1 

F-1 

1.2 

40% 

S-1 

4 

0% 

Ful 

N-4 

F-4 

1.2 

100% 

S-1 

4 

0% 

Ful 

F-6 

1.0 

10% 

■EM 

S-1 

4 

0% 

Fui 

■3E1 

F-5 

1.0 

40% 

S-1 

4 

0% 

FuB 

1 t±5 

F-1 

1.2 

40% 

■EEHM 

S-2 

6 

0% 

Full 

IK1 

F-1 

1.2 

40% 

■ITn'M 

S-3 

8 

0% 

FuB 

N-7 

F-1 

1.2 

40% 

■m 

S-4 

6 

10% 

Ful 

N-9 

F-1 

1.2 

40% 

■IT  ri  ■ 

S-6 

6 

40% 

Ful 

■si 

F-1 

1.2 

40% 

■llriM 

S-5 

8 

10% 

Ful 

■313 

F-1 

1.2 

40% 

■1  ■ 

S-7 

8 

100% 

Ful 

■3H 

F-2 

1.2 

10% 

mvmm 

S-6 

4 

0% 

None 

■3EH 

F-1 

1.2 

40% 

■ITrrM 

S-6 

4 

0% 

None 

N-1 

PSP 

F-1 

1.2 

40% 

1 6.379 

S-1 

4 

0% 

F-1 

1.2 

40% 

16.379 

S-8 

4 

0% 

ms 

This  table  summarizes  the  various  nozzle  configurations  tested. 
For  the  first  four  configurations,  nozzle  curvature  was  varied  with  the 
nozzle  baseline  length  being  held  constant.  Note  that  boattail  angle 
varies  for  these  configurations  from  1 1 .71 9°  to  23.437°.  Nozzle 
curvature  was  also  varied  for  the  nozzles  with  the  shorter  flap  length  as 
indicated  by  the  fifth  and  sixth  configurations. 

The  effect  of  varying  nozzle  sidewall  boattail  angle  and  curvature 
were  studied  with  the  remaining  configurations.  Sidewall  S-1  with  no 
curvature  and  4°  boattail  was  the  baseline  sidewall.  The  height  of  the 
sidewall  denoted  “full”  was  fixed  to  a distance  such  that  the  nozzle  flaps 
would  not  unport  with  the  nozzle  in  the  supersonic  cruise  position.  In 
this  position  the  nozzle  flaps  have  a boattail  angle  of  0°.  Some  tests 
were  conducted  with  a reduced  height  sidewall  to  determine  what  the 
penalty  might  be  for  having  these  large  sidewalls.  The  height  of  the 
sidewall  labeled  “none”  followed  the  contour  of  nozzle  flap  N1. 
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Experimental  Arrangement  


The  two-dimensional  propulsion  air-powered  simulation  system 
is  shown  in  this  figure.  This  model  is  composed  of  three  major 
components:  a nose-forebody  section,  a centerbody  section,  and  the 
nozzle.  The  nozzle  connect  station  is  at  model  station  50.90.  The 
nose-forebody  was  nonmetric;  that  is  it  was  not  connected  to  the  force 
balance.  This  simulation  system  is  equipped  with  a flow  transfer 
system  that  is  designed  to  minimize  the  transfer  of  axial  momentum 
across  the  force  balance. 

The  nozzle  is  shown  with  the  baseline  sidewalls.  Sidewall  height 
was  fixed  so  that  the  nozzle  flaps  would  not  unport  with  the  nozzle  at 
the  supersonic  cruise  position  (no  nozzle  boattail) 

The  nozzle  upper  or  lower  flaps  were  each  instrumented  with  two 
rows  of  pressure  tapes  with  25  taps  per  row.  Thus,  nozzle  pressure 
drag  was  determined  from  integrating  100  pressures.  On  one  of  the 
nozzle  sidewalls,  there  were  two  rows  of  pressure  taps  with  20  taps  per 
row. 

For  this  investigation,  nozzle  drag  is  defined  over  that  portion  of 
the  model  from  model  station  50.90  to  64.04. 
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CFD  - PAB3D  (version  1 3/ASM) 


• Code  Architecture 

- Modular  multi-block  structure  w/grid  sequencing 

- Multiple-to-one  and  patched  interfaces 

- In-code  calculations  of  integrated  forces,  moments, 
and  flux  quantities 

• Code  Performance  Statistics 

- Compact  memory  requirements:  23  words  per  grid 
point 

- 38  m-sec  (Cray-YMP)  per  iteration  per  grid  point 

- Compatible  with  most  workstations 


PAB3D  solves  the  three-dimensional  Reynolds-averaged  Navier- 
Stokes  equations  with  a finite-volume  formulation  on  structured  multi- 
block grids.  A grid  sequencing  scheme  allows  for  automatic 
assessment  of  grid  density  on  flow  solutions.  In  addition  to  memory 
management,  grid  sequencing  also  allows  for  quick  initial  solutions  and 
increased  convergence  rates. 

One-to-one,  multiple-to-one  and  general  patching  between  block 
interfaces  is  accepted  by  the  flow  solver  for  the  development  of 
complex  geometric  grids.  A conservative  patching  utility  is  used  to 
determine  the  communication  between  block  interfaces. 

The  in-code  performance  package  supplies  the  user  with 
integrated  forces,  moments  and  flux  quantities  in  output  file  formats  that 
are  compatible  with  many  standard  graphics  packages. 
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CFD  - PAB3D  (version  13/ASM) 


• Flow  Solver  Characteristics 

- Mixed  Roe  and  van  Leer  schemes 

- Local  time  stepping  and  upwind  biased  with  user- 
selected  limiter  options 

- Two-equation  k-e  or  ASM  Turbulence  Models 


The  code  allows  the  user  to  select  between  multiple  flow  solvers, 
limiters  and  boundary  conditions  at  code  run  time.  Generally,  the  Roe 
and  van  Leer  schemes  offer  improved  accuracy  and  quick  convergence 
rates. 

The  flow  solver  has  a robust  two-equation  k-e  turbulence  model 
and  several  anisotropic  algebraic  Reynolds  stress  models  (ASM). 
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CFD  Solution  Run-Times 

Nozzle 

Condition 

NPR=5 

Hours 

Cray-YMP 

Grid  | 
Level 

N1 

M=0.9  ! 

16.4 

cut  222 

N1 

M=0.9  | 

46 

base 

(211  sequence)  j 

N1 

M=l.ll  | 

15.9 

cut  222 

N1 

M=1.2 

12.2 

cut  222  j 

N3 

M=0.9 

13.3 

cut  222 

N3 

~ M=1  11 

7.9 

cut  222 

N3 

M=i.2 

7.2 

cut  222 

This  charts  illustrates  the  time  required  on  a Cray-YMP  to  develop  converged 
solutions  with  the  results  in  this  presentation.  This  charts  also  exhibits  the  benefit  of 
database  reduction  and  grid  sequencing. 

For  example,  the  first  two  records  represent  the  solution  of  nozzle  N1  at  M=0.9  and 
ivipp_5  a converged  solution  was  developed  within  16.4  hours  on  a cut  grid.  However, 
the  solution  was  developed  to  the  base  level  to  quantify  the  effect  of  grid  density.  Nozzle 
drag  decrease  a mere  0.2  of  a count  in  another  30  hours  of  Cray-YMP  time.  Since  the 
solution  appears  to  be  minimally  dependent  on  doubling  the  grid  density  beyond  the  cut 
222.  level,  the  remaining  solutions  were  developed  by  sequencing  on  the  cut  grid  only. 

This  allowed  for  quicker  solution  times  due  to  the  substantially  smaller  memory 
requirement. 

The  base  grid  is  a quarter  plane  representation  of  the  experimental  model  with  1 .57 
million  grid  points  in  9 blocks.  Using  a database  reduction  scheme,  a cut  222  grid  is 
generated  by  eliminating  every  other  grid  point  in  the  i,  j,  and  k directions.  This  cuts  each 
grid  dimension  by  2,  which  decreases  the  grid  count  to  207,437  and  substantially  reduces 
the  memory  required  to  run  the  flow  solver.  The  grid  can  be  sequenced  in  each  direction 
for  improved  convergence  rates  and  for  grid  assessment.  For  example,  the  flow  solver 
uses  alternating  points  in  the  i direction  and  every  point  in  the  j and  k directions  in  a 21 1 
sequence.  Generally,  a user  would  begin  sequencing  on  the  cut  222  grid.  A pattern  of 
222,  221 ,211,  and  then  1 1 1 , or  no  sequencing  might  be  used  to  assess  solution  behavior 
as  more  points  are  utilized  in  a particular  direction.  The  solution  is  developed  until 
convergence  requirements  are  met  at  each  level.  Once  the  solution  is  converged  on  the 
cut  222  grid,  the  solution  may  be  extrapolated  to  the  base  grid  and  sequencing  may  again 
be  utilized. 
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Typical  convergence  histories  for  Nozzle  3 are  shown  in  the  following 
three  figures.  Similar  results  were  obtained  for  Nozzle  1 . 

This  figure  shows  the  convergence  history  for  Nozzle  N3  at  M=0.9 
and  NPR=5.  Spikes  in  the  residual  history  exhibit  locations  of  solution 
extrapolation  to  a finer  grid  level.  Accordingly,  small  adjustments  in  drag 
coefficient,  Cd,  are  observed  at  these  locations,  also. 

The  solution  converges  at  the  222  grid  level  after  6500  iterations,  with 
no  further  change  in  drag  coefficient  as  the  residual  continues  to  drop.  Drag 
coefficient  increased  a mere  1 .2  percent  from  the  422  grid  level  to  the  222 
grid  level. 

Computational  drag  coefficient  was  calculated  with  a reference  area 
of  1 in2.  Therefore,  a reference  area  of  A=6824.407  in2  is  used  to  convert  to 
a scale  comparable  to  experimental  results. 
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This  figure  shows  the  convergence  history  for  Nozzle  N3  at 

M=1 .1 1 and  NPR=5.  Spikes  in  the  residual  history  exhibit  locations  of 
solution  extrapolation  to  a finer  grid  level.  Accordingly,  small 
adjustments  in  Cd  are  observed  at  these  locations,  also. 

The  residual  drops  3 orders  of  magnitude  on  the  444  grid  level. 
The  residual  then  flattens  and  Cd  remains  unchanged  after  500 
iterations,  suggesting  convergence  at  this  grid  level.  The  solution 
converges  at  the  222  grid  level  after  4500  iterations.  Drag  coefficient 
decreased  2.3  percent  from  the  422  grid  level  to  the  222  grid  level. 
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This  figure  shows  the  convergence  history  for  Nozzle  N3  at 

M=1 .2  and  NPR=5.  Since  minimal  change  in  drag  coefficient  was 
apparent  from  the  previous  solutions,  this  solution  was  developed  on 
the  222  grid  level  in  an  effort  to  minimize  resources. 

The  residual  drops  3 orders  of  magnitude  and  flattens  after  1 000 
iterations.  Drag  coefficient  remains  constant  after  1500  iterations. 

This  suggests  solution  convergence  at  this  grid  level. 
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Predicted  and  Experimental  Drag  Comparison 


Flap  1:  U /hm  = 1.2,  (if  = 16.38°,  rc/rc,max=  0.4 
Sidewall  1:  P*w  = 16.38°,  rc/rc,max  = 0 

NPR  = 5 


M=  0.9  Ms  1.11  M = 1.2 


Nozzle  N1  has  a 40  percent  flap  curvature  along  the 
boattail.  The  predicted  drag  is  in  excellent  agreement  with 
experimental  drag  at  M=0.9  and  is  within  0.7  of  a count  of  experimental 
drag  at  M=1 .2.  PAB3D  predicted  a nozzle  drag  of  1 1 .21  counts  at  M= 
1.11. 
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^Predicted  and  Experimental  Pressures  Along  Flap 


As  expected  from  the  agreement  between  nozzle  drag,  Cd(EFD)= 
9.8  counts  compared  to  Cd(CFD)=9.53  counts,  the  predicted  pressures 
are  in  excellent  agreement  with  experimental  data  at  M = 0.9  and  NPR 
= 5. 

Data  is  compared  along  the  flap  centerline  (z=0),  at  a flap 
outboard  station  (z=2.61  in)  and  at  a flap  ‘mid’  station  (z=1 .373  in)  in 
this  and  the  following  pressure  charts. 
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Predicted  and  Experimental  Pressures  Along  Flap 


Msl.11,  NPR  = 5 Nozzle  N1  - Flap  1 , Sidewall  1 


Predicted  pressures  exhibit  the  expansion  and  shock  along  the 
flap  at  three  stations.  Pressure  differences  between  the  stations 
suggest  highly  three-dimensional  flow  along  the  flap. 
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_F>redicted  and  Experimental  Pressures  Along  Flap 


Predicted  pressures  are  in  good  agreement  with  experimental 
data  at  M = 1 .2  and  NPR  = 5.  The  largest  difference  in  pressure  occurs 
near  x/L=0.8  at  the  outboard  and  mid  stations. 
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Nozzle  N3  has  a sharp  comer  (no  curvature)  leading  into  the 
flap.  The  predicted  drag  is  within  0.3  of  a count  at  M=1 .2  and  within  1 .6 
counts  at  M=0.9.  PAB3D  predicted  a nozzle  drag  of  9.5  counts  at  M= 
1.11,  which  is  an  improvement  compared  to  Nozzle  N1 . 
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Predicted  and  Experimental  Pressures  Along  Flap 


M = 0.9,  NPR  = 5 


The  pressure  recovery  along  the  flap  is  underpredicted 
compared  with  experimental  data  at  M = 0.9  and  NPR  = 5.  The  higher 
experimental  pressure  recovery  results  in  lower  drag  compared  with 
predicted  data,  Cd(EFD)=1 .97  counts  compared  to  Cd(CFD)=3.56  counts. 
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The  pressure  recovery  along  the  flap  corresponds  to  a drag  of 
^d(CFD)=®'^  counts. 
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Predicted  and  Experimental  Pressures  Along  Flap 


M = 1.2,  NPR  = 5 


In  general,  the  predicted  pressures  are  in  good  agreement  with 
experimental  data  at  M = 1 .2  and  NPR  = 5,  which  is  expected  with  the 
good  correlation  between  nozzle  drag  also,  Cd(EFD)=9.8  counts 
compared  to  Cd^p0)=9.53  counts. 
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Comparison  of  Pressure  Drag 


NPR  = 5 


12 


Nozzle  Nozzle  Nozzle  Nozzle 

N1  N1  N3  N3 

M=.9  M=1 .2  M=.9  M=1 .2 


m Integrated 
Pressure  Drag 

□ Balance  Measured 
Pressure  Drag 

□ CFD-PAB3D 
Pressure  Drag 


This  figure  compares  pressure  drag  for  each  nozzle  at  M=0.9 
and  M=1 .2,  in  an  effort  to  understand  the  larger  difference  between 
predicted  and  experimental  drag  for  Nozzle  N3  at  M=0.9,  compared 
with  the  other  solutions.  Three  methods  were  used  for  determining 
pressure  drag:  integrated  pressures,  balanced  measured  and 
computational  fluid  dynamics  using  PAB3D.  It  appears  as  though  there 
is  a discrepancy  using  the  integrated  method  compared  with  balance 
measured  and  PAB3D  predicted  pressure  drag  methods.  Work  is  in 
progress  to  discern  the  data. 
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Nozzle  Pressure  Drag  Breakdown  - Nozzle  N1 


NASA  Langley  16-Ft  Transonic  Tunnel  Test  477 

Flap  1 : U /hm  = 1 -2,  bf  = 16.38°  , rc/rc,max  = 0.4 
Sidewall  1i  b$w  — 16.38°,  Tc  /rc,max  = 0 


b Sidewall 
sRap 
b Total 


This  figure  illustrates  the  breakdown  of  pressure  drag  between 
the  nozzle  flaps  and  the  sidewalls.  Also  shown  is  the  effect  of  the  jet  on 
nozzle  pressure  drag.  Note  that  the  jet  effects  shown  are  typical  for  this 
class  of  nozzles  and  the  breakdown  of  the  pressure  drags  are  similar 
for  the  other  nozzle  tested. 
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Flap  Curvature  Effects  on  Flap  Pressure  Drag 

NASA  Langley  16-Ft  Transonic  Tunnel  Test  477 
Nozzle  Flap  Lf  / hm  = 1.2 


NPR  = 5 


M = 0.90  M = 0.95  M = 1 .20 


The  effect  of  varying  nozzle  boattail  curvature  on  the  nozzle  flap 
pressure  drag  is  illustrated  above  for  the  nozzle  with  the  baseline  flap 
length.  As  can  be  seen,  the  nozzle  N3  with  no  curvature  had  the  least 
flap  pressure  drag.  This  result  was  somewhat  surprising  because 
previous  experience  has  shown  that  nozzles  with  a sharp  shoulder 
generally  had  higher  drag.  However,  these  nozzles  had  shorter 
external  flaps  and  higher  boattail  angles  (usually  greater  than  16°).  In 
order  to  try  to  understand  this  result,  pressure  distributions  on  these 
nozzles  will  be  shown  in  some  subsequent  figures. 
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Flap  Curvature  Effects  on  Nozzle  Drag 

NASA  Langley  IS-Ft  Transonic  Tunnel  Test  477 


Nozzle  Flap  L,/  hm  = 1.2,  NPR  = 5 
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Similar  results  from  the  previous  figure  are  also  shown  above. 
Total  nozzle  drag  is  the  sum  of  the  flap  and  sidewall  pressure  drags 
plus  a skin  friction  drag  (calculated  as  a simple  flat  plate  friction  drag). 
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Flap  Curvature  Effects  on  Nozzle  Drag 

NASA  Langley  1 6-Ft  Transonic  T unnel  Test  477 
Nozzle  Flap  L, / hm  = 1 .2,  NPR  = 5 


Nozzle  tjxcmax  ft,  deg 


M = 0.90  M = 0.95  M = 1 .20 


Similar  results  are  obtained  when  the  total  nozzle  drag  measured 
by  the  force  balance  is  considered. 
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Pressure  distributions  along  the  center  row  of  the  top  flap  for 
nozzle  N1 , N2,  N3,  and  N4  are  presented  at  M = 0.9  and  at  a nozzle 
pressure  ratio  of  5.  An  NPR  of  5 was  used  in  the  CFD  installed  nozzle 
study.  Basically  what  is  shown  is  that  even  though  nozzle  N3  with  no 
curvature  had  the  greatest  expansion  of  flow  about  the  nozzle  shoulder, 
it  had  better  pressure  recovery  characteristics  than  the  other  three 
nozzles. 

It  should  be  pointed  out  that  these  pressure  distributions  do  not 
show  the  presence  of  shocks  on  the  nozzle  flap.  This  result  is  typical  for 
most  of  the  nozzles  tested.  This  apparent  lack  of  shocks  on  the  nozzle 
probably  results  in  the  low  boattail  angles  that  most  of  the  nozzles 
have. 
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Similar  results  to  those  previously  show  were  also  obtained  at 
M = 0.95. 
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Similar  results  to  those  previously  show  were  also  obtained  at 
M = 1 .20. 
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lap  Curvature  Effects  for  Nozzle  With  Short  Flap 


NASA  Langley  16-Ft  Transonic  Tunnel  Test  477 


Nozzle  Flap  Lf/hm  = 1.0 


NPR  = 5 


10 


0 

M = 0.90  M = 0.95  M = 1 .20 


This  figure  shows  flap  pressure  drag  for  the  two  nozzles  tested 
with  the  smaller  flap  lengths.  In  contrast  to  previous  results,  lower  flap 
pressure  drag  was  obtained  with  nozzle  N5  which  had  more  curvature. 
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Nozzle  Flap  Length  Effects 


NASA  Langley  16-Ft  Transonic  Tunnel  Test  477 
Nozzle  Rap  rc/rc  m«=  0.1 


The  effect  of  nozzle  flap  length  is  illustrated  in  this  figure  where 
flap  pressure  drag  is  compared  between  nozzle  N2  and  N6.  Each  of 
these  nozzles  had  a curvature  ratio  of  0.1 . As  can  be  seen,  the  nozzle 
with  the  longer  nozzle  flap  had  lower  drag  at  all  the  the  Mach  numbers 
tested.  This  was  probably  due  to  the  lower  boattail  angle  N2  had. 
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Nozzle  Flap  Length  Effects 


NASA  Langley  16-Ft  Transonic  Tunnel  Test  477 
Nozzle  Rap  rc/rcmax  = 0.4 

NPR  = 5 


M = 0.90  M = 0.95  M = 1 .20 


Results  similar  to  those  shown  in  the  previous  figure  were 
obtained  for  nozzles  N1  and  N5.  These  nozzles  had  a curvature  ratio 
of  0.4. 
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HSR  - 1997  Aerodynamic  Performance  Workshop 

Nozzle  Sidewall  Boattail  Angle  Effects 


NASA  Langley  16-Ft  Transonic  Tunnel  Test  477 
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This  figure  summarizes  the  effect  of  sidewall  boattail  angle  and 
curvature  on  sidewall  pressure  drag.  At  M = 0.90,  nozzle  N6  with  a 8° 
boattail  angle  had  the  lowest  sidewall  pressure  drag.  At  M = 1 .2,  just 
the  opposite  was  true.  At  this  mach  number,  nozzle  N6  had  about  2.3 
times  as  much  drag  as  nozzle  N1  with  the  4°  boattail.  Although  the  4° 
boattail  sidewall  was  considered  to  be  the  baseline,  it  is  probable  that 
the  boattail  angle  will  have  to  ber  as  much  as  8°  in  order  to  house 
nozzle  actuation  hardware. 

Putting  full  curvature  on  the  sidewall  with  the  8°  boattail  resulted 
in  a one  count  pressure  drag  reduction.  Sidewalls  with  full  curvature 
are  feasible  for  the  full  scale  aircraft  since  the  sidewalls  are  fixed. 
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Pressure  distributions  along  the  center  row  of  the  sidewall  for 
nozzle  N1 , N5,  N6,  and  N10  are  presented  at  M = 0.9  and  at  a nozzle 
pressure  ratio  of  5.  These  results  are  similar  to  thsose  already  shown 
for  the  nozzle  flaps.  Even  though  nozzle  N6  with  the  8°  sidewall 
boattail  had  the  greatest  expansion  of  flow  about  the  nozzle  shoulder,  it 
had  excellent  pressure  recovery  characteristics  such  that  sidewall  drag 
was  lower. 
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Similar  results  to  those  previously  show  were  also  obtained  at 
M = 0.95 
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At  M = 1 .2,  it  appears  as  though  pressures  are  lower  on  nozzle 
N6  with  the  8°  sidewall  boattail.  This  lower  pressure  is  acting  over 
most  of  the  sidewall  and  causing  the  higher  drag. 
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HSR  - 1997  Aerodynamic  Performance  Workshop 

Nozzle  Sidewall  Curvature  Effects 


NASA  Langley  16-Ft  Transonic  Tunnel  Test  477 
Nozzle  Sidewall  ps=  8° 


NPR  as  5 


M = 0.90  M = 0.95  M = 1.11  M = 1.20 


This  figure  illustrates  the  effect  of  sidewall  curvature  for  those 
sidewalls  with  8°  boattail  angle.  At  all  Mach  numbers,  the  sidewall  with 
a curvature  ratio  of  1 had  the  lowest  sidewall  pressure  drag. 
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HSR  - 1997  Aerodynamic  Performance  Workshop 

Nozzle  Sidewall  Height  Effects 

NASA  Langley  16-Ft  Transonic  Tunnel  Test  477 
Nozzle  Sidewall  P8  = 4°,  rc  /rc  max  = 0 


M = 0.9,  NPR  = 3.54 M = 1 .20,  NPR  = 6.65 


It  was  previously  stated  that  the  height  of  the  sidewall  was  fixed 
to  a distance  such  that  the  nozzle  flaps  would  not  unport  with  the 
nozzle  in  the  supersonic  cruise  position.  In  this  position  the  nozzle 
flaps  have  a boattail  angle  of  0°.  Some  tests  were  conducted  with  a 
reduced  height  sidewall  just  to  see  what  the  penalty  might  be  for  having 
these  large  sidewalls.  The  reduced  height  sidewall  followed  the  contour 
of  nozzle  flap  N1 . As  can  be  seen  in  the  above  figure,  not  only  was 
sidewall  pressure  drag  dramatically  reduced,  so  was  the  flap  pressure 
drag. 
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HSR  - 1 997  Aerodynamic  Performance  Workshop 

Summary 


• All  nozzles  exhibited  expected  pressure  distributions. 
At  most  test  condition,  there  were  no  shocks  or  shock 
induced  separation 

• Excellent  corelation  between  experimental  and  CFD 
results  were  obtained 

• The  nozzle  with  an  external  flap  with  a sharp 
shoulder  (no  curvature)  had  the  least  pressure  drag 

• For  the  nozzle  with  a sidewall  with  8°  boattail  angle, 
pressure  drag  increased  from  0.3  to  3.5  counts  from 
M = 0.9  to  1.2. 
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Abstract: 

In  early  1996  the  NASA-Industry  High  Speed  Research  Technical  Integration 
team  released  the  final  definition  of  the  HSCT  Technology  Concept  Airplane 
(TCA).  This  configuration  represents  the  integration  of  current  inputs  from  all 
technical  disciplines  into  a realistic  High  Speed  Civil  Transport  concept.  This 
paper  reviews  the  development  and  content  of  the  high  speed  aerodynamics 
inputs  to  the  TCA  sizing  and  flight  performance  predictions.  The  paper  also 
summarizes  subsequent  detailed  analysis  work,  CFD,  and  TCA  wind  tunnel 
test  data  that  are  now  being  used  to  assess  the  drag  levels  of  the  "status" 
airplane  (i.e.  without  projections).  A bottoms-up  assessment  of  the  high 
speed  drag  technology  projection  is  shown  to  identify  reasonable  sources  of 
drag  improvements  that  would  meet  the  target  levels.  Sources  of  uncertainty 
in  the  current  HSCT  high  speed  drag  predictions  are  outlined,  and  areas  for 
risk  reduction  in  future  performance  predictions  are  identified. 


* Lead  Engineer 

**Senior  Engineer 

HSCT  High  Speed  Aerodynamics 

Configuration  Development  Group 
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TCA  Flight  Drag  Polars  For  Airplane  Performance 


• High  speed  performance  drag  polar  build-up  process 

• Relative  magnitude  of  HSCT  drag  and  lift  components 

• Performance  polar  technology  projection  assumptions 

• TCA  “Status”  drag  and  uncertainty  sources 

• Summary 


Introduction: 

In  early  1996  the  NASA-Industry  High  Speed  Research  Technical  Integration 
team  released  the  final  definition  of  the  HSCT  Technology  Concept  Airplane 
(TCA).  This  configuration  represents  the  integration  of  current  inputs  from  all 
technical  disciplines  into  a realistic  High  Speed  Civil  Transport  concept.  The 
TCA  replaces  the  Boeing  " Reference- H"  and  McDonnell-Douglas  "M2.4-7A" 
designs  used  previously  as  parallel  HSR  baselines  (primary  and  alternate, 
respectively).  The  certainty  with  which  the  TCA’s  external  lines  and  aerodynamic 
parameters  correctly  reflect  the  characteristics  of  a production  HSCT  is  critical  for 
assessing  the  continued  viability  of  the  program,  and  for  correctly  down-selecting 
long  lead-time  items  such  as  propulsion  components. 

This  paper  reviews  the  development  and  content  of  the  high  speed 
aerodynamics  inputs  to  the  TCA  sizing  and  flight  performance  predictions.  This 
discussion  includes  a review  of  the  "technology  projection"  assumptions  and  a 
bottoms-up  assessment  of  potential  drag  improvements  required  to  meet  the 
target  levels.  Subsequent  detailed  analysis  work,  CFD,  and  TCA  wind  tunnel 
test  data  that  are  now  being  used  to  assess  the  drag  levels  of  the  "status" 
airplane  (i.e.  without  projections)  are  summarized.  Sources  of  uncertainty  in  the 
current  HSCT  high  speed  drag  predictions  are  outlined,  and  areas  for  risk 
reduction  in  future  performance  predictions  are  identified. 


709 


The  process  followed  in  building  up  the  high  speed  aerodynamics  inputs  for  the 
TCA  was  essentially  identical  to  that  used  for  airplane  level  Tl  trades.  This 
process  anchors  the  drag  polars  to  available  wind  tunnel  test  data  while 
accounting  for  increments  between  the  wind  tunnel  model  and  a full  scale 
airplane  in  flight.  These  increments  include  all  configuration  differences,  skin 
friction  changes  due  to  Reynolds  number,  trim  and  excrescence  drags,  and  a 
projection  of  future  drag  reduction.  It  can  be  seen  by  rearranging  the  terms  in 
the  first  equation  above  that  the  geometry  related  increment  is  equal  to  the 
computed  drag  of  the  full  TCA  plus  the  difference  between  the  measured  and 
computed  results  for  the  reference  test  model.  The  wind  tunnel  data  for  the  TCA 
polar  build-up  was  from  the  “Reference-H”  model.  The  skin  friction  at  flight 
Reynolds  number  was  computed  using  a strip-wise  summation  of  the  Karman- 
Schoenherr  flat  plate  skin  friction  with  the  Sommer-Short  compressibility 
correction,  assuming  an  adiabatic  wall.  An  allowance  was  made  for  excrescence 
("protuberance  drag")  approximating  that  of  current  subsonic  transports  at 
incompressible  speeds  (7%  of  skin  friction).  This  incompressible  value  was 
scaled  versus  Mach  to  account  for  compressibility.  The  high  speed  polars  used  a 
normalized  trim  drag  allowance  vs.  Mach  which  was  based  on  previous  detailed 
analyses  of  the  Ref.-H  and  B-2707.  Technology  projections  were  added  based 
on  a tops-down  L/D  assessment.  CFD  was  employed  as  a "numerical  wind 
tunnel"  on  the  final  geometry  to  confirm  that  additional  drag  had  not  inadvertently 
been  incurred  in  "real  flow"  during  lofting  of  the  linear  design.  (This  generally 
results  in  CFD-based  adjustment  of  less  than  +/- 1.5  counts  at  Mach  2.4.) 
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The  chart  above  illustrates  the  relative  magnitudes  of  the  various  components  of 
the  resulting  TCA  drag  polars  used  in  the  configuration  sizing  and  mission 
performance  analysis  for  the  Mach  2.4  cruise  condition. 


Several  facts  are  immediately  noted.  One  is  that  the  majority  of  the  cruise  drag 
is  in  the  areas  of  skin  friction  and  drag  due  to  lift.  Roughly  one  third  (1 6 counts) 
of  the  drag  due  to  lift  is  due  to  ideal  incompressible  induced  drag.  The  remainder 
is  vortex  drag  (span  load  efficiency),  and  wave  drag  due  to  lift.  While  non-linear 
optimization  can  significantly  reduce  the  drag  due  to  lift  term,  the  substantial  skin 
friction  level  is  a result  of  configuration  wetted  area  and  cannot  be  significantly 
reduced  without  the  application  of  laminar  flow  technologies  which  are  not 
included  in  the  TCA  definition.  The  large  size  of  the  skin  friction  term  relative  to 
the  total  high  speed  aerodynamics  technology  projection  means  that  accurately 
predicting  the  skin  friction  in  flight  at  combined  high  Reynolds  number /Mach 
number  conditions  is  of  enough  importance  to  warrant  more  attention.  That  is  to 
say  for  example  that  the  potential  drag  error  incurred  by  a 4%  uncertainty  in  the 
skin  friction  component  alone  is  equal  to  a 30%  shortfall  in  the  gain  projected  for 
direct  wing-body  non-linear  optimization.  Differences  between  the  various  flat 
plate  skin  friction  handbook  calculation  methods,  various  Navier-Stokes  CFD 
solutions,  and  wind  tunnel  measured  Cdmin  levels  easily  cover  a range  of  +/-  4 
to  5%  of  the  friction  drag. 


Another  fact  to  take  note  of  is  that  the  successful  nacelle  installation  design  of 
the  TCA  results  in  favorable  nacelle  pressure  field  wave  and  lift  interference 
nearly  canceling  the  nacelle’s  isolated  wave  drag. 
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TCA  High  Speed  Drag  Components 
Mach  1.1,  Ref.  CL  = 0.16 


The  above  charts  illustrates  the  equivalent  drag  component  break-out  for  the 
typical  transonic  acceleration  condition  at  Mach  1.1,  CL=0.16.  The  large  rise  in 
the  zero-lift  wave  drag  component  near  Mach  one  is  evident.  There  is  a 
corresponding  increase  in  the  excrescence  drag  level  and  the  nacelle  wave  drag. 
The  favorable  nacelle  interference  is  only  slightly  greater  than  the  supersonic 
cruise  value  so  it  no  longer  calcels  out  the  isolated  wave  drag  of  the  nacelles. 

The  technology  projection  near  Mach  1 is  scaled  up  to  remain  the  same 
percentage  improvements  in  wave  drag  and  drag  due  to  lift  as  in  the  supersonic 
cruise  tops-down  analysis.  While  the  resulting  projection  level  is  very  large 
(about  18  counts)  it  is  believed  that  a strong  favorable  spillage/bypass 
interference  term  will  be  able  to  make  up  nearly  half  of  this  amount.  This 
assumption  is  based  on  test  data  from  sting-mounted  nacelles  under  a Ref.H 
model  wing  in  the  ARC  1 1 Ft  tunnel.  Additional  thrust-drag  accounting  details, 
CFD  validation,  and  spillage/bypass  test  data  with  captive  nacelles  (l.e.  with 
wing-mount  and  diverters  present)  will  be  required  before  the  potential  large 
favorable  effects  are  understood  well  enough  for  them  to  be  bookkept  as  part  of 
the  “status”  flight  polar.  Another  potentially  favorable  Propulsion  Induced  Effect 
(“PIE”)  that  is  not  well  enough  quantified  to  justify  inclusion  in  the  Tl  polar  build- 
up is  the  nozzle  boattail  and  exhaust  plume  interference  pressures  on  the  closing 
aft  fuselage.  For  the  time  being  all  of  these  PIE  effects  are  lumped  in  with  the 
transonic  drag  technology  projection. 
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The  equivalent  drag  components  for  the  TCA’s  subsonic  “over-land  cruise” 
condition  at  Mach  0.9  are  shown  above.  The  reference  CL  of  0.18  was  chosen 
to  approximate  the  average  of  heavy  and  light  airplane  subsonic  cruise  legs, 
climbs,  and  descents.  The  zero  lift  wave  drag  terms,  nacelle  wave  and  lift 
interference  terms  and  PIE  terms  are  all  assumed  equal  to  zero  at  Mach  0.9. 
Nacelle  on/off  wind  tunnel  comparison  data  on  previous  configurations  has 
shown  that  a favorable  increment  in  lift  at  alpha  would  probably  be  present  on 
the  TCA  but  that  it  may  be  accompanied  by  a corresponding  increase  in  pressure 
drag  so  there  is  no  net  favorable  nacelle  interference  credited  to  the  drag  polar. 

Theoretically,  AER02S  or  similar  panel  method  code  could  be  used  to  obtain 
preliminary  design  (“PD”)  estimates  of  the  drag  due  to  lift  at  these  conditions,  but 
the  leading  edge  suction  levels  predicted  by  such  codes  tend  to  be  less  accurate 
at  locally  compressible  Mach  numbers,  especially  with  the  significant  bluntness 
of  the  TCA  inboard  airfoils.  In  lieu  of  a robust  calibrated  PD  drag  code  for  the  0.9 
to  1 .0  Mach  range,  drag  polars  for  the  TCA  were  created  using  Ref.-H’s  NTF 
wind  tunnel  data  adjusted  for  aspect  ratio  and  wetted  area  differences.  The  Ref.- 
H dragrise  shape  was  also  assumed  up  to  Mach  1.1,  beginning  at  around  Mach 
0.95-0.98.  The  TCA  polars  were  not  penalized  for  the  more  triangular  span  load 
of  the  TCA  relative  to  the  Ref.-H,  nor  were  any  possible  adverse  off-design 
aeroelastic  deflections  taken  into  account.  Leading  and  trailing  edge  flap 
effectiveness  equal  to  that  of  the  Ref.-H  was  also  assumed. 
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The  relative  proportion  of  drag  due  to  each  major  physical  component  is  shown 
above  at  all  three  Mach  numbers.  The  drop  in  skin  friction  and  increase  in  wave 
drag  with  increasing  Mach  is  obvious  when  plotted  in  this  manner  as  percents  of 

the  total  drag. 
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During  the  preparations  for  performing  non-linear  CFD  optimizations  on  the  TCA 
baseline  geometry,  the  choice  of  a reference  CL  as  the  “cruise  point”  for  the 
optimization  was  the  topic  of  considerable  discussion.  Numerous  questions  were 
raised  about  the  expected  differences  between  total  configuration  lift,  the  wing- 
body-nacelle  lift,  and  the  lift  contribution  of  the  TCA  wing-body  alone.  The  chart 
above  illustrates  the  large  amounts  of  lift  that  may  be  accounted  for  by  sources 
other  than  the  wing-body.  The  performance  and  size  optimization  program  used 
to  select  the  final  wing  area  and  engine  thrust  of  the  TCA  was  not  at  the  time 
coded  to  account  for  the  effects  of  geopotential  altitude,  orbital  velocity,  and  the 
downward  component  of  the  nozzle  gross  thrust  vector.  The  effective  “lift” 
contributions  of  these  terms,  while  small,  are  not  negligible  and  have  been 
included  in  the  chart  above. 
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HSCT  High  Speed  Aero  Technology  Projection 


Objective,  Tops-Down  Method, 
Doesn’t  Depend  On  “Status” 

“Optimum”  Is  Linear  Theory 
Physics  Based 

Configuration  Dependent  L/D, 
Showing  CDo  and  CDL  Terms 

“Target”  is  % of  Theory  Judged 
Achievable  With... 

- Fixed  General  Arrangement 

- MDO  Geom.  Constraints 


The  “aerodynamic  technology  projection”  is  a key  part  of  the  polars.  At 
supersonic  speeds  the  projection  is  split  between  zero  lift  drag  (CDo)  and  drag 
due  to  lift  (CDL)  improvements.  A fixed  projection  of  roughly  6%  L/D  is  applied 
subsonically,  and  an  excrescence  improvement  of  20%  is  assumed  at  all  Mach 
numbers.  The  supersonic  cruise  projection  was  formulated  based  on  a tops- 
down  method  (by  R.  M.  Kulfan)  that  uses  the  overall  geometric  parameters  of  a 
given  configuration  to  determine  a theoretical  "best  possible  LyD"  assuming 
complete  freedom  to  optimize  the  aerodynamic  shape  within  the  bounds  of  the 
wing-body  general  arrangement.  This  theoretical  optimum  is  then  adjusted 
downward  to  account  for  the  fact  that  additional  multi-disciplinary  geometry 
constraints  will  be  required  and  that  some  shape  changes  which  could  improve 
cruise  drag  will  be  negated  by  off-design  penalties.  With  input  from  the  HSR 
Configuration  Aerodynamics  (CA)  team,  the  TCA's  projected  performance  target 
level  was  fixed  at  a consensus  level  of  about  95%  of  the  theoretical  optimum. 

It  should  be  noted  that  using  this  approach  the  total  projected  performance  is 
defined  as  a specific  L/D  level  that  does  not  depend  on  knowing  the  performance 
of  the  "status"  (unprojected)  airplane-only  its  gross  geometric  parameters  (span, 
length,  volume,  cross-sectional  area,  nacelle  size,  etc.).  The  size  of  the 
"projection  increment"  itself  is  simply  the  difference  between  the  "target 
performance"  level  and  the  current  best  assessment  of  the  airplane's  status  at 
full  scale  flight  conditions. 
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TCA  High  Speed  Drag  Technology  Projection 
"Bottoms-Up  Assesment" 
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In  order  to  keep  a constant  check  on  the  assumptions  being  made  in  the 
committed  TCA  performance  levels,  Tl  has  maintained  a “bottoms-up”  check  of 
the  reasonableness  of  the  tops-down  projections.  The  projections  for  the  three 
most  critical  high  speed  Mach  numbers  are  shown  above.  Each  bar  is  divided 
into  identifiable  sources  of  potential  drag  improvements  that  would  sum  up  to  the 
total  projection  at  each  Mach.  In  recent  months,  the  CFD  non-linear  optimization 
tools  appear  to  be  well  on  their  way  to  providing  the  6.5  counts  improvement 
assumed  above  at  the  Mach  2.4  cruise  design  point.  It  is  believed  that  some 
combination  of  viscous  and  inviscid  non-linear  optimization,  analysis,  and  flow 
visualization  tools  will  enable  additional  improvements  of  1 .5  to  2.5  counts  in 
detail  design  at  Mach  2.4.  CFD-aided  detail  design  integration  and  multi-point 
design  of  body  area  rule,  nacelle  shapes  and  wing  camber  are  believed  to 
account  for  more  than  a third  of  the  Mach  1 .1  projection.  As  mentioned  earlier, 
the  other  large  piece  of  the  Mach  1.1  projection  will  be  in  understanding  the  PIE 
effects  well  enough  to  reliably  take  credit  for  favorable  transonic  interference 
terms  (or  eliminating  adverse  effects)  wherever  possible. 


Accounting  for  thrust  and  PIE  effects,  and  optimizing  the  wing-body,  nozzle,  and 
tail  contributions  to  trim  drag  at  Mach  2.4  could  account  for  2 to  3 more  counts  of 
the  cruise  projection.  Improved  validation  of  subsonic  CFD  methods  for  subsonic 
Reynolds  number  effects,  leading  edge  flap  optimization  and  trim  including 
aeroelastics,  are  assumed  to  account  for  most  of  the  subsonic  cruise  projection. 
The  projection  of  a 20%  reduction  in  excrescence  at  all  Machs  appears  to  be  a 
reasonable  goal  making  up  the  remainder  of  the  projections. 
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HSR  Technology  Integration 
‘Airplane  Design  Process”  Study 
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Based  on  experiences  with  the  downselect  to  the  TCA  configuration  in  1 995, 
it  was  recognized  that  the  HSR  program  would  benefit  from  improved  PD 
prediction  methods.  A Tl  “Airplane  Design  Process”  (ADP)  task  was 
implemented,  directed  toward  improving  the  Propulsion,  Noise,  Weights,  and 
Aerodynamics  groups’  inputs  to  airplane  level  trade  studies,  metrics,  and 
annual  baseline  updates.  Aerodynamics  discipline  tasks  under  this  on-going 
effort  were  divided  into; 

1)  understanding  differences  between  (and  improving)  the  various  linearized 
potential  flow  design  and  analysis  codes  used  by  NASA  and  industry,  and 

2)  improving  the  accuracy  of  the  build-up  from  the  potential  flow  analysis 
results  to  the  full  scale  flight  drag  “status”  polars. 

For  the  high  speed  performance  regime  (Mach  >0.6)  responsibility  for  the 
former  task  was  given  to  the  McDonnell-Douglas,  while  Boeing  (BCAG) 
undertook  the  latter,  with  NASA  LaRC  contributing  to  understanding  both 
areas.  The  BCAG  effort  includes  sub-tasks  addressing  excrescence  drag 
assessment,  aft-body  / trim  drag,  transonic  flap  predictions,  and  Tl’s  use  of 
corrected  wind  tunnel  data  in  the  drag  polar  build-up  process.  The  potential 
flow  analysis  methods,  excrescence  drag,  and  an  alternate  approach  to  wind 
tunnel  trip  drag  corrections  are  the  subjects  of  other  papers  in  this  Workshop. 
ADP  work  planned  for  1997  includes  updating  Tl’s  performance  polars  for 
CA’s  technology  projections  and  TCA  wind  tunnel  test  data,  as  well  as 
improving  capability  for  drag  predictions  of  multi-surface  control  concepts. 
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The  “ADP”  task  must  take  an  end-to-end  view  of  the  performance  prediction 
process  in  order  to  identify  sources  of  uncertainty  and  prediction  improvements. 
Accurate  aerodynamic  performance  prediction  is  important  for  all  types  of 
aircraft.  For  most  military  aircraft  and  civil  light  or  utility  aircraft,  the  most  critical 
performance  parameters  are  often  those  not  directly  dependent  on  cruise  drag 
(e.g.  thrust-to-weight  ratio,  turn  capability,  takeoff  distance,  roll  rate,  maximum 
payload,  or  range-at-any-cost).  On  long-range  commercial  transports  where  the 
economic  margin  for  error  is  very  small,  the  ability  to  very  accurately  predict 
absolute  cruise  performance  is  essential  for  success.  The  cruise  performance 
parameters  that  directly  impact  the  commercial  viability  of  a transport  are  the 
payload  carried,  the  distance  flown  in  a given  time,  and  the  fuel  burned  to  get 
there  — all  of  which  are  relatively  easy  to  measure  in  flight  on  a given  airplane  to 
within  1%  or  less.  To  sell  a new  airplane  type,  the  manufacturer  must  typically 
guarantee  that  a deviation  of  no  more  than  a few  percent  will  exist  between  the 
predicted  performance  and  that  measured  in  flight  tests.  Unfortunately,  the  pre- 
flight predictions  can  only  be  calculated  from  the  constituent  lift,  drag,  thrust, 

SFC,  and  mass  properties,  each  of  which  have  their  own  uncertainties 
(represented  by  Probability  Density  Functions,  PDF).  The  ability  of  the  Aero  and 
Propulsion  disciplines  to  accurately  establish  their  components  of  performance 
can  only  be  assessed  in  “hindsighf  by  back-calculating  the  corresponding 
components  from  the  measurable  flight  test  quantities.  (Physical  understanding, 
testing  and  CFD  technologies  and  accounting  methods  can  then  be  developed  to 
minimize  unpleasant  surprises  in  the  future.) 
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High  Speed  Database  Uncertainty  Elements 


Data  Components: 


Basic  Truncated 
Wing/Body/Nac/Div 
Forces  & Moments 


Evaluation  of  Full  Scale  ^ 
Airplane  Effects 

• Propulsion  induced  effects  “PIE” 
(Jet/Bleed/Bypass/Spillage/Boattaill 

• Control  surface  drag 

• Flight  Rn  skin  friction 


Corrections  Full  Scale 
Airplane  Geometry 

• Aftbody  closure/up wseep 

• Tail/Canard  trim  drag 

• Excrescence  drag 


While  most  aerodynamic  technology  development  and  optimization  work  can 
focus  on  incremental  CFD  and  test  databases,  the  importance  of  establishing 
accurate  absolutes  as  early  as  possible  is  obvious.  The  net  impact  of  every  one 
count  cruise  drag  miss  (A  CD=0.0001 ) integrated  over  a complete  flight  results  in 
a 10,000  LB  increase  in  the  design  Maximum  Takeoff  Gross  Weight  (MTOW) 
required  to  fly  the  mission.  Once  the  engine  size  is  frozen  and  MTOW  is  fixed 
the  airplane  must  loose  6 seats,  or  find  a one  ton  savings  in  structural  weight  just 
to  break-even  with  a one  count  cruise  drag  increase. 


At  this  early  stage  of  the  HSCT  program,  there  are  really  two  main  areas  where 
integrated  assessments  of  “absolute”  drag  levels  are  needed.  One  is  in  defining 
the  conceptual  airplane  used  as  a confirmation  of  continued  program  viability  and 
a common  baseline  airplane  geometry  for  all  disciplines  to  use.  (The  TCA  now 
filling  this  role,  along  with  periodic  baseline  concept  updates  will  be  the  basis  of 
several  key  program  decisions  in  the  next  few  months  and  years.)  The  second  is 
to  accurately  calibrate  the  CFD  tools  and  processes  that  are  being  developed  to 
enable  the  design  and  optimization  the  future  production  HSCT. 

From  the  preceding  discussion  of  the  TCA  drag  polar  build-up  it  can  be  seen 
that  the  high  speed  drag  polars  for  performance  modeling  are  a complex 
combination  of  increments  and  absolutes  from  different  sources.  Each  of  these 
inputs  has  its  own  inherent  precision  and  bias  errors  — some  easier  to  determine 
and  control  than  others.  For  the  purposes  of  the  “flight  polar  build-up”  ADP  sub- 
task, data  used  in  the  TCA  polar  build-up  was  assessed  for  potential  sources  of 
bias  and  precision  (repeatability)  error. 
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• Test-to-test,  and  tunnel-to-tunnel  differences  (1 .5ct) 

• Wind  tunnel  model  corrections: 

- Trip  drag  / laminar  run  (2.0ct) 

- Model  fidelity  and  aeroelastics  (2.5ct) 

- Nacelle  internal  forces  (2.0ct) 

- Installation  T&l  (~2ct  ? trans.  / 0 super.) 

• “Full  scale  airplane”  drag  increments: 

- Excrescence  (2ct) 

- Empennage  / trim  (6ct) 

- Propulsion  induced  effects  (lOct  trans.  / 2 super.) 

- Reynolds  number/skin  friction  effects  (4ct) 


Significant  sources  of  uncertainties  in  the  absolute  drag  levels  quoted  for 
performance  use  were  identified  based  on  existing  discrepancies  between 
various  semi-empirical,  CFD,  and  wind  tunnel  sources.  Uncertainties  in  the 
ability  to  obtain  improvements  with  non-linear  optimization  technology,  and  any 
biases  related  to  the  basic  linearized  potential  flow  design  and  analysis  tools 
were  not  considered. 

Depending  on  the  exact  data  sources  compared,  varying  uncertainty  levels  can 
be  assessed.  Those  shown  in  the  chart  above  are  believed  to  be  the  most 
representative  relative  to  the  specific  components  in  the  TCA  polars.  All 
quantities  are  +/-  the  value  shown. 

As  part  of  BCAG’s  portion  of  the  ADP  task  in  1996,  the  areas  of  excrescence 
drag,  empennage  drag  (upsweep,  closure,  tail  trim),  and  wind  tunnel  database 
corrections  and  averaging  were  chosen  for  more  detailed  investigation.  CFD 
solutions  (Overflow  and  TLNS3D-MB)  were  used  to  get  a better  assessment  of 
the  the  empennage  drag  terms  and  Reynolds  number  effects.  It  should  be  noted 
that  the  wind  tunnel  database  uncertainties  generally  get  incorporated  into  the 
status  polars  both  directly  through  the  “test-theory  increments”,  and  indirectly 
through  the  CFD-based  geometry  increments  as  the  CFD  codes  were  validated 
using  wind  tunnel  data. 

(The  potential  uses  of  CFD  data  in  the  assessment  of  excrescence  drag  and 
wind  tunnel  trip  drag  are  the  subject  of  additional  BCAG  papers.) 
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TCA  Drag  Build-Up  “Test-Theory  Increment” 


The  tunnel-to-tunnel  repeatability  as  seen  on  the  Reference-H  configuration 
(roughly  +/- 1 .5  count)  has  been  a concern  for  several  years.  The  “Benchmark 
Models  Program”  co-op  testing  program  begun  by  BCAG  and  NASA  prior  to 
HSR-2  was  intended  to  uncover  and  resolve  such  discrepancies  to  reduce 
uncertainty  in  absolute  level  drags.  (The  Ames  Unitary  tunnel  was  unavailable  to 
complete  the  model-to-model  portion  of  the  program  and  the  remainder  of  the 
study  was  out-prioritized  by  HSR-2  and  BCAG  activities.)  ‘Test-theory 
increments”  to  be  applied  to  linearized  potential  flow  results  are  shown  in  the 
above  left  bar  chart  based  on  various  Ref.-H  wind  tunnel  sources.  Error  bands 
about  the  individual  wind  tunnel  levels  are  shown  representing  80%  confidence 
levels  for  measurement  repeatability  and  uncertainties  in  the  data  corrections 
(primarily  trip  drags  and  model  fidelity).  Using  just  the  Ref.-H  wind  tunnel 
database  would  have  indicated  that  the  TCA  drag  prediction  would  need  an 
average  1 .5  count  test-theory  increment  at  Mach  2.4.  As  shown,  this  was 
confirmed  by  the  Ref.-H  Overflow  Navier-Stokes  results.  Pre-test  CFD  results  for 
the  TCA  however,  indicated  a level  about  1 count  below  the  uncorrected  potential 
flow  results.  As  it  was  believed  that  the  higher  sweep  angles  and  more  slender 
body  of  the  TCA  might  have  more  favorable  “real  flow”  effects  than  the  Ref.-H,  a 
test-theory  increment  of  -1  count  was  chosen  for  the  TCA  performance  polars. 

The  Ref.-H  wind  tunnel  database  test-theory  increments  was  still  used  as  a guide 
in  fairing  out  the  TCA’s  -Icount  across  the  Mach  range  (see  upper  right  graph). 
TCA  transonic/subsonic  test  results  from  LaRC  16T  are  still  being  analyzed. 
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Initial  TCA  test  results  show  that  the  -1  count  test-theory  increment  ag 
the  measured  drag  polar  of  the  “clean  wing”  model  if  one  assumes  that  the 
adjustment  for  trip  drag+laminar  run  should  correctly  be  2.5  to  3.0  counts. 

Unfortunately  the  statistical  uncertainty  in  trip  drag  correction  due  to  curve  fit 
questions,  basic  data  scatter,  and  laminar  run  correction  unknowns,  is 
siqnificantly  larger  than  the  test-theory  correction  so  a definitive  answer  is 
impossible  at  this  time.  A method  using  CFD  results  and  excrescence  drag 
calculation  methods  to  reduce  these  correction  uncertainties  shows  promise  but 
requires  calibration  with  additional  very  high  quality  wind  tunnel  data. 

The  other  fact  that  can  be  gathered  from  the  above  test  data  plot  is  that  the 
impact  of  model  fidelity  (including  model  excrescence)  cannot  be  ignored ■ 
the  test-to-test  repeatability  of  the  TCA  wing  which  included  flap  cut-outs  was 
excellent  (fractions  of  a count),  it  clearly  has  a 1 .5  count  model  fidelity  penaty 
relative  to  the  second  model  which  was  built  with  a “clean  wing, 
fidelity  penalty  grows  to  nearly  2 counts  for  the  nacelles-on  case. 


Subsonic  and  transonic  TCA  force  data  obtained  recently  in  the  URC  16R  tra 
tunnel  are  still  being  analized  at  the  time  of  this  wntmg_  No  ofr:d®s'9h  supf[? , 

(i.e.  Machs  front  1 .3  - 2.1)  were  taken  on  the  TCA  performance  model  so  the  large 
test-theory  increment  predicted  from  the  Ref.-H  data  cannot  be  confirmed. 
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Several  areas  were  identified  where  viscous  CFD  analysis  might  provide  additional 
data  to  reduce  the  uncertainty  levels  in  HSCT  performance  polars.  The  use  of 
Navier-Stokes  results  in  the  resolution  of  trip  drag  corrections  and  excrescence 
drag  estimates  has  already  been  mentioned.  Navier-Stokes  results  were  also 
used  to  provide  an  assessment  of  empennage  drag  and  possible  Reynolds 
number  effects  not  capytured  by  the  simplified  PD  methods  used  in  the  current 
polar  buildups. 

Considerable  time  and  effort  went  into  using  TLNS3D-MB  for  empennage  analysis. 
Seven  different  aft-body/  vertical  tail  geometries  were  run,  plus  three  tail  angles  on 
the  baseline  TCA  body.  While  the  results  of  the  trim  drag  portion  were 
encouraging,  the  aft-body  upsweep  and  closure  increment  study  was  inconclusive. 
Grid  problems  with  viscous  flow  continuity  at  block  boundaries  and  the  unexpected 
complexity  of  the  aft-body  viscous  flow  caused  significant  errors  in  both  the  skin 
friction  and  viscous  pressure  drag  of  the  various  bodies.  This  problem  is  currently 
being  investigated  further. 
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CFD  Assessment  of  Empennage  Drag  and  Rn  Effects 


• TLNS3D-MB  study  of  empennage  drag: 

- 7 aft-body  / tail  geometries 

- 3 tail  angles  on  baseline  TCA  aft-body 

- Horiz.  is  unported  at  best  trim  angle 

- Status  trim  o.k.  (+1ct  downwash,  -let  vert-horiz.  interf.) 

- Upsweep  / closure  drag  results  inconclusive 

• Ref.  H and  TCA  wing-body  @ w.t.  and  flight  Rn 

- Overflow  (dense  grid,  Bladwin-Barth) 

- Both  showed  dCD/dRn  slope  2ct  less  than  flat  plate  calc. 

- -let  vise,  pressure  drag  at  flight  on  Ref.-H 

- +2%  vise.  CL  at  flight  on  TCA 


The  trim  drag  results  showed  that  the  baseline  tail  is  unported  at  the  optimum  tail 
angle.  If  the  tail  is  “resealed”  to  the  body  wiping  surface,  the  large  effective 
downwash  (4  degrees)  and  poor  installed  tail  polar  cause  a one  count  trim 
penalty.  The  total  status  tail  drag  is  about  equal  to  that  predicted  using  the 
normalized  trim  approximation  as  there  is  also  about  one  count  favorable 
horizontal-vertical  tail  volume  interference.  Similar  results  were  obtained  across 
the  Mach  range.  The  tail  polar  and  sensitivity  to  downwash  and  body  shape 
indicate  that  the  projected  additional  reduction  of  1 to  2 counts  is  probably 
achievable  with  non-linear  optimization.  Of  course  wind  tunnel  data  is  still 
needed  to  validate  the  empennage  CFD  results  for  future  designs  and 
predictions. 

Overflow  solutions  on  the  TCA  and  Ref.-H  wing-bodies  were  studied  to  identify 
trends  in  skin  friction  and  viscous  pressure  drag  at  cruise  Mach  conditions.  Both 
geometries  showed  a CDf  trend  from  wind  tunnel  to  flight  Reynolds  number  that 
is  about  2 counts  less  than  predicted  from  standard  fait  plate  skin  friction 
methods.  In  addition,  the  Ref.-H  showed  a potential  1 count  pressure  drag 
reduction  in  flight.  The  TCA  showed  a potential  2%  increase  in  lift  at  flight 
Reynolds  number.  Again,  test  data  and  additional  CFD-to-flat-plate  theory 
comparisons  are  needed  to  reduce  uncertainty  in  this  area. 
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Summary  Of  Drag  Polar  Build-Up  For  High  Speed  Performance 


• “Performance  Polars”  = Projected  Level  from  tops-down  method 

“Tech.  Projection”  = (Projected  Level)  - (Status  polars  @ flight) 
“Status”  =(W.T.  database)  + (potential  flow,  PD,  CFD-based  A’s) 

• Uncertainty  bands  around  Status  can  add  significant  program  risk 
(e.g.  viability  assessment,  technology  metrics,  downselects) 


In  addition  to  CFD  analysis  and  optimization  tools,  HSCT  testing 
and  accounting  methods  must  be  validated  to  a high  confidence 
level  before  production  go-ahead,  e.g — 

- Empennage  drag 

- Propulsion  induced  effects  (PIE) 

- Rn  effects  on  supersonic  pressure  drag  and  skin  friction 


This  paper  has  reviewed  the  various  components  of  the  TCA  high  speed  drag 
polars  used  for  airplane  sizing,  mission  performance  analysis,  and  Technical 
Integration  (Tl)  airplane  level  trade  studies.  The  level  of  cruise  L/D  in  the  TCA 
polars  has  been  set  equal  to  the  target  cruise  performance  determined  by 
concensus  of  the  joint  NASA-industry  Configuration  Aerodynamics  (CA)  team. 

An  unprojected  "status"  level  TCA  polar  was  developed  based  on  a combination 
of  available  wind  tunnel  data,  geometry  changes,  and  "real  airplane"  details.  The 
"technology  projection"  is  the  drag  improvement  required  to  cover  the  difference 
between  the  concensus  "tops-down"  based  performance  L/D  and  the  status 
polar.  A projection  formula  based  on  the  percent  change  in  cruise  wave  drag, 
and  drag  due  to  lift  was  then  applied  at  all  supersonic  Mach  numbers.  A 20% 
excrescence  reduction  at  all  Machs  and  a 6%  improvement  at  subsonic 
conditions  are  also  applied. 

Potential  errors  in  the  absolute  level  of  predicted  drag  can  have  an  adverse 
impact  on  the  determination  of  program  viability  or  major  concept  downselects. 
Several  sources  of  uncertainty  in  the  current  polars  have  been  described.  Some 
efforts  under  the  Tl  "ADP"  task  to  further  quantify  or  reduce  several  uncertainty 
elements  have  been  shown.  The  use  of  calibrated  CFD  solutions  appears  to 
offer  benefits  in  several  areas.  Future  wind  tunnel  tests  to  extend  the  existing 
wind  tunnel  databases  and  provide  CFD  validation  must  be  carefully  planned  to 
obtain  maximum  data  quality  if  significant  progress  is  to  be  made  in  eliminating 
uncertainty  sources. 
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Task  2.1  - Technology  Integration 

Sub-task  2.1.5  Aircraft  Design  Processes 
Sub-Sub-task  2. 1.5.1  High  Speed  Aerodynamic  Processes 
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Mach  1.1  Analysis  of  BRISTOW  Designs 
CD  Induced  Summary 
Preliminary  Conclusions  and  Future  Plans 
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ACD  Trim  -1  Count  Assumption  tbd-time 

ACD  Techno.  f(CD,LE  Break)  Assumption  permitting 


For  Technology  Integration  studies,  all  methods  will  be  corrected  with  a test  to 
theory  increment,  calculated  from  a common  database  developed  by  Boeing’s 
parallel  task.  Therefore,  getting  the  correct  change  in  drag  between  designs  is 


CD  >> 

X ~ 


T3 

C <D 

a 7D 
G CS 

<D  ^ 
"O  <J5 

a c 3 
3 £ 

o 60 

H c 

13 
— • -o 

<D  ■ O 

j>  £ 

CD 

60  X 
P 

" C/3 

73 

o O 
G O 

3 a 

« SP 

X i-H 


<D  C T3 
60  — O 

_2  - x 

O D 

e 

<D  u_ 
^ "O  n 


c 0)  CL, 

•o  c 2 

U ts  t 
a a £ 
G-l  O 

O O on 


g X 

73  ^ 

p T3 
h C 


60-0 
G o 

* »H  - \ 


c 

cs  c 

X CD 
*->  <o 

£ 

c3  P 

g 

O on 
Oh  CD 

§ § 
CD 

Uh 

2 2 
O G 

S 73 


G ft  3 
O g C 

4-^  4-J 

o e 

a -a  o 

C -G  <D 

o *■* 

O 4-»  5/5 

d O <- 
o X!  CO 

*g  G < 

a £ z 

c 

G o 

ffl  G . 

2 73 
• CL  CD 

s«  00 

2 •*-*  03 

T3  cd  (D 
O C g 
0^-0 


T3  C D 
O *-H  X 

0-0  “ 

■§  % <=? 
> G .g 

Jg  ^-o 
- 60  g 

G O Ch 

--h 

w *§  Q 
® o < 

c *5 

O O ►— 1 

42  S g 

L-c  (D 

* C G 

U g c 

-£  D 

X>  c3 
^ P 

fe  o cn 


. W C 

o GO  P 

B5  3 ^ 

— t C .<D 


G ’G 


UJF 

13  ^ 

"§  o 

on  *73 
p (G 


=3  fa 

73  £ £ 

o ^ 

O CD  - 

60  60  G 

2 *2  *o 

3 w)  § 

.S  « 

c g g 

<D  03  G 
<D  <D  g 

X ^ O 

a .s  73 

X C/3 

C/3  £ G 

w 33 
<"*!  C ^ 

5 U C/3 

SEW 

(U  ID  Q 

> O M 

0 c u 

1 d 7 

o-  x a 


OJ  C/3 

X 3 
*■*  O 

x '> 

a ^ 

§ a 

S OD 

2 £ 

"O  c 

CD  O 

0/3 

> ^ 
53 

^ cd 
£ X) 

P C/3 

£ ~ 
k.  bO 

^ G 

• fH 

O .G 
H g 
m G 


PQ  § 

c £ 

■ '“H  w 

C^3  /IV 

(D  ON 
C 00 
•—  cn 

o ^ 

G 00 
t>0  ^ 

2 q 

s 


. <D 

a £ 

g g 

G <D 

<D  u_i 


<D  "O  C 
crt  o T3 

03  C O 


8 » n .§  S ^ 


c3  X 
c 


cb  X 

-4— * 

5 <D 


c3  <D 


bO  <D 

•o  ^ 

T3  c 
<D  G 

•*— > Uh 

o o 

G c 
£ | 
CL 

£ X 


(D  CD 

g oc 

cd 

w T3 
Vh  c 

O § 

CL 

E S 

C/3  CLH 

CD  G 
G G 

C3  h 

P Oh 
bO  ^ 

C4-4  o 

o fc; 

2 < 
68  l- 

00  J. 

-o  G; 

C <N 


(D  C 
O <D 

C 

’33  ^ 

4—> 

^ CD 

>,  x 


X CD 

a <4_h 
*0 

< T3 


733 


compared. 


high  sna  hisuhch 


C/3 

G 

60 

• ^ 

C/3 

<D 

73 

g 

<d 

<D 

£ 

15 

X) 

oo 

Q 

55 

W 

Pi 

H 

+-> 

o 

<D 

fc 

o 

o 

C/3 

• 1— I 

cd 

O 

60 

v-H 

=tt= 

o 

X 

H 


c^> 

CD 

C/3 

cd 

o 

4— > 

C/3 

CD 

Tin 

s 

• PH 

c/a 

Q 

JS 


2 

<4-1 

T3 

CD 

w 

Uh 

cd 

+-> 

c/3 

<D 

& 

c/a 

C 

O 

c/a 

‘C 

cd 

Oh 


O 

U 


C/3 

73 

G 

<D 

Vh 

H 

>> 

73 

O 

CO 

73 

a 

cd 

60 

G 


cd 

E 


00 

c 

*5  C 
3 D 
3 D 

O £ 
C £ 

a)  -D 

G3  O 

t:  ti 

>*  o 
o _, 

C c 

S o 

+Z  •!-! 

c« 


C/3 

c 

o 

o 

c 


J3 

o 

& 

d 


—<  c 


73 

D 

o 

s 

73 

H 

aS 

J3 

■*— * 

<L> 

TD 

cd 

e 

CD 

u 

D 


c/a 

CD 

W) 

c 

cd 

pC 

o 

00 

c 


73 

O 

a 

73 

c 

< 


73 

(D 

O 

S 

73 

CD 

ip- 

73 

C 

cd 

'd 

CD 

C 

CD 

i-i 


,D 
4-h 
i— i 

o 

w 

CD 

00 

G 

•i-N 

C 

•PH 

73 

CD 

<D 

W) 

CD 

iS 

O 

13 

o 

Vh 

C/3 

#s 

a. 

<+H 

<D 

£ 

i-H 

£ 

o 

H 

CD 

X> 

00 

h-H 

c/a 

<D 

C* 

73 

o 

CQ 

D 

II 


734 


Wing  Design 

Wing/Warped-Plate-Body  Analysis  of  Best  Designs 


Before  this  task,  the  drag  of  a design  could  vary  2 to  7 counts  at  Mach  2.4  and 
was  different  by  as  much  as  10  counts  at  Mach  1.1.  Modeling  and 
interpolation  processess  were  made  more  consistent,  and  paneling  density  was 
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FLAT  WING  AND  BODY  - NO  THRUST 


736 


WINGDES  SDAS/A389  BRISTOW 


All  three  codes  use  the  same  method  to  calculate  full  leading  edge  thrust,  so 
the  differences  in  thrust  increments  are  due  to  differences  in  the  pressures  used 
to  extrapolate  the  leading  edge  pressure.  The  consistent  drag  increments, 
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WINGDES  SDAS/A389  BRISTOW 


Attainable  thrust  increments  also  agree  within  1 count,  but  notice  that  the 
SDAS/A389  attainable  thrust  is  less  limiting-relative  to  WINGDES  and 
BRISTOW.  Probably  due  to  differences  in  the  reference  7 versus  reference  3 
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FLAT  WING  AND  BODY  - ATTAINABLE  THRUST 
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WINGDES  SDAS/A389  BRISTOW 


NASA,  Boeing,  and  MDC  each  use  different  codes  that  calculate  linearized 
potential  flow  using  a vortex-lattice.  Even  though  all  codes  employ  the  same 
theory,  each  code  uses  a different  method  to  design  optimally  cambered  wings. 
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CD  Induced  Codes 
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In  FY95,  a great  deal  of  effort  was  expended  to  understand  what  the  induced 
drag  increment  should  be  between  consistent,  well  optimized  Arrow  and  Ref. 
H planform  designs.  Wind  tunnel  data  of  CFD  optimized  designs  were 
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Drag  levels  are  shown  below  for  designs  optimized  and  analyzed  by  their  own 
code.  Since  WINGDES  is  the  only  code  that  can  design  for  Leading  Edge 
Thrust,  two  sets  of  WINGDES  designs  are  shown:  with  LET  and  without 
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MACH  2.4  DESIGN  OPTIMIZATION  COMPARISON  - NO  THRUST 
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WINGDES  WINGDES  SDAS/A389  BRISTOW 

DESIGN  @ DESIGN  DESIGN  @ DESIGN  @ 

CL=0.10  withLET  @ CL=0.10  CL=0.10 

CL=0.10 


MACH  2.4  DESIGN  OPTIMIZATION  COMPARISON  - FULL  THRUST 
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WINGDES  WINGDES  SDAS/A389  BRISTOW 

DESIGN  @ DESIGN  DESIGN  @ DESIGN  @ 

CL=0.10  withLET  @ CL=0.10  CL=0.10 

CL=0.10 


MACH  2.4  DESIGN  OPTIMIZATION  - ATTAINABLE  THRUST 
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WINGDES  WINGDES  SDAS/A389  BRISTOW 

DESIGN  @ DESIGN  DESIGN  @ DESIGN  @ 

CL=0.10  withLET  @ CL=0.10  CL=0.10 

CL=0.1 0 


WINGDES,  SDAS/A389,  and  the  BRISTOW  design  trends  are  within  1 
count  of  each  other  and  are  within  1 count  of  the  established  trend  when  the 
CD  induced  tech,  projection  is  included  (tech,  projection  = -7.3  counts  for 
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MACH  2.4  DESIGN  TREND  COMPARISON  - ATTAINABLE  THRUST 
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WINGDES  WINGDES  with  LET  SDAS/A389  BRISTOW 


Qualitative  Assessment: 
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the  LET,  nacelle  and  fuselage/thickness  effects  unique  to  each  code,  will  be 
investigated  in  future  studies. 


Ref.  H Planform  Optimization 
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SPANWISE  LOCATION  Y,  FT 


Of  course,  these  trends  were  based  on  each  code  analyzing  its  own  designs. 
To  ensure  consistency,  the  lowest  drag  BRISTOW  designs  were  choosen  to 
be  analyzed  by  the  other  codes  to  see  if  their  results  would  match  the 
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BRISTOW  DESIGN:  ANALYSES  COMPARED  - NO  THRUST 

Mach  2.4,  CL=0.10 
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WINGDFS  SDAS/A389  BRISTOW  DESIGN  @ CL=0.10 


BRISTOW  DESIGN:  ANALYSES  COMPARED  - ATTAINABLE  THRUST 
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WINGDES  SDAS/A389  BRISTOW  DESIGN  @ CL=0. 10 


Predictions  compare  well  with  and  without  thrust  estimates  at  CL=0.10; 
however,  at  the  minimum  drag  CL,  SDAS/A389  predicts  several  counts  more 
thrust-resulting  in  negative  drag  (at  Mach  2.4  and  1.1).  Further,  SDAS/A389 
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BRISTOW  DESIGN:  ANALYSES  TREND  - ATTAINABLE  THRUST 


S 

1 


LU 

Q 

0 


§ 


CO  LO 

(sjunoo)  0000 1 


Tf  CO  CM 

paonpuiQO  V113Q  V - a N\Tld 


o 


ii 


758 


Mach  2.4  Analysis  of  BRISTOW  Designs  - Full  Thrust 
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TCA6  LET  Prediction  Comparison  at  Mach  2.4 
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Here  SDAS  is  shown  without  its  LET  estimate,  since  that  is  the  way  Boeing 
uses  the  code  at  Mach  2.4.  Interestingly,  the  trend  is  a better— much  flatter- 
comparison  versus  CL.  (The  trend  mimics  the  BRISTOW  fuselage  and 
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TCA6:  METHOD  minus  TEST  at  Mach  2.4  ■ ADRAG  VERSUS  CL 

Linear  CDwave  = 0.00187 


Poly.  (CFL3D  Euler) 

Poly.  (BRIS) 

-0.0005  0 Poly.  (WDES) 

Poly.  (BRIS(Fuselage  & Thickness)) 


HIGH  SPEED  RESEARCH 
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BRISTOW  Designs:  Mach  1.1  CDinduced  at  CL=0.15 
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WDES  SDAS  BRIS  WDES  SDAS  BRIS 

Arrow  Planform  Ref  H Planform 


A plot  of  the  spanwise  leading  edge  thrust  prediction  from  WINGDES  and 
BRISTOW  were  compared  to  examine  why  thrust  predictions  differed. 
WINGDES  had  to  reduce  is  paneling  to  19  columns  at  Mach  1.1.  BRISTOW 
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SPANWISE  THRUST  DISTRIBUTION-ARROW  at  MACH  1.1 
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